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Summary 

A further survey of southern hemisphere meteor activity, extending from 
1956 December to 1958 August, has been made with the 67 Mc/s narrow-beam 
equipment. ‘The faintest detectable meteor trail has an electron line density 
of 1:3 10'' electrons/em. The potentialities of radio equipments for the 
detection of weak meteor showers are demonstrated by an objective statistical 
analysis of echo counts. Weaker showers can be detected by the present 
equipment only if they transit within 30 of the zenith. .The statistical 
analysis reveals an inherent clustering in the ecliptical component of sporadic 
meteors, and there is some evidence for very weak or short-lived showers 
which do not satisfy all the criteria laid down for an accepted shower. Radiants 
and echo rates for fifteen accepted and seven possible showers are tabulated 
and discussed. Echo rates measured at different equipment sensitivities for 
four showers indicate in all cases a deficiency, relative to sporadic meteors, in 
the numbers of fainter meteors in these streams. 





1. Introduction.—This paper deals with the meteor showers detected during 
a radio survey of meteor activity at Adelaide. With the exception of a few short 
breaks, the survey was continuous from 1956 December to 1957 December 
inclusive, and again from 1958 April to August. No records were obtained 
during 1958 January-March. ‘These observations continue the survey of 
Southern Hemisphere meteor activity which has proceeded intermittently since 
mid-1952 (Weiss (1, 2) hereafter referred to as Papers | and I1). 

The potentialities of radio equipments in the identification and measurement 
of weak showers is a subject which has received little attention in the literature. 
For this reason a considerable part of this paper is devoted to a discussion of the 
factors which determine whether or not a given shower is detectable with the 
narrow-beam equipment in use, and to a description of the method developed 
to minimize the subjective element in the identification of the weaker showers. 
This is followed by factual data on all showers identified during the survey. 

The application of an objective method of analysis of data, collected with an 
equipment whose inherent resolving power is superior to that of the previously used 
27 Mc/s wind equipment, has resulted in the recognition of several new showers, 
which were either missed or regarded as uncertain in previous surveys. 

2. Equipment.—The 67 Mc/s radiant equipment is basically the same as 
that used in 1953 and described in Paper I, except that the peak pulse power 
has been increased to 50kW. ‘The discrimination afforded between noise 
impulses and meteor echoes by transmission of double pulses 100-1 sec apart 


29 
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was supplemented in 1957 by a jitter on the time-base of the recording tube, 
which resulted in the quadrupling of each echo. ‘This however was found to be 
unnecessary, and ample discrimination was obtained in 1958 by retaining the 
time-base jitter and transmitting a single pulse 25-sec wide. The time of 
occurrence and slant range of each echo were recorded. 

The axes of the two identical beamed aerial arrays are directed at elevation 
g° along azimuths 35°-5 S of E (N aerial) and 14° N of E (S aerial). Half-power 
beam widths are +3°°5 in azimuth and +4°-5 in elevation. The collecting 
sectors of these two aerials, which contain the radiants of all the detectable 
meteors, are narrow sectors stretching from horizon to horizon along great 
circles perpendicular to the beam axes. The daily rotation of the Earth sweeps 
these collecting sectors across that part of the celestial sphere which is visible 
from Adelaide, so that each portion of the visible celestial sphere is sampled 
twice per 24 hours, once by each sector. Each collecting sector crosses a point 
source in from 1 to 2 hours, except when the source lies close to the south celestial 
pole, where the times of passage are much longer. The time-difference between 
reception of echoes from a given source in the two aerials is a unique function 
of the declination of the source. Curves for the determination of meteor shower 
radiants from measured times of reception of echoes from the two sectors are 
given in Paper I. ‘The finite beam widths of the aerials result in a variation of 
echo range with time, as the collecting sectors pass across the source; typical 
range-time envelopes have already been published (3). 

3. Calibration of the equipment.—The average sporadic echo rate is about 
700 per day in each channel. 

The electron line density in the faintest detectable meteor trail has 
been estimated by two indirect methods. A limiting line density of 
1:2 +0°2 x 10! electrons/cm has been deduced for 1957 October to December 
from a comparison of echo rates at different equipment sensitivities (4). The 
second method, based on shower meteor fluxes, will be described in detail. 

Following the procedure already applied to the 27 Mc/s wind equipment (5), 
the echo rate N, for the radiant equipment has been determined as a function 
of the meteor particle flux F, at those times when the direction of maximum 
equipment response lies in the echo plane, i.e. when the direction of flight of the 
meteor particles proceeding from an individual radiant is normal to the direction 
of maximum response. ‘This direction lies in the azimuth of the beam axis 
and at elevation 10-5; the slight displacement from the direction of the axis 
of the aerial beam is due to the variation of returned echo power with slant range 
of the reflection point. Assuming an inverse power law for the differential 
distribution of meteor masses, the integral flux function F becomes* 


F(2,) = (c/(s—1))2,!-* (1) 


where «, is the maximum electron density in the trail at the limit of detection in 
the direction of maximum response. The ratio N,/F (km)? is plotted in Fig. 1 
as a function of radiant zenith angle x for three values of s._ If the flux F can be 
measured independently, x, can be found at once from (1). The Geminid 
meteor stream is most suitable for this calibration; the maximum echo rate is 
remarkably constant from year to year, and the mass distribution is well 


* It is assumed that the ionizing efficiency of a meteor is independent of its mass. 
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Fic. 1.—The echo rate N for the radiant equipment, as a function of the zenith angle of the 
radiant at the time when the direction of the meteor flight is normal to the direction of maximum 
equipment response. F is the integral meteor flux, and s is the mass distribution parameter (see 
equation (1)). The total echo rate (broken curve) is the relative number of echoes received whilst 
the radiant lies within the collecting sector of the aerial; it is drawn for sporadic meteors, s=2:0. 


determined (s~1-5 in the brightness range which contributes most of 
the echoes). Now the 27 Mc/s wind equipment has been calibrated 
absolutely, and equation 6(a) of reference (§) gives the relation between N,, 
and F(x,,) for this equipment when x=70 and «,.=3°5 x 10" electrons/cm. 


With s=1°5, 
2/N.\2 
%, =3°5 x 104 (5) (=) electrons/cm, (2) 


N,,=24/hr (Paper Il). The ratio f= N,/F(«,) is to be taken from Fig. 1 at the 
zenith angle of the Geminid radiant at the time of detection at maximum response, 
and N, is the instantaneous echo rate at this time, which may be ascertained 
from the range-time plots. For the N aerial, y=75 , f=4050km*, N,=54/hr 
(1957), leading to 2,=1-21 x10" electrons/cm. For the 5S aerial, y=70°, 
f=5230km’, N,=66/hr, «,=1°35 x10! electronsicm. These results are 
insensitive to a small error in the value of s, and to the finer details of the aerial 
polar diagrams. The calibration, of course, applies only to 1957 December, 
but the shower echo rates for 1956 and 1957 were almost identical. 

A mean value of «,=1-3 x 10!! electrons/cm agrees satisfactorily with the 
above-quoted independent estimate of 1-2x10'!. The fact that this latter 
calibration applies to a different period lends considerable weight to the 
supposition that the equipment sensitivity was essentially constant over the 
whole of the survey. 

29* 
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The limiting line density «,=1-3 x 10" electrons/cm is the line density of 
the faintest trail segment detected by the equipment; it corresponds to a meteor 
some three magnitudes below the visual threshold. It is not at the same time the 
maximum line density in the trail which is detected most frequently. Because 
the effective aperture of the aerial increases with the line density at the reflection 
point, then in the range of maximum line densities from amax tO %max + 4%max 
the equipment collects a maximum number of trails when %max~22,. This 
result may be derived as follows. Again assuming the elementary integral 
distribution (1), the actual differential distribution of the maximum line densities 
of all trails detected by the equipment, regardless of the position of the reflection 
point on the trail, may be written as 


N damax = W(amax)%max~* dxmax. (3) 


The weighting function W is proportional to the volume of space within which 
a trail, whose maximum line density lies in the range amax to &max + d%max, May 
be detected. It depends on the polar diagram of the aerial, and has been 
calculated by an extension of the numerical approach adopted by Clegg (6) 
in determining the collecting area for a point radiant. The extension of the 
weighting factor W from a surface to a volume has been made because with the 
fixed aerial system in use the collecting area appropriate to a point radiant moves 
across the whole of the aperture of the aerial as the radiant is carried across the 
collecting sector by the motion of the earth. In other words, the characteristics 
of shower echoes are determined as an average over the whole of the range-time 
envelope. Differential distributions for the echo rate for three values of s are 
given in Fig. 2, along with the curve for a meteor distribution of the form 


N damax = Wamax''5+0"453 108 (a asl? damax 
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Fic. 2.—Differential distributions of the numbers of meteor trails with maximum line densities 
in the range maz tO %maz +d%mar, as detected by the radiant equipment. Full lines are for constant 
values of the mass distribution parameter s. The broken line is for s = 1-5 +0°453 log (%maz/ lim): 
Xtig 18 the equipment threshold. 
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which, apart from the particular numerical values, is a possible representation 
of the mass distribution data for such showers as the Geminids and the Perseids. 
It is seen that the location of the maximum in the echo rate distribution is not 
strongly dependent on the mass distribution of the incident meteors, and lies 
about 1 magnitude above the equipment threshold. 

The computations leading to Fig. 2 apply strictly only to a radiant which 
transits at the zenith. However, collecting areas examined during the preparation 
of Fig. 1 were found to be quite insensitive to the zenith angle of the radiant, 
and it is expected that the curves of Fig. 2 will apply with good accuracy to all 
radiants other than those which transit at low elevation. 

The polar diagrams of most conventional aerial systems, if reduced to constant 
half-power widths, are very similar in the vicinity of the beam axes. The above 
conclusion, that the position of the maximum in the differential echo rate 
distribution as a function of radio brightness lies about 1 magnitude above the 
equipment threshold, should therefore not be limited to the Adelaide radiant 
equipment. It is a general result which will apply to the majority of radio 
equipments used in meteor studies. 

4. The recognition of shower activity.—A detailed discussion of what is 
understood by the term ‘‘ meteor shower ’’ has already been given elsewhere (3). 
This question will not be pursued here, as the context of this paper should make 
clear the significance attached to this term. 

Shower activity is detected and measured against a background of sporadic 
meteor activity, whose fluctuations set the ultimate limit to recognition of shower 
activity. ‘The background activity has strong diurnal and seasonal variations, 
on which are superimposed the essentially random fluctuations which are 
examined in detail in Section 5. The sporadic rate ranges from a maximum 
of 100 echoes/hr near mid-day in March to a minimum of 2 echoes/hr near 
sunset in September. The corresponding range in the absolute value of the 
statistical fluctuations. in the sporadic echo count precludes any constant 
threshold for shower detection; the chances of discovery are heavily weighted 
in favour of radiants with longitudes within + 60° of the antapex. Most showers 
in fact are detected from 22 hr through midnight to 10 hr L.T., when the sporadic 
rate is also high (average 40/hr), and the limiting echo rate at which a shower 
is detectable is 20~-30/hr. 

The method of searching the records for shower activity is based on an 
objective statistical analysis of the echo counts (3). The problem, as already 
mentioned, is one of disentangling rate increases due to shower activity from 
random fluctuations in the background count. Briefly, the echoes detected in 
each half-hourly period are counted, and the probability of obtaining any 
particular count is ascertained on the assumption that the background count 
is random and follows Poisson’s distribution. The validity of Poisson’s 
distribution has been established already for periods up to 10 min duration (3). 
During the present survey additional evidence, not described here, has been 
obtained that this distribution is valid also for echo counts in 30 min periods. 

Periods of high count, which may include shower activity, are further 
scrutinized by plotting individual echoes in a range-time plot. Only those 
periods giving recognizable range-time envelopes are accepted as certainly 
including shower activity, unless the existence of the shower has been 
demonstrated independently by other observers. This reservation is made 
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because the clarity of a range-time envelope depends on the diffuseness of the 
radiant, and is also severely prejudiced by intermittent activity, which seems 
to be characteristic of some of the weaker showers, although this patchiness 
can be overcome to some extent by integration of the records for successive days. 

In addition to considerations of meteor particle flux, time of detection, and 
radiant structure, the probability of detection of a shower is influenced strongly 
by the declination of the radiant, through a dependence of echo rate on zenith 
angle of the radiant at the time of passage through the collecting sector. Fig. 1 
illustrates this dependence for those special instants when the direction of 
maximum equipment response lies in the echo plane. This is not the most 
useful representation, as the total number of echoes received whilst the radiant 
lies inside the collecting sector depends also on the shape of the range-time 
envelope and the duration of the passage, both of which are determined by the 
declination of the radiant. ‘The concentrated Sun and anti-Sun sources of 
sporadic meteors (7, 8) afford the most reliable indication of this variation of 
total echo rate with zenith angle of the radiant. The mean total echo rate from 
these sources, as found in the present survey, is sketched in Fig. 1. Since this 
refers to sporadic meteors, the value of s is probably close to 2-0, but there is 
no reason to doubt that a similar strong variation of total echo rate would apply 
also to showers whose s-values may differ from 2:0. 

It is therefore to be expected that the probability of detection of a shower 
will increase as the elevation of the radiant at transit increases from the northern 
horizon to the zenith. Beyond this to the south, the probability of detection will 
fall off rapidly, not only because of decreasing echo rates for a given flux, but also 
because the rapidly increasing duration of the passage of the radiant across the 
collecting sector may result in more, rather than less, confusion with the 
background count. 

The tendency for the radiants of the weaker showers of 'l'able III and of the 
possible showers of ‘Table VII to group in the vicinity of 5=—35 is a clear 
indication that the radiant type of equipment is most sensitive to showers whose 
radiant declination is approximately equal to the latitude of the observing station. 
A similar, spurious grouping of declinations near 6= —35 could arise from a 
subjective bias in pairing periods of high activity on the two aerials to determine 
possible shower radiants, in that there is a natural tendency to group together 
periods of high activity which are close together in time. ‘The more objective 
type of analysis adopted here makes such a bias very improbable. 

5. The statistics of echo counts.—Levels of significance were set at 5 per cent, 
1 per cent and o-1 per cent, and a detailed analysis made of the occurrence of all 
cells (of }-hr period in time) whose count was either significantly high or 
significantly low. Data for both aerials have been pooled. For further details 
of the method of analysis, see (3). 

The result of this analysis for the 5 per cent level is illustrated in Fig. 3 (a). 
All cells in which recognized shower activity occurred have been excluded. 
The unexpected shape of Fig. 3 (a) is therefore not due to shower activity which, 
as previously noted, occurs mainly from 22 hr to 1ohr L.T.; its implications 
are mentioned below, after discussion of other possible causes. Data for the whole 
period of the survey have been combined, as there is no apparent seasonal 
influence on the frequency of significant counts. The echo rate, Fig. 3 (c), is 
also an average over both aerials. The diurnal variation in occurrence of cells 
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significant at the 1 per cent and the o-1 per cent levels is substantially the same 
as at the 5 per cent level. ‘Their absolute numbers are of course lower (‘Table 1), 
but not as low as would be expected from proportionality to the level of 
probability. Although many of the cells occupied by showers give a highly 
significant excess count, the overall contribution of showers to the total number 
of significant cells only becomes appreciable at probability levels <1 per cent. 


CELLS SIGNIFICANT 
AT SZ LEVEL 








ECHOES PER HOUR 


~ ———< 


—- ——s 


NOISY CELLS 








10 i2 i4 
LOCAL TIME 


Fic. 3.—The diurnal variation in the percentage of half-hourly cells whose echo counts are 
significantly high (O > E) or significantly low (O < E) at the 5 per cent level, compared with some other 
characteristics of the cells. Cells whose echo count is possibly reduced by interference are classified 
as noisy. 


The dip in the number of cells per half-hour around 11 hr L.T. (Fig. 3 ()), 
and a general tendency for the number of cells to be lower during the day than 
the night, are caused by a combination of a regular mid-morning film change, 
and essential equipment maintenance. Cells occupied by showers amount to 
less than 5 per cent of the whole and their exclusion from the statistical analysis 
is unnoticed in Fig. 3 (6). This is a further striking example of the insignificant 
contribution of recognizable shower activity to the total meteor echo rate 
integrated over the basic period of one year. 

One of the difficulties encountered in the application of statistical methods 
to the study of echo counts springs from cells whose counts may be in error due 
to extraneous causes. Of these causes of loss of rigour in the statistics, interference 
by noise impulses is by far the most serious. Contamination of the records by 
man-made interference, against which echo discrimination in the form of doubling 
and even quadrupling the meteor echo at times proved powerless, was a 
continuing problem during the survey, but fortunately the effect on the echo 
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statistics is determinable. All cells possibly, but not necessarily, affected by 
noise were noted; the counts for these ‘‘noisy’’ cells may either be correct or 
low, but never high. 

The average distribution of noisy cells throughout the day is shown in 
Fig. 3(d). ‘Their incidence differs markedly from that of the cells of significantly 
low count. There is also no correlation between the total number of noisy cells 
for a given month, and the total number of cells of significantly low count for that 
month. Remembering that the echo count in a cell affected by noise is not 
thereby necessarily incorrect, the absence of correlations in these two instances 
suggests that noise interference is of minor consequence to the statistics. This 
is of course borne out by the general similarity between the incidence of cells 
of high and of low count (Fig. 3(a)). It is only when we turn to the last column 
of Table I that we find clear evidence of meteor echoes lost by interference. 
Amongst noisy cells, there is some loss of significantly high counts, and a much 
stronger gain in significantly low counts which also increases with the level of 
significance. 


TABLE I 


Observed probabilities of occurrence of cells significant at different levels of expectation 


Including Excluding Percentage of 


Category showers showers 


per cent 


per cent 


“€ 


per cent 
6°4 
2°! 


0°60 


noisy 


cells 


per cent 
10 
I! 
9 


: 6:2 
I 2°0 
orl o's! 
‘Total number of cells 20380 


Whilst a high order of accuracy cannot be claimed, it seems reasonable to 
conclude that Fig. 3 correctly depicts the general pattern of the diurnal variation 
in significant cells, that this pattern is identical for high and low counts, but 
that there is a small excess, say 10-20 per cent, of high counts over low counts. 

The extent to which this excess of high count cells is indicative of permanent 
shower activity may be judged by examining the records for repetition of 
significant cells at identical times from one year to the next. ‘This comparison 
is made in Table II between 1957 and 1958 for the months of May and July. 
As before, times of accepted shower activity have been omitted. A coincidence 
is scored only if the significant cells occur at identical times in the two years. 
Expected numbers of coincidences have been calculated on the assumption 
that significant cells are uniformly distributed in time throughout the day. 
Since in fact most significant cells are grouped into about 16 hr (Fig. 3) the 
expectations in Table II could well be increased by a factor of 2, in which case 
the number of coincidences between events in 1957 and in 1958 is no more than 
would be expected on a chance basis. The fact that coincidences between cells 
of low count are no less frequent than those between cells of high count supports 
this conclusion. 
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Tasie II 


Observed and expected coincidences between significant cells for corresponding months in 
successive years 


Category 1957-58 May 1957-58 July 
Expect Observe Expect Observe 
per cent 
O>E 5 5°3 1°6 
I 06 o'09 
orl 0°06 0'005 
O<E 5 31 1°7 
I 03 0°09 
ovr 002 ©'000 


Despite this negative result, the small excess of high count cells is probably 
real, and as such must originate in minor shower activity. Regardless of the 
reality of this excess, its very smallness implies that any attempt to identify all, 
or even a major portion of, periods of high count with shower activity is 
misleading. Lists of minor showers compiled on this basis will be virtually 
meaningless. At the 5 per cent level of significance, the fluctuations in the total 
echo count, after accepted showers have been excluded, are essentially random. 
At the 1 per cent and o-1 per cent levels the fluctuations are greater than can be 
attributed to a purely random influx of meteors. ‘The excess fluctuations, which 
are strongest when the concentrations in apparent radiants of sporadic meteors 
towards the Sun and the anti-Sun lie in the collecting sectors, are presumably 
to be associated with some inherent clustering amongst the ecliptical component 
of sporadic meteors. 

6. The accepted showers.—Table III lists those showers which gave definite 
range-time envelopes, or are otherwise regarded as certain. 


TaBLe III 
Accepted showers detected from 1956 December to 1958 August 


Date of Duration Radiant Maximum 
Shower maximum (days) x r) hourly 
activity echo rate 
n-Aquarids May 8 + 2 30 
Sagittarids June 5 —35 30 
Arietids } 
f-Perseids June nad - une 60 
B-Taurids 
Phoenicids July 30 
5-Aquarids July 
Capricornids July 
Pisces Australids July }r00 
Cetids July 
Sextantids/Leonids September 29 30 
Orionids October 21 ; 8 
Puppids/Velids December 6 50 
Phoenicids December 5 20 
Geminids December 14 8 60 
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(a) n-Aquarids.—Radiant coordinates and estimates of radiant area are given 
in Table IV. There is some intermittent activity prior to May 3 and after May 10, 
which is probably attributable to this stream. The diffuseness of the radiant 
increases progressively as the Earth moves through the stream. Echo rates 
are plotted in Fig. 4; maximum activity occurs at O = 45". 


ECHO RATE 
a 
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~ 
9 








a7 re) ‘i “i 


3. 4.-——-Echo rates (shower echoes falling within range-time envelopes) for the y-Aquarid 
shower, as a function of solar longitude. 


The »-Aquarids is one of the more prominent visual showers. In comparison, 
the echo rate of 30/hr, at a zenith angle at tranist of 30°, for an equipment whose 
threshold is «= 1-3 x 10!! electrons/cm (visual brightness M,= ~7-0), seems 
rather low. However, after rejection of 40 per cent of long-enduring echoes 
on the grounds that their amplitudes showed obvious irregularities, a minimum 
of 20 per cent of the echoes detected in the 1957 apparition had durations 
exceeding 0:36sec, corresponding to «~1°4 x 10!% electrons/em, M,~ 2°0. 
The visual prominence of this shower is therefore not surprising, and the 
comparatively low echo rates obtained at the higher sensitivity can only mean 
that this stream is deficient in small meteors, i.e. the mass distribution parameter 
5 €2°0 for «max < 10! electrons/cm. 


TABLE IV 


The radiant of the n-Aquarid shower 


Radiant co-ordinates Radiant area* 
1957 1958 1957 1958 
ie s VC 
339, —2 338, —1 , VC 
338, —3 , vc 

339, —3 an c 

339, +2 f D 

340, +2 sea D 


340, +3 rene D 


* C= compact, D = diffuse 





No. 5, 1960 Radio-echo observations of meteor shower activity 397 


(6) Sagittarids.—This is a newly discovered shower of which there is no trace 
in the 27 Mc/s wind equipment records for 1952, 1953, 1954 0r 1956. ‘The radiant 
coordinates are: 1957, 307°, —35°; 1958, 301°, —36. The radiant is compact 
and the echo rate about 30/hr. Fig. 5, which is typical of the range-time plots 
generated by the weaker showers, serves to illustrate the degree of contrast 
between the background echoes and those belonging to a shower which is 
approaching the limit of resolution of the equipment. 


if june i988 
NORTH AE RIAL 


600 
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H 400} 








LOCAL TIME 


Fic. 5.—A typical range-time plot of echoes for the Sagittarids shower 


(c) June daytime streams.—The records for 1957 and 1958 June are marked 
by strong activity—echo rates up to 60/hr in excess of normal—from 06-09 hr 
L.T. on the N aerial and from 10-13 hr on the S aerial. Range-time plots 
are disappointing in that no shower envelopes can be delineated. The times 
of activity, however, require 6= +20° to +30, and this activity is therefore 
identified with the June daytime showers: the Arietids, ¢-Perseids and later in 
the month the 8-'Taurids. Although the first two showers are active simultaneously, 
the radiants differ in R.A. by 20°, and if compact their resolution should be 
within the capacity of the equipment. Failure to resolve them, and the absence 
of well-marked range-time envelopes, suggest that the radiants are diffuse rather 
than compact, especially as the Geminid shower, whose radiant is compact and 
which transits still further to the north, gives reasonably good envelopes. 

Activity near the end of June, associated with the 8-Taurid stream, was strong 
in 1958 but weak almost to undetectability in 1957. Similar variability is evident 
in previous surveys (see Paper I). 

(d) Phoenicids.—The 1958 ephemeris for this previously unreported shower 
is plotted in Fig. 6, in which each individual radiant is joined to the point on the 
mean radiant path corresponding to the longitude of the Sun at which it was 
determined. The radiant coordinates are 


%=31 ‘1+1°05( O— 109 °6), b= — 47 -9+0°54( O — 109°'6). 


For 1957, a mean radiant of 31°, ~44° was determined for July 12-17. The 
hourly echo rates of Fig. 7 indicate a marked asymmetry in activity across the 
stream, with a fairly sharp maximum at ©=112°. The equipment sensitivity 
was identical for the two years. ‘The activity is patchy and the radiant tends to 
be diffuse. 

In the 1953 survey with the radiant equipment, activity of up to 13 echoes/hr 
was observed on the S aerial on July g, 10 and 13 between 07 and 10 hr, and 
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again on the N aerial on July 16 from 08 to 11 hr. The fragmentary nature 
of the recordings did not at that time warrant the assumption that this activity 
proceeded from only one radiant. With the confirmation now available, the 
earlier records have been re-examined, and a radiant at 30°, —43° has been 
determined. Further the 27 Mc/s wind equipment recordings for 1954 and 
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Fic. 7.—Hourly echo rates for the July Phoenicid shower, as a function of solar longitude. 


1956 show an increased activity, up to 8 echoes/hr, centred on July 13, which 
could be due to a radiant which transits near o8 hr. This is sufficiently close to 
the transit time of from 06-49 to 07-12 found from the radiant equipment to 
identify this activity also with the Phoenicid radiant. 

Having now been detected during each of the five surveys made between 
1953 and 1958, this shower may be classed as permanent. 
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(e) 5-Aquarids.—Radiant measurements for 1953, 1957 and 1958 have been 
combined. A least squares fit gives for the radiant coordinates over the period 
July 20—August 10: 

& = 342°°2+0°77( O — 126°:0), 
d= — 16°-6+0°19( O — 126°-0). 


Maximum activity of 110 echoes/hr occurs at O=124°:5. The echo durations 
contain some evidence for a deficiency of more massive particles near the centre 
of the stream. 

The three minor showers (Capricornids, Pisces Australids, Cetids) whose 
presence was suspected in 1953, have been confirmed. Their combined activity 
is approximately equal to, and rises to maximum at about the same time as, that 
of the 5-Aquarids. 

A detailed report on this group of showers will be made in a separate paper. 

(f) Sextantids-Leonids.—At the time of detection of this new shower only 
the N aerial was in use, and resort to matching of theoretical envelopes to the 
echo range-time plot was necessary in order to measure the declination. For 
5~o this is, unfortunately, not a particularly accurate process. The radiant 
of 155°, O° is considered accurate to only +8° in right ascension and + 10° in 
declination. The fitted envelope is shown in Fig. 8, superimposed on contours 
of constant echo rate. The radiant area is evidently compact. Peak activity 
of about 30 echoes/hr persisted from September 29 to October 3. 

No trace of activity due to this shower could be found in previous surveys. 
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Fic. 8.—Contours of constant echo rate at the time of the 1957 Sextantid/Leonid shower. 
The numbers on the contours are the total numbers of echoes (shower plus background) received in 
10 min intervals of time for 100 km intervals of slant range, summed over the 9 days from September 
26 to October 4. The broken curve is the range-time envelope for a radiant at a=155°, =o". 


(g) Puppids/Velids.—If there were not clear evidence for this shower in 
previous surveys, and support for this view in an unpublished survey by Ellyett 
(private communication), there would be some hesitation in ascribing this 
activity (December 1-g) to a shower. The activity is intermittent, with a 
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maximum of about 50 echoes/hr near December 6, but the radiant is so ill-defined 
that range-time envelopes cannot be delineated, nor can an accurate radiant be 
determined. A tentative radiant at 140°, —50 differs considerably from the 
estimate of 124°, —36 reached from the wind equipment data (Paper II). 

(h) Phoenicids.—This shower, which made an isolated appearance on the 
night of 1956 December 5, has already been described (g). It has not recurred. 

(1) Geminids.—Five years’ observations, extending from 1952 to 1957, 
have been summarized elsewhere (10). It is worthwhile noting here that a mean 
radiant at 113-4, +31 ‘4. at © =260°-2, and a maximum activity at © = 260°°8, 
agree with the results of Bullough’s analysis (11) of the extensive Jodrell Bank 
radio measurements. 

7. The echo rates.—\t is a striking feature of the echo rates listed in ‘Table III 
that most of the accepted showers should be so close to the detection threshold 
of 20-30 echoes/hr. ‘Two reasons may be advanced for this. Firstly, the majority 
of the stronger showers have radiants to the north of the ecliptic and are either 
invisible or unfavourably situated for detection at Adelaide. ‘The second factor 
is illustrated by the comparisons of echo rates in ‘Table V. Sporadic echo rates 
are daily averages for the months in which the showers occurred. ‘The percentage 
of echoes with x>1-4x 10 electrons/cem is measured as in Section 6(a). 
Although some uncertainty is introduced by comparing data for different years, 
Table V indicates that for all the showers tabulated the increase in frequency 
with increasing magnitude is less rapid than is the case for sporadic meteors ; 
this rule applies also to the Quadrantids and the Perseids (1§). In fact, the 
only shower which seems to violate this rule is the daytime Arietids, but if we 


TABLE V 
Comparison of shower echo rates at different equipment sensitivities 


Hourly echo rates Percentage of echoes 
Shower (1) (2) (1) with durations 
> 0°36 sec. 
Shower Sporadic Shower Sporadic Shower Sporadic 
Geminids 60 33 24 39 
5-Aquarids 23 30 3°6 
7n-Aquarids 4 25 van si 
Phoenicids (July) 23 13 3°6 : - 
Years of survey 1957-58 1957-58 


Yim (electrons/cm) 3x 10" 1*4. % 30" 


TaBLe VI 
Comparison of Adelaide and Jodrell Bank echo rates 


Arietids Geminids Sporadic 
Hourly Elevation Hourly Elevation echo 
rate at transit rate at transit rate 
Adelaide 55 32 60 24 5° 
Jodrell Bank 75(a) 60 70(a) 68 8(d) 
Mass parameter s («) 2°7 1°6 2'0 
(a) Rates for 1952-53, reduced to 1950 sensitivity; references (11) (12), (13). 
(b) Reference (14). 
(c) 7>M, >5 Browne et al. (15). 
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may assume even approximately constant shower activity from year to year, then 
Table V1 suggests that in this case also the echo rate at the higher sensitivity of the 
present survey is less than would be expected from the s-values. The relative 
sensitivities of the two equipments are given roughly by the ratio of the sporadic 
rates. ‘The corrections required to reduce the echo rates of Table VI to the 
zenith are difficult to assess, but a low rate for the Arietids at the Adelaide level 
of sensitivity seems very likely. 

8. The possible showers.—It has been shown in Section 5 that the fluctuations 
in the echo rate, after subtraction of accepted shower activity, are substantially 
random, with a small excess of periods with high counts. A selection of the 
more prominent times of high echo count is listed in Table VII, as possible 
showers. ‘This tabulation includes only those events in which high counts 
appeared on both aerials with a time-difference of less than 3 hr, persisted for at 
least 1 hr in each aerial, and contained at least one cell significant at the 1 per cent 
level. These are just the characteristics expected of shower activity, but none 
of these events is classified as a certain shower because in no case could range-time 
envelopes be discerned in the plots of individual echoes. 
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Fic. 9.—Echo rates (excess over normal background count) for 1958 May 31, illustrating a possible 
shower. 

This list of possible showers is both provisional and incomplete. Some of the 
events may be chance coincidences of periods of high count which happen to 
meet the conditions stated above. Again, comparisons of records for 1957 and 
1958 for the months of May and July resulted in inclusion in this list of some 
periods which in one year alone did not meet all of the above conditions. From 
1957 August to November the equipment operated with only one channel, 
which automatically excluded from the list any high counts noted over this 
interval. 

As an indication of the type of event listed in Table VII, the echo counts 
(excess over normal) for 1958 May 31 are illustrated in Fig. g. If this is a shower 
it lasted for less than two days, as the echo counts on the following and preceding 
days were normal. 
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TasBLe VII 


List of possible meteor showers 


Approximate Approximate Maximum 
transit radiant echo rate 
.7T.) a r) per 1/2-hr 


Dec 5-6 21°42 36 — 39 12 
Mar 31-—Apr 8 10°07 340 — 32 34 
Apr 3-4 21°42 153 —39 + 
Apr 22-23 02°20 241 —44 22 


May 3-1 > 34(@) 


May 1-12 10°12 13 —18 21 
May 31 21 
May 3! 02°58 289 —18 23 (b) 


July 11 14 

July 10-11 10°20 80 — 44 9 
(a) Accompanied by an abnormally large number of long-enduring echoes, 
(6) See Fig. 9. 


Some sixty periods of high count, which persisted for at least 1 hr and 
contained at least one cell significant at the 1 per cent level, but which appeared 
in one aerial only, have also been noted. Judging by the reclassifications already 
made as a result of comparing 1958 with 1957 records, some of these events 
might be regarded as possible showers if recordings in other years were available 
for comparison. 

We may conclude that the detection of minor showers is well within the 
potentialities of the type of radio equipment used in the present survey. Their 
identification will require an objective, but perhaps not necessarily statistical, 
approach to the data and a survey extending over several years. Even when these 
requirements have been fulfilled, an effective search for minor showers is sharply 
restricted by the aperture of the system—in our case minor showers are virtually 
undetectable if their radiants transit at zenith angles greater than 30°. ‘The 
effective aperture for the major showers is of course much larger. This 
limitation as to aperture applies in some form or other to all radio equipments, 
indeed in most cases more strongly than to visual or photographic observations. 

The less rapid increase in echo rate with decreasing line density of electrons 
in the trail, for shower meteors as compared with sporadics, may suggest that the 
efficiency of a search for weak meteor showers would be improved by operation 
at low sensitivity. However, at the low echo rates so obtained the statistical 
fluctuations would be high and could easily offset the gains to be expected from 
the relatively higher shower echo rates. 

Acknowledgment.—The author wishes to thank Mr J. W. Smith for his 
great assistance in the uninterrupted operation of the equipment over the period 
of the survey. 
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THE LUNAR SURFACE AS AN IMPACT COUNTER 


E. F. Opik 


(Received 1959 October 19) 


Summary 


The dust cover of the Moon sticks to mountain slopes and cannot possess 
any degree of fluidity. It is probable that only a limited migration of the 
dust from the slopes into adjacent valleys takes place, caused by meteor 
bombardment. ‘The material of the micrometeors cannot escape from the 
Moon and is expected to form a layer of the order of 40 gr/cm?, accumulating 
over the maria and protecting them from further erosion by small meteors. 

Large meteorites are not stopped by the dust and produce meteor craters; 
at a velocity of 20 km/sec the ratio of crater to projectile diameter is estimated 
to be about 20. 

812 crater diameters, measured over an area of 465,000 km? in Western 
Mare Imbrium, show a frequency distribution which is close to that predicted 
from the present population of interplanetary space and a time interval of 
4°5 x 10" years. It is concluded that since the formation of Mare Imbrium 
the population of interplanetary space has not changed appreciably; a tentative 
explanation is given. 





1. The nature of the lunar surface and crater formation.—The craters of the 
Moon can be regarded as imprints testifying to the last stage of the building-up of 


this, and probably other members of the terrestrial group of planets. Their 
origin by way of impact stands practically beyond doubt, as contrasted by a few 
small cones pointed out by Kuiper (1954, and verbal communication 1959) whose 
relief appears to be similar to terrestrial volcanoes, and which may be either the 
relics of early volcanic activity, or the outcome of secondary processes following 
the impact of bodies of asteroidal dimensions. The eruption of gases from the 
central peak of Alphonsus as described by Kozyrev (1959) may well be an example 
of such long-enduring after-effects. On the other hand, the topography of the 
Moon lacks any indications of orogenic processes of the terrestrial type, such as the 
still continuing Alpine foldings and concomitant volcanoes along the tectonic arcs. 
Nor could such activity be expected to take place on the Moon, on account of the 
practical absence of erosion; the orogenic cycles on Earth are the result of 
unbalancing by large-scale erosion and sedimentation, absent from the Moon. 

Of course, the possibility of large-scale erosion has been suggested, in the 
form of creep or flow of lunar dust (Gold 1955). A critical consideration of 
observational evidence indicates, however, that large-scale flow actually does not 
take place. The presence of dust (probably due to battering by meteors) is 
revealed by bolometric observations of temperature changes during lunar eclipses 
(Pettit and Nicholson 1930, 1940), which indicate a layer at least 5 cm thick of low- 
conducting material which can only be dust in vacuo, or a similar porous substance 
(Wesselink 1948). The observations of 1927 (Pettit and Nicholson 1930) were 
made on a spot in the rugged region 48” from the southern limb; for purely 
geometrical reasons, only mountain slopes turned toward the observer could have 
contributed to the measured radiation. ‘The presence of the dust layer on these 
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slopes, as revealed by the bolometric curve, refutes its fluidity postulated by Gold ; 
whatever the cause, the dust sticks to its place of origin (as indicated by some 
experiments, dust particles in vacuo become welded together at the points of contact, 
which may offer an explanation), being disturbed only by meteoric bombardment ; 
and this can cause displacement but on a very limited scale, a gradual diffusion from 
the mountain slopes into the nearest valleys. It can be estimated that, with the 
actual content of micro-meteoric material in interplanetary space, and a ratio of 
mass eroded to meteoric mass of about 100, a dust layer of the order of 20 metres 
could have been battered out of an exposed rock surface during the lifetime of the 
Moon (assumed to be 4°5 x 10° years). ‘The estimate follows from a theory of 
meteor impact (Opik 1936, 1958a) with 5 km/sec as the impact velocity, and an 
accretion of 0-01 gr/cm? per one million years, or one-quarter that estimated for the 
Earth (Opik 1956). However, a layer of 0-2 gr/cm? would already protect the 
underlying surface against further erosion by micrometeors, so that, if the dust is 
not removed, subsequent erosion could have been caused only by the larger and 
faster meteors whose total mass is relatively small; an upper limit of some 100 cm 
for the eroded layer can be set in this case. 

The dust may be supposed to consist of a mixture of locally pulverized surface 
rock, and of the products of destruction of micrometeors. These, identified as 
zodiacal particles, at a velocity of impact of the order of 5 km/sec (Opik 1956), 
may become partly or completely vaporized at impact, but the velocities of the 
vapour molecules will be mostly less than the velocity of escape (2:3 km/sec), 
because most of the kinetic energy will be absorbed by the medium or go into the 
heat of vaporization. ‘The vaporized material will return to the lunar surface and 
be deposited there, providing thus another means of cementation of the lunar dust. 
A layer of meteoric material of the order of 40 gr/cm* may have been produced in 
such a manner during the lifetime of the Moon. ‘This layer would suffice to cover 
the Moon with a uniform coating, efficiently protecting the rock surface from 
further erosion. 

The differences in reflecting power of different portions of the lunar surface 
indicate, however, that such uniformity has not been achieved. If we assume 
that, as the result of meteoric impact, the mountain slopes of the continents are 
currently relieved of the excess of dust which collects in the adjacent valleys, there 
will be intense erosion of the slopes and the dust will mainly consist of local material, 
with but a moderate proportion of meteoric substance intermingled; this will 
refer also to the valleys between the mountains where the material from the slopes is 
collected. On the contrary, in the extended plains of the maria the dust remains in 
situ, erosion is prevented by the protective layer, and the amount of deposit of 
meteoric origin will greatly exceed that originally ground off from local rock. 
If the meteoric deposit is dark, the rock powder bright, the difference in color- 
_.ation between the continents and the maria can be explained. We note that the 
, conclusion about the continuing erosion of mountain slopes necessarily follows if 
a considerable fraction of the meteoric material is retained by the Moon. 

As to the ray systems of Tycho, Copernicus and others, they may be due partly 
to spray of a light-coloured material, partly to a ruggedness of their surface, 
perhaps to cracks whose slopes are still subject to erosion. 

In the eclipse of 1939 a central region of the Moon was observed bolometrically 
(Pettit and Nicholson 1940). The temperature curve indicated thermal properties 
of the lunar soil that were similar to those of the previous observation. About 25 


30* 
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per cent of the observed area was covered by a mare surface. As pointed out by 
Gilvarry (1958), 25 per cent of bare rock (lava) surface would have considerably 
influenced the temperature variation, making it less steep; hence he rightly con- 
cludes that the maria are covered with dust. However, contrary to his further 
conclusions, this does not contradict the lava theory of the maria; a lava surface, 
like that of any other rock, if exposed to the vacuum and meteorite impact, will 
ultimately be covered by a dust layer. The eclipse observations indicate only the 
existence of the thermally insulating layer of dust, without producing evidence 
relating to the nature of the underlying rock or rubble. 

It appears, thus, that the whole surface of the Moon is covered by dust due 
mainly to the bombardment by meteors. On mountain slopes the layer must be 
thin, in the maria much thicker but hardly more than the equivalent of some 
40gr/cm®. ‘To meteorites exceeding 100cm in diameter such a dust presents an 
insignificant obstacle, and craters will be produced as at impact into solid rock. 

TABLe | 
Depth of erosion and crater diameter, in units of the spherical equivalent diameter 
of the projectile 
(Assumed density of surface rock=2-5; of stone meteorite=3°9; of iron=7-8) 
Stone Meteorite Iron Meteorite 
Velocity 20 km/sec 5 km/sec 20 km/sec’ 5 km/sec 
Depth of erosion 2°5 2°2 40 3°6 
Crater diameter 18 9°5 20 10°5 
From the theory of meteorite impact into solid surface (Opik 1958 a), the data as 
giyen in Table I are derived. ‘The crater diameter (measured from rim to rim) as 
determined by the volume of crushed rock is not supposed to depend much on 
gravity, but for large craters more rubble will fall back and remain inside the crater 
which will therefore appear shallower than small craters. We note that there 
is but a very slight dependence of the crater to projectile diameter ratio on the 
density of the projectile. This is because the increase in crater volume, due to 
increased density of the projectile, is achieved mainly at the expense of depth, and 
only to a minor degree of the area of the crater. It is admissible, therefore, in the 
statistics of lunar craters, to deal solely with diameters, disregarding the densities 
or masses of the projectiles. 

The figures are supposed to represent a close-order estimate, perhaps to 20 
per centin the linear dimensions. From a special analysis (Opik 1958 b) it appears 
that crater-producing collisions are at present expected to take place chiefly with 
the members of the Appollo group of asteroids, and with comet nuclei, predomin- 
antly non-periodic comets, both groups being represented in almost equal numbers. 
The average impact velocity may be estimated roughly to be 20km/sec for the 
first group, and higher for the second. We may assume 20km/sec in Table I, the 
higher velocity of the comets being compensated by the smaller density of their 
nuclei, and conventionally accept a ratio of crater diameter to projectile diameter 
equal to 20, with the understanding that the mass of the projectile should be 
calculated with a density of meteoric stone, i.e. of 3-9. 

The continents are much more densely covered with craters than the maria. 
According to current concepts this means that the contintents correspond to an 
earlier stage as compared with the maria. These latter are believed to be caused by 
impacts of large projectiles toward the end of formation of the Moon’s surface 
features; in this process, the previously formed topographic features (still partly 
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recognizable in the so-called ghost craters) have been levelled out. The craters of 
the maria are thus formed during the period starting with the end of formation 
of the Moon’s surface and ending with present time, therefore reflecting the history 
of the Moon under conditions of surrounding space which were more or less 
similar to those prevailing now—contrary to the more chaotic conditions of the 
preceding primitive stage. The statistics of craters in the maria can thus be 
directly compared with what is known of the frequency of stray bodies in the 
vicinity of the Earth and, among other things, can be used as a quantitative check 
of the mechanical theory of formation of meteor craters. Needless to say, quanti- 
tative tests in such matters can be relied upon only as to the order of magnitude. 

2. The distribution of crater diameters in Western Mare Imbrium.—Rim to rim 
diameters of craters were measured and counted by the author on a print of the 
well-known photograph made with the roo-inch Mt Wilson reflector on 1919 
September 15, (scale 1mm=4:32km); by courtesy of Mt Wilson Observatory, 
the print was reproduced in the Irish Astronomical Fournal (Fitzgerald 1953). 
Craters down to 1:2 km diameter appear to be completely recorded on the original 
print; on the half-tone reproduction in the Journal, made to three-quarter size, 
the limit of completeness is about 2-4km. The limit of recognition of small craters 
is uniform all over the photograph, irrespective of the distance from the terminator 
(which passes just beyond the western edge of Mare Imbrium). The diameters 
were measured with a millimetre scale, in the direction of the major axis of the 
projected crater image. Allowing for the conditions of illumination, subjective 
judgement had to be applied, to determine the position of the crater rim ; it was felt 
that no systematic errors of importance could have originated from this source, as 
well as from resolving power ; the latter conclusion was confirmed by re-measuring 
the diameters on the half-tone print. The probable error of the scale readings was 

+011 km or + 0"-05, and this also may be assumed to be the probable error in the 
diameters of the small craters. For larger craters the error may be greater, due to 
judgement errors caused by irregularities in the rim topography. Altogether 812 
craters ranging in diameter from 0-6 to 73-4km were recorded. The counts were 
made along twelve strips stretching in an east-west direction, the southernmost 
strip No. 1 running just north of Eratosthenes, the northernmost strip No. 12 
south of Plato and bordering on Sinus Iridum but excluding the latter. 
‘* Epicontinental’’ areas such as that south-east of Archimedes and the rugged 
region around Autolycus were also excluded but the floor of Archimedes with 
eight counted craterlets was included. Areas completely immersed in shadow 
(Cassini) were also excluded. The areas quoted in the second column of Table 
II refer to the portions of Mare Imbrium in which counts were actually made. 
In the third column, cos « is the ratio of projected to true area of the strip. In the 
fifth column, p, the recorded surface density of craters of all sizes, shows a 
systematic decrease from south to north which far exceeds the sampling error. 
The southern edge of Mare Imbrium is more thickly populated by craters than 
the rest of the surface of this mare; this may be due to a greater proportion of 
objects of the pre-mare stage having survived there. There is no significant 
difference in the relative number of small and large craters (below and above 2°5 
km diameter) in the southern portion of the mare, as compared with the rest. 

Table III summarizes the distribution of crater diameters for the entire region. 
From the general run of the counts it appears that the data are complete down to a 
diameter of 1-19km. 
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TaBLe II 
Number of all craters recorded 


n=observed number; S=area of strip, in square kilometres; ‘COS a= average projection 
ratio; p+p.e.=number per 10,000 km® and its sampling error + p/\/2n 
Strip No. S, km* cosa n prp.e. Remarks 
I 20,000 0°93 «88 44+3°5 Southernmost strip 
21,300 o9g2 64 30+3°0 Contains half-destroyed crater wall, diam. 
24°2 km (counted) 
18,300 o°92 40 22+2°5 
24,500 o89 62 25+2°5 
28,600 o88 58 20+2:0 Contains Timocharis and Lambert. 
39,800 0°87 78 20+1°5 Contains Archimedes and Autolycus. 
54,500 o85 71 13+1°0 Middlestrip; contains Aristillus. 
51,700 o'82 50 10+1°0 
9 44,600 o80 84 I9gtI°'5 
10 49,100 7o 14+1°5 Contains Piton; Cassini excluded. 
II 49,800 69 14+1°5 Contains Leverrier and Helicon. 
12 62,500 78 12+1°0 Northernmost strip, with Pico. 
Total 464,700 81217°5+0'4 


Tasie III 
Distribution of crater diameters in Mare Imbrium 


area = 465,000 square kilometres; »=number within limits; N= cumulative number 


Diameter, km n N Diameter, km n N 
x ° 5°40 
19 
70°6 4°10 


21 


55°5 3°35 
36 
97 
135 


75 


1°19 
7°99 0°97 


6°70 0°64 


5°40 35 0°00 812 


The number dn within the limits of diameter of the craters, or the projectiles 
which caused them, from x to x + dx can be assumed to be defined by the conven- 
tional frequency function 


dn=Cx~? dx, (1) 
where the exponent p or the ‘‘ frequency index ’’ is normally a positive number. 
The expression is an idealization which cannot be conveniently applied to small 
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counted numbers, on account of sampling fluctuations. Its integral from 
infinity to x can be used instead to represent the cumulative number, N, with a 
satisfactory degree of stability even in the case of very small numbers, 


N=Nox-* (2) 
where s, for obvious reasons always positive, can be called the population index *. 
In the case of strict validity of equation (1) we have evidently 


ds ate (3) 
an idealization which, however, seldom applies in real cases. It is preferable to 
define the population index directly by eqn. (2), even with s variable according to 
circumstances. Usually but a slow variation of s withx would ocur. The popula- 
tion index reflects the history of the population and its mode of origin. Thus, 
fragmentation will tend to increase the number of small objects, and to yield a 
large value of s (meteors) ; condensation, on the other hand, leading to the building 
up of large masses (planets, stars), would tend toward small values of the popula- 
tion index. The population index of the Mare Imbrium craters (or their projec- 
tiles), as compared with that of other populations (from various sources, partly 
unpublished), is given in Table IV. 


TaBLe IV 


Population index (s) of diameters (x) 
Population x stp.e. 
Mare Imbrium craters 1°14- 2°48 km 1°70 + 0°07 
2°48- 5:40 km 2°24 0°15 
§°40-12'7 km 1°49 + 0°22 
12°7 -34°3 km 1°22+0°35 
34°3 -70°6 km 1°52+0°50 
1°14-34'3 km 1°64+0°10 
Micrometeors 10-4107? cm 1°6-1°8 
Visual Meteors 
In cometary orbits 2x 107! cm 3°3 
In asteroidal orbits { —- I cm 472 
Asteroids 2-90 km 1°59 
90-360 km 2°5 
Appollo Group asteroids 0°4-2°8 km (2°7:) 
Comet nuclei >6km 21+0°2 
Meteorites 50-10‘ cm 3°3 
In the Solar System 
Planets and 200- 4,000 km o's 
Satellites \ 4,000-64,000 km o's 


The population index of the lunar craters is 1-6, or nearly equal to that of the 
smaller asteroids; the range of the projectile diameters is, however, very different, 
from 0-06 to 1-7 km for the hypothetical lunar projectiles, thus linking up, at their 
upper limit, with the smallest observable asteroids. 

3. Comparison with the contemporary population of the Earth’s vicinity.—From 
an analysis based on observed data bridged over by inter- and extrapolation (Opik 
1958 b) the frequency of collisions of the Earth with objects of different size has 
been estimated. Ata relative velocity of about 20 km/sec the condensing action of 


* The population index, under the name of “ gradient’’, has been successfully used by 
Vaucouleurs (1956) in the analysis of stellar populations. 





410 E. }. Opik Vol. 120 
the Earth’s gravitational field amounts to an increase in the number of impacts in a 
ratio of only 1-3 toone; for the Moon the frequency of impacts per unit area may be 
thus assumed equal to 80 per cent that for the Earth. In Table V is given the 
estimated cumulative number of impacts (Ntneor) upon a lunar area of 465,000 km? 
over a period of 4,500 million years, corresponding to the present-day frequency 
(Opik 1958 b); this is compared with the cumulative number (Nops) of craters in 
Mare Imbrium registered over the same area; the projectile diameters are 
assumed equal to one-twentieth of those of the craters. 


TABLE V 


Cumulative number of impacts on Mare Imbrium (465,000 km*) in 4°5 X 10° years 
Diameter of projectile, 
lower limit, km 0060 O'124 0270 0°64 1°72 
Nobs 733 208 35 10 3 
Nuheor 1050 202 35 5'0 0°44 
Ratio Nopg/Ntneor 0°70 1°03 I'o 20 7 


The observed numbers are in close agreement with the theoretical figures 
except for the large craters whose numbers are in excess of the expectedones. The 
agreement is closer than could have been expected from the character of the 
estimates. The assumed crater to projectile diameter ratio, as well as the fre- 
quency of meteoritic bodies in the diameter range 0:06-0:27 km are based on rough 
' theoretical and very scanty observational data ; the very close agreement can, there- 
‘fore, only be a matter of chance. Nevertheless, the conclusion appears to be 


justified that the craters in Mare Imbrium can be accounted for by bombardment 
at a more or less uniform rate during 45 x 10° years by stray bodies of the same kind, 
and of the same space density, as those populating the solar system in the vicinity 
of the Earth’s orbit at present. The observed excess in the numbers of large 
projectiles, or large craters in Table V can be ascribed to survival of the larger 
formations, which may have withstood the events connected with the formation of 
the mare; they belong thus, probably, to the pre-mare stage. 

The conclusion is somewhat unexpected; the stray bodies are continually 
swept away by the planets, and their number is expected to decrease exponentially 
with time; in the beginning there should have been very much more stray bodies 
than at present, and the number of impacts should have been considerably greater 
than that calculated from the present-day density. 

However, a plausible explanation can be given to this puzzling result. The 
accretion of the planetary masses most probably proceeded at the expense of 
debris and gas placed more or less in the region where the planetary orbits 
ultimately stabilized. The projectiles which battered the surface of the Moon in the 
pre-mare stage must have belonged to the same category, circling the sun in orbits of 
small eccentricity and inclination, say e= 0-1, sini=0°1, a 1 astr. unit. Formulae 
of collisional probabilities (Opik 1951) yield in such a case a collisional lifetime 
with respect to the Earth of the order of 2 x 10° years. These ‘‘ planetesimal ”’ 
fragments must have been removed rapidly long ago, and no trace of them remains 
at present; the surrounding space was swept clean of them already at the time of the 
major collisions which produced the lunar maria. [These collisions may have 
been caused by moonlets which, like the Moon itself, originally, may have been 
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formed as independent planets on the Earth’s orbit (cf. Gerstenkorn 1955) 
drifting very slowly toward each other (because the periods of revolution were 
almost equal) and colliding after the other fragments were swept away. ] 

The lifetimes of the objects of the Appollo group, of the order of 10% years, are 
also too short for them to be regarded as remnants of the primordial fragments ; 
they rather are continually replenished from outside, to replace those swept away, 
and an equilibrium density of them has established itself in our surroundings (Opik 
1951). A source of these objects must exist which is not easily depleted. The 
source may be partly in the asteroidal region, whence cumulative perturbations 
in close approaches to Mars will deflect some of them toward the inner portions 
of the solar system. Another, more inexhaustible source is in the cloud of comets 
surrounding the solar system to distances of the order of 10° astron. units (Oort 
1950); from the outskirts they are efficiently removed by stellar perturbations 
(Opik 1932), but in their bulk they are very stable. A comet of a= 10 astron. 
units which happens to enter the inner regions of the solar system has a collisional 
lifetime of the order of 10!* years, thus very much longer than the solar system ; 
the lifetime of the Mars asteroids is of the order of 10 years. With such lifetimes 
the supply of stray objects toward the inner portions of the solar system may have 
remained constant, indeed, during all the time of the existence of the solar system 
since its rapid condensation. 

If this is so, the surfaces of the lunar maria can be regarded as a counter which, 
over the ages, has registered the impacts of stray bodies and from which conclusions 
can be made on the present-day distribution of diameters and number density of 
these bodies. 


This research was supported in part by the Geophysical Research Directorate 
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Summary 


A review of existing methods of determining the strength and direction 
of magnetic fields in sunspots, as a function of position and depth within the 
spot shows that these have given fragmentary and somewhat discordant 
results. In particular, little is known of the direction of the magnetic 
vector. 

A new method of determining this quantity is described. By means of a 
Babinet compensator, polarization fringe systems are produced in the 
components of a normal Zeeman triplet. From the phase difference of the 
o-component fringes, the constants of the elliptical polarization at each point 
in the spot are determined, and the orientation of the lines of force deduced. 
The method minimizes the effect of instrumental polarization; the reduction 
is simple and fairly short. 

An experimental test on four sunspots has been carried out, using the 
Oxford 35m solar tower and 12m spectrograph with a Babcock grating. 
Measurements of field direction were made, using the triplet 6173 (Fe 1); 
the associated field strengths were determined from the same line and from 
A6149, a weak line of Fett. The resulting field-strength variation is in agree- 
ment with Broxon’s parabolic law, but in contrast with Joy and Nicholson’s 
measures of direction, 40 per cent of the flux is found to be confined within a 
cone of unit solid angle. Comparisons of field-strength at A6173 (Fe 1) and 
A6149 (Fe 11) suggest that the gradient in depth does not exceed 1 oersted/km. 





Introduction.—It has long been recognized (3) that the physical processes 
operative in sunspots present a problem in magneto-hydrodynamics, for the 
correct investigation of which it is necessary to have an accurate observational 
picture of the magnitude and direction of the magnetic fields and the velocity 
fields in sunspots, as a function of position and depth within a typical spot. 
It is now just fifty years since Hale first interpreted the splitting of sunspot 
lines in terms of the Zeeman effect, and so laid the foundations of quantitative 
magnetic measurements of sunspot fields (7, 8), and since Evershed (4) discovered 
the effect named after him, in which he interpreted variations in the wave-length 
of the same line in different parts of a spot in terms of radial and tangential velocity 
streams. During this time much work has been done in the development of 
techniques and in systematic measurements and yet it cannot be said that the 
present position regarding these basic data is at all coherent or satisfactory. 

1. Existing methods and results.—The methods which have been used almost 
universally for determining the magnetic fields of sunspots have been derived 
from procedures described in a classical paper of Hale and Nicholson (9). 


* Skynner Research Fellow in Astronomy, Balliol College, Oxford. 





No. 5, 1960 The spatial analysis of magnetic fields in sunspots 413 


To determine the field strength, the observation consists essentially of the 
measurement of the splitting of a Zeeman triplet, from which the total magnetic 
field in which the line originates is determined by the relation 


A\=gCHX (1) 


where A denotes the wave-length, g the splitting factor, C is a constant, and 
H the field strength. The practical difficulty in this basically simple measurement 
is due to the fact that the splitting for average spot fields is of the same order of 
magnitude as the line width, so that even with high dispersion and resolving power 
the components are only properly separated at the centre of large spots. Hale 
solved this problem by reducing to a practical form a procedure first used by 
Zeeman to detect small magnetic splittings in laboratory spectra. In the case 
of a normal triplet, viewed along the field, the central 7-component is absent 
while the outer o-components are circularly polarized with opposed directions 
of rotation. Consequently one or the other can be extinguished when the triplet 
is viewed through a quarter-wave plate and a Nicol whose plane of transmission 
makes an angle of + 45° to the wave-plate axis. In this way the difficult measure- 
ment of an unresolved doublet can be replaced by the determination of the 
relative positions of two undisturbed lines. Hale constructed a mosaic of 
quarter-wave strips so oriented that alternate components appeared in adjacent 
strips. ‘This method permits rapid visual or photographic measurements, and 
forms the basis of a long series of field measurements by Hale and Nicholson 
and their successors at Mount Wilson. In an important paper (31), von Kliiber 
refined the method of Hale and applied it to extensive photographic measure- 
ments, and his techniques form the basis of a series of determinations of field 
strength profiles by himself (31), by Schréter (23) and by Mattig (18), and 
of a programme of regular maximum field measures which is now in progress 
at Potsdam (17). Other methods of determining the Zeeman splitting have been 
suggested, including the accurate photometry of the blended triplet (31) and the 
use of Savart fringes to replace the multiple quarter-wave plate (19), but no 
measurements with these techniques have yet appeared. 

The problem of determining the law of variation of the field with depth 
has been fully discussed by Houtgast and van Sluiters (10). Measurements of 
field strength made from a strong and a weak line in the same spot, or by means of 
the same line, when the spot is near the centre of the disk, and again when it is near 
the limb, provide, in principle, means of comparing fields which exist at different 
optical depths. On reasonable assumptions concerning the physical conditions 
in the sunspot, Houtgast and van Sluiters determined the mean linear depths 
at which particular lines are formed, and hence provided two methods for 
determining the gradient of the field. Yet another method is suggested by them, 
which assumes knowledge of the direction of the lines of force. If we assume 
that the lines of force are straight over the observable depths of the sunspot 
(an assumption already implicit in assigning a unique value to dH/dh), the value 
of the gradient can be immediately calculated from the divergence of the lines 
of force. 

The classical work on the determination of this law of direction is due to 
Nicholson (g), who used two quite different procedures. The first depended 
on a comparison of the intensities of the Zeeman components of a triplet, when 
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observed with the quarter-wave plate and Nicol. According to Seares (26), the 
intensity of these two components is then given by 


,=}(1 — cos yo)? 

a= }sin* y, (2) 
o,=}(1+ cosy)” 
where +, is the inclination of the magnetic vector to the line of sight. Assuming 
the validity of this law, if it were possible to measure the intensities, corrected 
for blending, by accurate microphotometry, the directions of the lines of force 
could be deduced. Nicholson noticed however that when y)=9g0°, o,=0,. 
By observing the point in a spot at which these components were of equal 
intensity, and repeating the observation as the spot moved towards the limb, 
he was able to synthesize the directions of the vectors over a considerable part 
of the spot. 

His second method made use of the fact that the 7-component is always 
plane polarized in a plane containing the line of sight and the line of force. Near 
the limb of the Sun, this implies that the direction of polarization of the 
m-component coincides with the direction of the lines of force. Nicholson 
determined this direction at different points in spots near the limb by extinguishing 
the central component with a Nicol and rotatable half-wave plate, and so extended 
his observations to near the edge of the penumbra. His pioneer work in this 
field has not been repeated. 

A source of systematic errors which has to be reckoned with in applying the 
methods developed at Mount Wilson is introduced by the oblique reflections 
which occur at the coelostat mirrors. The instrumental polarization thus created 
(26) results in a partial reappearance of the unwanted o-component in the method 
of determining field strengths, and a consequent displacement in the apparent 
centre of the measured component. The effect is even more serious however 
in the methods of determining field direction since these depend directly and 
critically on the intensity relations of the c-components, and on the form of their 
polarization, both of which quantities are distorted in a rather complicated 
manner by the instrumental polarization. 

Summary of results.—The present position of the problem may be seen by 
placing together the main results of the application of the techniques described 
above, which constitute the elements of our present picture of sunspot fields. 
Broxon (2) published in 1942 a much-quoted law of field intensity, which is 


believed (3) to represent unpublished results obtained by Nicholson. This is 
the parabolic law: 


H(r) = Ho(1 — (r°/6?)) (3) 


where H(r) denotes the field strength at distance r from the spot centre, Hy is 


the central field, and 6 is the penumbral radius (about 2-4 times the umbral 
radius). 


The measures of Mattig yield a somewhat different result; the field was here 
measured very carefully to the edge of the spot, and was found to disappear 
more gradually than the parabolic law indicated, being better represented by 


H(r) = H,(1 — 4 [b*) er, (4) 





No. 5, 1960 The spatial analysis of magnetic fields in sunspots 415 


The normally quoted values of the vertical field gradient are mainly those of 
Houtgast and van Sluiters (10), who discovered the centre-limb variation in 
field intensity in the Mount Wilson results (9g), and who applied the measurements 
made by R. B. King (14) on Mount Wilson plates on lines of different strength. 
They find, from the centre-limb effect dH/dh=5-7 oersted/km; from the 
comparison of weak and strong lines dH/dh=2-5 oersted/km, while from the 
divergence of the lines of force according to Nicholson (g) (see equation (5) 
below) they deduce a field of 0-5 oersted/km. In an independent investigation 
of lines of different strength Tanaka and Tagaki (29) found a barely detectable 
gradient. 

Nicholson’s investigation of the directional properties of the lines of force 
by the study of the positions of equal intensity of the o-components gave a graph 
which can fairly be represented by 


Yo=0°75 4n7/b (5) 


while his determination from the orientation of the plane of polarization of 
m-components in spots near the limb is well represented by 


Yo= hn r/b. (6) 


Chapman (3) suggested this second formula as representing roughly the combined 
results of Nicholson’s measure, which are the only published data on line 
directions. 

It is clear that where independent methods and observations are available, 
there is considerable disagreement in the resulting values. This fact has given 
rise to various criticisms and suggestions. Hale shows that the quarter-wave 
plate method of field strength measurement, in which the assumption is made 
that the central component is absent, can be used validly, for fields where the 
value of y, does not exceed 20°. Now it is clear that if Nicholson’s law of 
directions is even a first approximation to the real condition, the method can only 
be applied without further correction to a region near the centre of the umbra 
in spots near the centre of the disk. The situation is complicated where the field 
direction is irregular, or where the third component is provided by scattered 
light from the photosphere round the spot, a factor always present in observations 
near the limb. On grounds of this kind Thiessen (30) has criticized the Potsdam 
results, which all suggest the exponential form of equation (4), and believes that 
a correction for these effects tends to restore the parabolic form. More radical 
criticism has been made by Stepanov (28) who considers that Seares intensity 
law (equation (2)) is not valid, in the case of a line as broad as (6173, for field 
strengths below 800 gauss, but must be replaced by a more complicated law 
derived by detailed consideration of the transfer problem. Recent measures of 
spot fields by Severny (27) using a solar magnetograph, reveal moreover fine 
structure and irregularities of polarization in spots. It is a little difficult at 
present to assess the physical meaning of these irregularities, but it seems possible 
that variations of field direction are involved, a quantity which a magnetograph 
is better equipped to detect than to analyse. If these irregularities are typical, 
the conditions under which the quarter-wave plate technique is valid are further 
delimited. It must however be admitted that the results achieved with this 
method have provided a fair general picture at least of the umbral field. 
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The disagreement, amounting to a factor of 10, between the values of the 
vertical field gradient, is rather disappointing, particularly when we note that 
the intermediate value is given little weight by Houtgast and van Sluiters on 
account of the discordant results of Tanaka and Tagaki (29) who found large 
differences between lines of different elements. Since it seems most unlikely 
that the gradient computed from the law of directions should be too small by an 
order of magnitude, one is forced to consider that the centre-limb variation is 
influenced more than was suspected by scattered light and by the presence of the 
central component. 

If this is the case, then we depend on the accuracy of the law of directions for 
our knowledge of the vertical field gradient. Now it is true that there is broad 
agreement between the two types of measurement performed at Mount Wilson. 
But it must be remembered that the measures are extremely difficult ones, 
involving measures near the limb of the Sun and near the edge of the spot where 
the components are imperfectly resolved, and that they are subject to systematic 
errors due to instrumental polarization. The most serious objection to these 
results lies however in the difficulty of reconciling the law of direction with the 
results of the Evershed effect. Studies of this effect by Evershed (4, 5, 6), 
Abetti (1) and others and notably a recent study by Kinman (15, 16) indicate 
that though irregular transverse components are sometimes present, the 
velocity-field is mainly radial. Now Jensen (11) has pointed out that if, as 
seems most likely, the motions in sunspots are constrained to the lines of force, 
these results cannot be reconciled with Nicholson’s law of direction. 

It appears from this survey of existing methods and results that the basic 
data of sunspot fields are insecure at almost every point. More observations are 
required in which the limitations of the various methods are carefully respected 
and the development of other methods, free from the peculiar systematic 
limitations of the techniques developed at Mount Wilson, is desirable in order 
to provide an independent check. In this paper, we now present an account 
of the development and application of a new method of determining the law of 
direction. 

2. The Babinet fringe method of field direction determination.—The principle 
of the new method is a simple one. Consider the circularly polarized components 
of a normal Zeeman triplet. Viewed along the field, the path of the amplitude 
vector is in fact a circle. Viewed normally to the field, one of the orthogonal 
components of the circularly polarized light is in the line of sight and ineffective, 
so that the light observed is plane-polarized in the plane of the field (transverse 
Zeeman effect). It is easy to see, qualitatively, that at any intermediate angle, 
the light will be elliptically polarized, and that the ratio of the major to minor 
axes of the ellipse will determine the angle y, between the sight-line and the 
field, while the orientation of the axes will determine the azimuth of the projection 
of the field on a plane perpendicular to the sight line. Hence a complete 
determination of the two constants which describe the shape and orientation of 
the ellipse also determines the orientation in space of the magnetic field. 

In order to express this quantitatively and to devise a method of performing 
the analysis, we recall some basic relations between equivalent analytic descriptions 
of polarized light. Consider light whose state of polarization may be described 
by an ellipse of major and minor axes A and B, the major axis being inclined to 
arbitrary orthogonal axes at an angle y. (y is wrongly indicated in Fig. 1. It 
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should be measured from the major axis.) It is shown, for example by Schuster 

(25), whose notation we follow, that such light may be represented by 

perpendicular vectors a,, a, directed along the arbitrary axis of reference, and 

differing in phase by 5. If we set B/A=tan'Y, a,/a,=tany, then the following 
relations hold: 

tan 6=tan 2‘ /sin 2y, (7) 

cos 24s = cos 2y cos 2, (8) 

sin 2‘Y = sin 2 sind, (9) 

tan 2y = cos 6 tan 2y, (10) 

cos 2'¥ sin 2y = sin 25 cos 4. (11) 

Note that in the particular case in which the axes of reference coincide with 
the axes of the ellipse (y=0), by equation (7), 5=47. Schuster (25) shows 
that the sign of 5 governs the direction of rotation of the ellipse. 

To obtain the quantitative relation between the constants of the ellipse and 
the direction of the field, consider Fig. 1, in which the circle of polarization is 
represented in perspective, OE is the sight line, OH the field direction, 
EOH=yYo; OH is normal to the plane of the circle. The circular polarization 
may be defined by a vector OD=Asinwt, in the line of intersection of the 
plane EOH with the circle, and a vector OF = A sin (wt + $7) in the plane of the 
circle, perpendicular to OD and therefore to OE. Now since the vector OF 
has no component along OE, it is observed unchanged for all values of yo, 


Fic. 1.—Relation between polarization and direction of the lines of force. 


whereas the effective component of OD in the plane perpendicular to OE (the 
plane tangential to the celestial sphere) is 


OG = A cos yg sin wt. 


These two perpendicular components OF, OG will produce at E light elliptically 
polarized, in such a way that 


COS yp = tan ‘Y (12) 


while the direction of the minor axis OG in the plane perpendicular to OE 
measured from a suitable origin (the North Point) defines the azimuth of OH 
with respect to this origin. 

2.1. Determination of the elliptical constants.—In order to determine the 
quantities tan '¥, y, which describe the shape and orientation of the ellipse, we 
perform an analysis by means of a Babinet compensator and a polarizer. The 
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method here described differs somewhat from the conventional application of 
a Babinet compensator, since it avoids any measurement which involves 
assessment of intensity or fringe visibility, and yields the quantities tan'Y, y, 
more directly. 

The Babinet compensator consists of two similar quartz wedges, of equal 
angle, cut with their crystal axes normal to the direction of incident light, and 
respectively parallel and perpendicular to their prism axes. The prisms are 
cemented together to form a rectangular block (see Fig. 2). We may represent 
incident elliptically polarized light by vectors in the directions of the crystal axes. 


a, cos © sinwlt 

+a, 5.n6 sin (wt+ $+ cx) 
POLARISER 
TRANSMISSION 4 
PLANE. 
BABINET t ail 
PRISM AXIS ee w 


\ ey sin(wt +5+ cx) 





———_ Qa, SLN (wt + 5) 
Fic. 2.—Analysis by Babinet compensator. 


One of these components is the ordinary ray in the first prism and the extra- 

jordinary ray in the second; the prisms act the opposite way for the other 
_ component. Consequently, the action of a compensator of prism angle « 
(in radians) ona ray normal to the prisms and distant x from the point of equal prism 
thickness is to introduce between the two components an additional phase 
difference 


=(47/)( tuo m)ax=cx (13) 
where ji, », are the indices of the prism and c is a constant at a given wave-length. 
Let the vectors representing the incident polarized light, in and perpendicular to 


the prism axes of the Babinet, be a,sinwt, a,sin(wt+5). The emergent 
vectors at any point x along the Babinet will be 


a, Sin wt, a, sin (wt + 5+ cx). 


A plane-polarizer, whose transmission plane makes an angle @ to the prism axes, 
will transmit the components of these vibrations in its transmission plane. The 
resultant plane vibration is 


a= 4a, cos @sin wt + a, sin 8 sin (wt +5 + cx). (14) 


The intensity a? of the resultant vibration will vary periodically as x is increased, 
having maxima for (6+cx)=0, 2n7, and minima for 6+cx=(2n+1)m. For 
incident plane polarized light (8=0) maxima and minima will occur for 
cx=0, 2mm, and cx=(2n+1)m. Hence a phase difference 5 between the incident 
components produces a shift 5/27 of the distance between successive fringes. 

Now consider the following experiment. Incident elliptically polarized light 
is transmitted by the Babinet and polarizer, and the phase 5, determined (for 
instance by comparison with plane polarized light). The Babinet and polarizer 
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are now rotated through an angle 7 and the phase 4, is determined, corresponding 
to the vectors in the new axes. Referring now to equation (7) we note that ‘¥ 
has remained unaltered in both measurements, but that we have altered y by r. 
Hence by (7) 

tand,=tan2/sin2y; tand,=tan2'Y/sin 2(y +7). (15) 
Hence we easily derive 

-— tan 6; 
tan 2y =sin 27 l= 5. cos 27) (16) 

and by (7) 


tan 2‘Y = tan, sin 2y. 

The equations (16) determine the required orientation and axis ratio of the ellipse 
of polarization in terms of the measured phases 5,, 5,. In the particular and most 
practical case in which by choice t= 45°, (16) reduces to 

tan 2y = tan6,/tan5,, 7} 

tan 2¥ =tan6, tand,//(tan®S, + tan®8,). | (17) 

Visibility of the fringe system.—Referring to equation (14) we note that the 

ratio of minimum to maximum intensity of the fringe system is 


a, cos 6 —a,sin@\? 
R={— 3 ‘ 
(jeareame) (18) 





To obtain maximum visibility a particular value of @ is required for each state 
of polarization. If however we set @=45°, we obtain a system which will show 
fringes for all configurations except light plane-polarized in or perpendicular to 
the transmission plane of the polarizer. We note that the value of @ does not 
enter into our calculation of the phases 4,, 55. 

Determination of 5 from Zeeman o-component fringes.—In observing a Zeeman 
triplet, we obtain by this technique sets of adjacent fringes in both the 
o-components, as well as in the 7-component of the line observed (see Plate g). 
Now we know that the ellipses of polarization of the o-components are identical 
except for the direction of rotation of the vector. This provides a very convenient 
means of determining 5 without recourse to standardizing fringes. For under 
these conditions, the phase difference of the vector-components in and perpen- 
dicular to the Babinet prism axes are equal in magnitude but opposite in sign in 
the two Zeeman o-components. The phase differences between corresponding 
fringes in the red and violet components are therefore 25, and this may be directly 
measured. (To satisfy oneself that this is the case, we note that in accordance 
with the comment following equation (11), we may represent the ellipses by 
vectors in their principal axes differing in phase by + 7/2 and — 7/2 respectively, 
and by (7) we see that the value of tan2‘Y must also differ only in sign. But 
since tan2‘¥ is independent of y, the phase difference of the red and violet 
components referred to any axes will by (7) always be equal in magnitude and 
of opposite sign.) 

The advantage of this method of measuring 6 lies in the fact that only small 
differences of adjacent fringes photographed simultaneously and under exactly 
similar conditions, have to be determined, and the resulting quantity is twice 
the value required, which reduces accidental error. Another important 
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advantage is that the differential method much diminishes the effect of 
instrumental polarization. This may be seen simply by considering the case 
of a single instrumental metallic reflection, which introduces a spurious phase 
difference between the vectors in and perpendicular to the plane of reflection. 
In the case in which this plane coincides with one of the axes of the Babinet, 
the same systematic error is thereby introduced in the values of 5 of both the 
o Zeeman components, but their difference remains unaltered. There is a 
residual effect in the difference when the plane of reflection is not coincident 
with or normal to the Babinet axes, but this is small in most practical cases. 

The proposed method of measuring the inclination of the magnetic lines of 
force consists, therefore, in setting up before the slit of a large spectrograph, 
a Babinet compensator followed by a plane-polarizer transmitting in a plane 
which makes 45° to the axes of the Babinet. The resulting fringe system in a 
Zeeman triplet is photographed: the Babinet and polarizer are rotated 
conjointly by 45° and the new fringe system photographed. ‘The double phase 
difference 25,, 25, between corresponding fringes in the o-components is then 
measured, and the values of y, tan'Y of the ellipse of polarization is determined 
from equation (17); from tan 'Y, the inclination to the line of sight y, is calculated 
from (12). The values y, yp (respectively the angle measured from the North 
Point in the image plane along the projected direction of the line of force, and 
the angle between the magnetic field and the sight line) completely determine 
the direction of the field at each point in the spot. By a suitable transformation 
the results may be referred to a coordinate system with axes in and perpendicular 
to the solar surface. 

3. Experimental arrangements.—In the summer of 1957 some experiments 
were made at Oxford to explore the practical possibilities of the Babinet fringe 
method ' of field direction measurement. The Oxford 35m solar telescope (21) 
was used in conjunction with the 12m Littrow spectrograph, equipped for this 
purpose with a grating of 20cm aperture ruled with 600 lines/mm by 
H. W. Babcock, and used in the fourth order, where the theoretical resolving 
power, which these gratings closely approach (20), is 480000. The dispersion 
at A6173 is 3°83 mm/A. We are much indebted to Dr H. W. Babcock for 


Fic. 3.—-Apparatus for spatial analysis of sunspot fields: A, objective; B, Cassegrain 
hyperboloid; C, first coudé; D, guiding system; E, first half-wave plate; F, Babinet compensator; 
G, Rochon prism; H, second cmdé; J, re-imaging lens; K, second half-wave plate; L, inter- 
fe rence filter; M, Spectrograph slit. 
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The Oxford 35m solar telescope and spectrograph equipped for the Babinet fringe 
method of sunspot field analysis. 


P. J. Treanor, The spatial analysis of magnetic fields in sunspots 
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6173 A 


Babinet fringes in the a and m-components of the Zeeman triplet 6173, for spots 
12580a (upper photograph) and 12579 (lower photograph). 
P. J. Treanor, The spatial analysis of magnetic fields in sunspots. 
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providing this fine grating, which is ideally suited to work of this kind. The 
arrangement of the accessory equipment is shown in Fig. 3 and Plate 8. The 
35°2m Cassegrain telescope has a hollow polar axis, down which the beam is 
deflected by a coudé mirror. At the lower end of the polar axis, in the 
spectrograph basement, the beam is deflected into an east-west direction by a 
second coudé mirror to the slit of the spectrograph where the image is normally 
formed. In our observations, the focus was advanced to form the solar image 
in the polar beam, before the reflection at the second coudé ; the analysing 
apparatus was placed here, so that the polarization produced by the second 
coudé does not affect the observations. The analysing arrangement consisted 
of a quartz Babinet compensator, of clear aperture 40 x 25 mm with a prism 
angle of 7°30’. The plane polarizer was a Rochon prism 40 x 40mm, of 
fluor-crown glass and calcite. This prism produces an undeflected achromatic 
plane-polarized image, and a chromatic image, polarized orthogonally to the 
first, which is deflected through several degrees and does not enter the spectro- 
graph. ‘The large aperture, high transmission, and optical quality of the Rochon 
gives it a considerable advantage over other polarizers. ‘These components 
were mounted in calibrated rotatable frames on a rigid optical bench supported 
on a frame of ‘‘ Dexion’”’ and erected in the polar axis. A lens, placed after the 
second coudé, re-imaged the spot and fringe system on the slit. 

The guiding system.—It is important to be able to guide the spot on the slit 
accurately during the required exposures, of 5-15 minutes. Since the guiding 
system must be placed in front of the Babinet, it must be as free as possible 
from polarizing effects, a condition not well satisfied by tilted-plate guiding 
devices. ‘The system used consisted of three spectacle prisms. ‘Two of these, 
with their axes arranged orthogonally, could be moved up and down the polar 
optical bench by cords controlled by the observer, who watched the image of the 
spot on the slit through a small telescope. If the powers of the prisms in prism 
dioptres (deflection in cm/m) are D,, D,, and their distances from the image 
plane 2, 2, they produce linear image displacements x= D,2,, y= D2. ‘This 
provides the guiding requirement, but produces an unwanted angular decollima- 
tion of the beam, equal to the vector sum of the angular displacements of the 
prisms. This is compensated for by a fixed third prism of power (D,?+ D,?)** 
placed close to the image plane, its direction of deflection being opposed to 
that of the decollimation. By choosing prism powers in the ratio 3:4:5 
irrational powers were avoided. ‘To avoid vignetting, in re-imaging the Babinet 
on the slit, a field lens is required in the first image plane. The small positive 
power required for this was incorporated in the optician’s prescription for the 
third prism. 

The remaining components shown in Fig. 3 are the interference filter used 
to isolate the region round A6200, and two half-wave plates. The second of 
these, placed in front of the filter plays a simple but important role. The light 
received by the slit is always plane-polarized. The polarizing action of the 
spectrograph is considerable and would nearly extinguish the incident light at 
certain orientations. ‘The half-wave plate is used to rotate the plane of 
polarization to the most favourable position. The first half-wave plate represents 
a slight modification of the experiment as described in the previous section. 
Instead of rotating the Babinet and polarizer conjointly, with respect to the 
incident polarized light, one can rotate the incident polarized light in the 
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opposite direction, the two effects being optically equivalent. A half-wave plate 
may be used for this purpose, the effective rotation of the light being twice the 
angle through which the plate is turned. The action of a half-wave plate in 
connection with the rotation of elliptically polarized light may be understood 
from Fig. 4. In this OX, OY are incident vectors in the plane of the axes of the 
Babinet, Q is the advancing axis of the half-wave plate, at an angle — to the 
OX axis. OX’, OY’ are the transmitted vectors. If OP is the major axis of the 
incident ellipse, making an angle y to the OX axis, the axis of the emergent ellipse 
will be OP’ (the mirror image of OP in the plate axis, as can be seen from the 
symmetry). ‘The action of the half-wave plate has been to rotate the ellipse 
through an angle —2(8+y). The orientation y* of the ellipse OP’ which we 
determine in the experiment in this form is simply related to the orientation of 
the required ellipse by 


y= —(2B+y*). (19) 


| 





é p’ 
Fic. 4.—Rotation of elliptically polarized light by a half-wave plate. 


The advantage of this system is that the orientation of the fringes on the slit, 
and hence their scale in the spectra, remains unchanged. 

Subsidiary measurements.—The values of the retardations and the uniformity 
of the half-wave plates were determined by a method worth mentioning on 
account of its simplicity. A single-element Lyot filter consisting of a piece of 
selected mica, about 1 mm thick, placed between crossed polaroids, was set up 
between a white-light source and the slit of a Hilger direct-reading spectroscope, 
the axis of the mica being oriented at 45° to the transmission planes. The 
spectrum shows under these conditions about half a dozen Lyot fringes. The 
experiment consisted in introducing the wave-plate into the beam between the 
polaroids, so that its retardation (a) supported, (6) diminished the retardation 
of the mica sheet. From the position of the Lyot fringes under these two 
conditions it is an elementary matter to calculate or obtain graphically the 
retardation of the plate as a function of wave-length. We found for the first 
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half-wave plate a retardation of 0-485, at 46173 (it behaves as a true half-wave 
plate at 46330). ‘This slight deviation in specified performance leads to quite 
trivial phase errors never exceeding 5°:5. 

It was originally intended to determine the polarization produced by the 
first coudé at the time of observation. Owing to an unavoidable interruption 
in the observations this was not possible, and a few measures obtained much later 
using artificial light proved rather uncertain. However the order of magnitude 
of the effect can be estimated from the theory of metallic reflection and the known 
optical constants of aluminium. For an angle of incidence of 55°, corresponding 
to a solar declination + 10°, and using the recent values of v and « for evaporated 
aluminium films derived by Schultz (24) we find for rays in and perpendicular to 
the plane of incidence a phase difference « = 20°-6 and a ratio of reflectivities 1-11. 
Drude’s much older values for solid aluminium (4) give a phase difference 24°°3 
and ratio of reflectivities 1-22. Seares (26) has shown that the effect of the 
different reflectivity is negligible ; this is at least true in the case of good evaporated 
films. ‘The phase difference is by no means negligible in itself, but its effect 
enters only slightly into the differential form of the method. It is worth noting 
for future work however, that when, as happens with us, the only significant 
source of instrumental polarization is the phase difference produced by a single 
metallic reflection, this may be accurately compensated by means of a wave-plate 
of retardation « (or 27 +«, a value more easily realized, since such a plate is a 
whole-wave plate for a specifiable wave-length A, A), with its axis in the plane 
of incidence. An alternative possibility would be to set one Babinet axis parallel 
to the plane of incidence (the condition for complete independence of instru- 
mental polarization), and to place the first half-wave plate cemented between 
optical flats, before the reflection at the first coudé. These are alternatives 
to the procedure of determining the instrumental effect, and applying corrections 
to the individual measures. 

4. Observations.—During the months 1957 July-Sept. observations were 
made by means of the equipment described above on a number of sunspots. 
For the field direction determination, the well-known normal triplet of Fe1, 
A6173 was used. Many plates were obtained in which all three components 
of this line were crossed by the polarization fringe pattern. The normal procedure 
was to expose in immediate succession on a single plate, at two different settings 
of the half-wave plate, corresponding to a rotation of the ellipses of 45°. ‘This 
was followed by a third exposure containing standardizing fringes in plane 
polarized light (which were used mainly as a convenient means of obtaining the 
scale of the fringes) and by a fourth exposure in which the Babinet was removed 
and a normal spectrum of the spot obtained for control purposes. The mounting 
of the Oxford telescopes endows the image with a rotation at the diurnal rate. 
By taking the related fringe spectra in immediate succession, it was possible to 
ensure that the associated spectra related to virtually the same section on the spot. 
Advantage of the rotation was taken to obtain different sections of the same spot 
by repeating the observations when possible after a period of some hours. (‘The 
image rotates 15°/hr.) On a few plates, an additional fringe spectrum was 
obtained for an intermediate position of the half-wave plate, corresponding to 
an angle r=22}°. From the collected material, 18 plates were finally selected 
referring to spots in four groups, of which some data are assembled in Table I 


(33). 
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C.M.P. Mt Wilson Mt Wilson 
Aug. 1957 Lat. No. Max H (oer.) 


219 + 16°0 12563 2600 
I'l —30°0 12579 2800 
31°9 +25°0 12580 2300 


319 +25°0 12580a 2300 


All the groups had been observed in previous rotations. ‘The spots studied 
were the principal spots in groups 12563 and 12579. ‘The principal spot in 12580 
was an active elongated object from which attention was eventually transferred 
to a smaller regular neighbouring spot 12580a in the same group. 

‘The method of determining the field direction was essentially that described 
in Section 2. The fringe positions in the o-components were measured on a 
Hilger measuring machine. Settings on individual fringes were usually consistent 
to about o-o1 mm, the distance between successive fringes being 0°33 mm. 
Means of pairs of adjacent fringes were used to reduce the accidental errors. 
The position of the spot centre, which is required to relate the corresponding 
o-values in the pairs of spectra required for each determination of angle, were 
determined from microphotometer tracings across the nearby continuous 
spectrum of the spots. 

Plate gis a reproduction of the fringe systems in A6173 for the spots 12580a 
and 12579, The quality of the spectra was fairly uniform for the four selected 
spots, though these two had a more regular appearance. ‘Two noteworthy 
features in Plate g are the presence of the fringed central component, and the 
extension of the field beyond the limits of the penumbra. Though this latter 
effect is somewhat enhanced by the contrast of the reproduction, comparison of 
the field extent with the penumbral diameters of these spots on copies of the 
daily photographs provided by the Royal Greenwich Observatory, show that 
our field diameters exceed the penumbral diameters by a factor of 1-2. This 
excess may be characteristic of rather small spots in disturbed areas. 

Table II gives some typical results of the determination of the inclination 
of the field to the line of sight, as a function of the distance from the spot centre 
for the same two spots 12579, 12580a. 

‘Two diametral sections of each of these spots was observed which intersected 
at about 45° and 75° respectively for 12579 and 12580a, when they were not far 
from the central meridian (Cols. 1-4). 12579 was again observed when quite 
close to the limb (Col. 5). In Cols. 1-4 the characteristic increase in the angle 
of incidence with distance from the centre is evident, and the agreement at the 
point of intersection near the spot centre is good. The symmetry of the field 
in the isolated spot 12579 (compare Cols. 3 and 4) is remarkable, and also the 
change in the direction of the lines of force with respect to the sight line in the 
limb observations (compare Cols. 4 and 5). ‘The most interesting feature of these 
values, and one which is common to all our results, is the limited range of angle w 
of the observed sections of the cone of lines of force, which in no case approaches 
the value of 180° embraced by the complete field according to Nicholson’s law. 
Since our measures of direction do not extend to the limit of the penumbra, we 
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can only arrive at a quantitative estimate of this concentration of the flux by 
correlating these results with the field strength measurements, which will now 
be described. 


Taare Il 


Variation of field direction for spots 12579, 12580a 
mm 12580a 12579 


—4'0 a = ag 65°58 
— 3°5 ene eee oer 58-3 
—3°O ve go'o 45°6 48°7 
—2°§ 64°5 79°9 43°! 440 
—2:0 57°7 68:9 39°0 43°! 
— 1°s 46°38 611 43°1 44°00 
= Ie 33°9 63°0 33°3 35°9 
—0'5 11‘o 54°85 27°7 231 
roo) 41°3 47'1 27°! 29°4 








+ 0'§ 49°9 318 517 51° 

+ 1:00 49°4 25°9 55°3 44°4 
oe §3°1 

47°1 

oe ove vee 55°3 eee 

w 986 69°5 46'1 74°9 32'0 


+ 1°5 5671 43°! 
+2°0 éae os 
+ 2°5 


Field strengths.—These were determined from the same spectra by direct 
measurement of the distance between the o-components. It was found that 
the grating spectra clearly resolved the Zeeman triplet \6173 at the centre of the 
spots where the fields were of the order of 2000 oersted, and the central fields 
were determined in this way. Attempts to extend the measurements to the edge 
of the spot were hampered by failing resolution of the components, due almost 
entirely to their intrinsic width and to the presence of the central component. 
All the plates however contained the line A6149 (Fett), a line of which the 
physical and magnetic properties mark it out as one of exceptional interest. 


TABLE III 


Intensity 
Source Disk Spot E.P. Transition gA* Pattern 
6149°255 Fe Il 2 —I 3°873 b‘D, .—zP, /.° 50°4 Doublet 
6173°348 Fe | 5 3 2213. a®P,—y'D, 95‘2 Triplet 


Table III1 compiles data from Rowland (32) and Russell and Moore (22) 
relative to this line, and the corresponding comparative data for 46173. The 
line A 6149 is weak in spots, it is an ionized line, and has a high excitation potential. 
These factors all suggest that it is formed in the deeper and therefore hotter part 
of the spot, at a level not very different from that in which the continuous 
spectrum of the spot originates, which according to Houtgast and van Sluiters 
(10) is about 120 km lower than the effective depth of 46173. The Zeeman 
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splitting factor is lower than that of 6173, but the line has a narrower profile. 
Its most interesting feature, however, from the point of view of measurement, 
is that it has no central component. The simple theory of the Zeeman effect 
predicts a quadruplet in which the o- and 7-components have the same splitting 
factor and are therefore superimposed. One would not expect this to be exactly 
verified in the complicated spectrum of Fem. In fact the line was observed in 
the spark spectrum by A. S. King (13) who found it to be quadruplet, being 
made up, however, of very close doublets separated by only 1/10 of the main 
separation. ‘The value of gd? in the table is theoretical. We may compare this 
with King’s value in various ways. Applying the theoretical value to King’s 
measurements, using the mean separation of his 7- and o-components, we find 
that the implied field strength is 16 490 oersted. King estimated his field strength 
in various ways and gives a value of 16 000 + 160 oersteds. The ratio of splitting of 
A 6173 (which King measured with difficulty) and A6149 is found to be 2-049 
which is to be compared with the theoretical value of 1-89. The agreement to 
within a few percent of the theoretical values with these old and not very consistent 
values of King justifies our present use of this value, but if the line is to be further 
used, new laboratory measurements are much to be desired. 

It should be pointed out that the complex polarization of 46149 makes it 
unsuited to the quarter-wave plate technique. This is of no consequence in 
direct measurements however, and it was found that these could be carried out 
over a large part of the spot. The partial failure of resolution near the edge of 
the measurable area produces a systematic error, for which however an order 
of magnitude correction can be applied, on the assumption that the photographic 
profiles of the lines may be approximated to a Gaussian form. The correction so 
deduced rose from about 50 oersteds for a measured field of 1000 oersteds, to 
200 oersteds for a field of 700 oersteds, which represented both theoretically and 
practically the lowest field that could be measured in this simple way. However, 
it was possible to estimate with a consistency of about o-1 mm the point at which 
the field disappeared in both 46173 and A6149. ‘This measurement of the total 
extent of the field provided two additional points of lower weight, which completed 
the field variation curves. 

Although the observations of the four spots were well scattered over the disk, 
analysis of the central intensities did not reveal a significant dependence on this 
quantity, a result to which not much weight can be attached on account of the 
statistical limitation of the material. However, in view of their homogeneity 
in this respect, the observations were used to determine the law of variation of 
field strength over the spots. ‘To do this, standard field profiles were determined 
by expressing the intensity at given fractions r/r, of the field radius r, as a fraction 
I/I, of the central intensity. The results are summarized in Table IV, in which 
they are to be compared with the variation predicted by a parabolic law, differing 
from Broxon’s formula only by the substitution of r,, the actual field extent, 
for the penumbral radius 6. With this adjustment, the agreement is very good. 
Our observations tell us little of the exact form of the curve near the edge of the 
field, and cannot be alleged as positive evidence against Mattig’s formula, but 
they are adequately represented by the simpler parabolic form. 

The law of direction.—Using the above determination of the distribution of 
field strength as a function of distance from the centre of the spot, it is a simple 
matter to calculate the percentage of total flux contained within the limits of the 
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TaBLe IV 
Field-strength distributions in sunspots 


Spot 12579 12580 12580a 
No. of 
spectra 14 9 10 
rir H/H, H/H, H/H, Mean s.d. Parabolic 
of mean. law 
—o'8 37°9 39°6 30°5 34°4 22 36°0 
—o'6 63°7 64°4 56°0 59°5 2°2 64°0 
—Oo'4 84°3 84:2 77°0 789 29 80-0 
—o'2 95°7 97°8 96°6 94°6 19 96°0 
0.0 100°0 100°0 100°0 100°0 oe 100°0 
+0o0°2 99°3 96-2 g2'1 97°0 1°7 96-0 
+0°4 87-0 81°8 73°77 80°7 2°4 80-0 
+0°6 68-2 60°0 55°4 , 61°6 2°3 64:0 
+o°8 41°79 31°8 38°3 37°8 1°7 36-0 


field over which we have determined the variation of direction. In this way 
we find that 40 per cent of the total flux of the spot is funnelled into unit solid 
angle. We can make a similar computation for Nicholson’s measurements by 
combining Chapman’s representation of the latter (equation (6)) with Broxon’s 
law, and in this case find less than 20 per cent of the flux in unit solid angle. It 
is clear that our results indicate a more strongly radially directed field than those 
of Nicholson. 

A partial funnelling of this kind is not altogether unexpected. Indeed, as 
was suggested in equation (5), observations by the neutral line method already 
indicate a degree of concentration intermediate between that implied by 
equation (6) and our present results, while the interpretation of the observations 
based on the extinction of the central component is influenced by a single point 
of low weight near the edge of the penumbra. If the general validity of this 
phenomenon can be established it may well eliminate the inconsistency already 
referred to between the magnetic field measurements and the direction of the 
velocity-fields determined from the Evershed effect. 

The variation with depth.—It has already been stated that no statistically 
significant variation in centre-limb field strengths was found, so that no evidence 
for a variation of field with depth can be deduced in this way from the present 
observations. Comparison of the central fields of 46173 and A6149, yielded, 
somewhat unexpectedly, a systematic negative vertical gradient. However, this 
value is less than the uncertainty of a few per cent in the relative splitting factors 
of the two lines, a quantity which affects all the relative measures. Even allowing 
for the effects of these limited systematic uncertainties, and for the accidental 
errors, it seems however impossible to reconcile the present measures with a 
large positive gradient. An upper limit of 1 oersted/km is suggested. It will 
be recalled that this is consistent with the gradient of 0-5 oersted/km calculated 
from simple geometrical considerations of the divergence required by Nicholson’s 
law of directions. The law of directions indicated by our present observations 
suggests an even smaller vertical gradient. 
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It may be recalled that Tanaka and Tagaki also found only slight evidence of 
a small vertical gradient. Their observations, made with a 14 m solar telescope 
and a dispersion of 2-5 mm/A, refer to a large number of lines, ranging in intensity 
from —3 to 28, in a single large spot. The value of their results has been 
obscured and the interpretation complicated by the large systematic differences 
of the measured fields in lines of different elements, though the agreement between 
lines of the same element is good. The differences must surely be due to the 
use of the splitting values tabulated by Kiess and Meggers (12). ‘These are 
theoretical values derived from Landé’s formula which assumes L—S coupling, 
a condition not strictly applicable to many of the measured lines. The 
recomputation of the homogeneous data provided by these observations, using 
accurate laboratory measurements of the Zeeman splitting of these lines is most 
desirable. 

Conclusion.—Although an attempt has been made to use the first fruits of the 
Babinet method to obtain a pre-view of the complete spatial analysis of a sunspot 
field, it would be premature to put too much weight on these initial and limited 
results, obtained from a random collection of spots, in weather far from ideal. 
The main purpose of the investigation was to demonstrate the practicability of 
the method, and this much at least we believe has been achieved. 

Thanks are due to a number of skilled hands on whom this work has been 
very dependent: to Messrs Hilger and Watts who produced the Babinet 
compensator, and to Mr D. S. Brown of Grubb, Parsons, who constructed 
the Rochon; also to Mr H. P. Wiggins, for several large half- and quarter-wave 
plates used in these and subsidiary experiments, which satisfied the most 
searching tests of quality we could devise, and to Mr H. G. East, who constructed 


various mechanical parts. To Professor Plaskett I owe a special debt of gratitude 
for his unfailing interest and encouragement. I should also thank Dr Hindmarsh 
for valuable theoretical discussions. ‘The project was carried out during the 
tenure of the Skynner Research Fellowship at Balliol College, and the writer 
wishes finally to record his indebtedness to the College. 


University Observatory, 
Oxford: 
1959 July 13. 
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Summary 


The available data for 205 subgiants of type G5 and later are listed in 
Table I. Most of these stars have visual luminosities between +2 and + 4°5 
as judged by spectroscopic criteria. On the basis of current evolutionary 
theory this list probably contains some of the oldest stars in the solar 
neighbourhood. 

The stars for which the most reliable trigonometric parallaxes are 
available, shown in Fig. 1, indicate an (My, B—V)-distribution similar to 
that found for the y Leonis and 61 Cygni groups and the cluster M67. The 
lower limit to the luminosity observed for the subgiants indicates an upper 
limit of ~ 1-5 x 10'° years for their ages. 

The solar motion, derived from the motions of approximately half the 
subgiants, is near 35 km/sec toward a= 19", 6= + 44°. 

The masses derived from the orbits of five visual binaries are discussed. 
First orbits are derived from three of these binaries including HD 194433 
which may contain the oldest stars for which a mass determination is available. 


The masses of the subgiants probably lie between 1 and 1°5 


{ 





1. Introduction.—A class of objects falling in the H—R diagram intermediate 
between the G-type giants and dwarfs was first isolated and labelled ‘‘ subgiants ”’ 
by Stromberg (1930). The two-dimensional spectral classifications at the Mt 
Wilson (Adams et al. 1935), Victoria (Young and Harper 1924), and Yerkes 
(Morgan 1938) observatories, as well as trigonometric parallaxes for a few of the 
brighter members of the class have confirmed the intermediate luminosities of these 
stars. In current theories of stellar evolution the late-type subgiants have been 
recognized as representing an intermediate stage in the evolution of stars with near- 
solar mass. Since evolutionary changes in stars with masses ~1°5© are not 
expected in less than ~ 10° to 10" years, the subgiants are easily recognizable 
examples of old stars. Therefore it is of some interest to collect the available 
information for known subgiants (Section 2) for a discussion of their luminosities 
(Section 3), motions (Section 4) and masses (Section 5). Subgiant components 
of close eclipsing binaries are not included because both the masses and the 
evolutionary rates of such stars may be influenced by their close companion. 

2. Catalogue:—Table 1 contains a catalogue of all stars classified G5 or later 
which, on the basis of spectroscopic luminosity criteria, are suspected to lie between 
M,=+2 and +4°5 (Adams et al. 1935) or classified in luminosity class IV or 
IV-V (MK) and for which radial velocity determinations are available. Stars of 
luminosity class I1J-I1V(MK) are not included. A few stars observed at Mt 
Wilson but not included in the published list of spectroscopic absolute magnitudes 
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have been subsequently classified as subgiants by Joy and so listed in the radial 
velocity catalogue (Wilson 1953). ‘Table I contains the following information: 

HD: The number in the Henry Draper Catalogue. 

Name: The name, BD or CoD number. 

V,,B-V, U-B: The magnitude and colour. The colours are from various 
sources and were all determined photoelectrically. The one decimal magnitudes 
are from a great variety of sources. 

Sp, S: The spectral type and luminosity classification and its source; 
W = Adams et al. 1935; DD = Heard 1956; C= Evans, Menzies and Stoy 1957, 
1959; MK=Johnson and Morgan 1953; F=Boulon, Duflot and Fehrenbach 
1959, R= Roman (1955), MW = M. L. Woods (1955), MP = Moore and Paddock 
(1950), H = Halliday (1955). 

Ha», 5: The components of the proper motion derived from the following 
sources ; GC = Boss General Catalogue ; N30 = Morgan N30 catalogue; Y = Yale 
Zone Catalogues; GAZ=Greenwich astrographic zones; CP = Cape photo- 
graphic catalogues; CR-=relative proper motions from Cape parallax series. A 
plus sign in the column ‘‘S”’ indicates that modern meridian positions have 
been used to improve the catalogue proper motion. 

p,Q: The radial velocity in km/sec from the General Catalogue of Radial 
Velocities (Wilson 1953) and its quality,Q. A few determinations made since the 
compilation of that catalogue have been incorporated into the table. 

Notes: Additional information for several stars is included in the notes; 
alternative classifications by KK =: Keenan and Keller 1953; absolute magnitude 
estimates by V = Young and Harper 1924, and WB = Wilson and Bappu 1957. 

With only a few exceptions the stars in Table I are of spectral classes G5 to K1, 
inclusive. Among all stars known from trigonometrical parallaxes to be within 
20 parsecs of the Sun and classified as G5 to K1, about 10 per cent have lumin- 
osities in the range M,,+2 to +4°5 as can be seen from the following tabulation : 


M, No. (G5-K1) 
<+2 5 
+2to +4°5 13 
> +45 93 


II! 


3. Luminosity.—Mean trigonometric parallaxes greater than 0”-02 based on 
results from at least two observatories are available for 51 of the stars in Table I. 
These stars are listed in Table II where, in addition to the HD number, the 
following information is given: 

m,n: The mean trigonometric parallax (unit =0"-o01) and the number of 
observatories contributing to it. 

B-V: The colour; in the few cases where photoelectric observations are not 
available the colour has been inferred from the spectral type and is followed by a 
colon. 

M,: The absolute visual magnitude resulting from the mean trigonometric 
parallax. 

M,(W): The Mt Wilson spectroscopic absolute magnitude (Adams et al. 
1935) when available. 
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Taste I 
Subgiants of spectral type later than Gs 


HD Name Ver B-V U-B Sp S B, Ms Ss p 


© 


m m « 

28 33 Psc , +1°06 wo “See W -—o013 +0096 GC,Z, N30 — 61 
249 +25°5073 ° rae a Gene DD +0121 ~—o103 Y +15'5 
502 +18°4 . esi we G4IV MP +0018 ~—oos5 G,Y + 2.0 
877 —23°46 . os .. sg@Gs W +0005 +0018 GC,Y —16'0 
1037 —15°32 ° aad .. 8gG8 Ww +o 110 —oo2zg GC, Y —29°8 
1605 +30°43 nel Ki IV H +oo51 +0°057 Y +10°2 
2589 +76°10 , +084 sg Go Ww +0°336 —oo020 GC, GAZ +19°0 
3141 +41°87 ° i KolV MP —o189 —0o°088 GC + oO 
3846 «And +088 sg G3 WwW ~~ O'231 o251 GC, N30 —83°6 
4372 +30°113 ne Ki IV H +O°012 0054 GC,Y +13°0 
4398 HR 210 ‘ akg sg G6 Ww +or19s oo04 GC, Y —14°7 
4744 +29°141 ' +107 G8 IV R +or214 0028 GC,Y — 162°8 
5286 36 And ° 1°03 sg Ki Ww +o'128 —o032 GC,Y + 1°5 
5395 28 Cas ‘ 1°03 i” ae Ww o'092 —or042 GC+ 47°° 
6296 HR 300 0°94 cae. SE Cc O°125 oors GC + 56°5 
6473 HR 312 ° sale sg G6 Ww o'024 ~0'035 GC,GAZ 27°4 
6734 29 Cet 0-87 sg G5 Ww +0133 ~0'418 GC,Z —95"4 
8375 +33°220 a G8 IV H +o226 +0118 GC + 28 
8512 OCet +- 1°00 sg Ko Ww 0078 0216 GC, N30 +17°2 
10015 +28°271 + 0°83 KolV-V R +o'421 —o'257 Y — 2°9 
10486 HR 495 sg K2 WwW +0143 ~oo15 GC +42°2 
11937 x Eri +o8s —< “eee MW +0°676 +0300 GC, N30 - 63 
15590 ~— 42°834 + 0°66 Gs lV ; +0099 —o'065 GC,CP + 36°8 
16141 79 Cet 0°66 Gs IV —~O'1§2 0430 GC,Y -§3'0 
16708 +2°412 a aor sg Gs +0°058 +0°007 Y +65°3 
18907 « For o'79 oo GaSue ; +0°270 0423 GC,Y +39°9 
19618 +14°524 +orB4 ‘st KolIV-V o'069 —0'297 Y —23°0 
19743 —61°167 . + 0°93 a. “eae ; to110 +0106 GC 65 
20277 HR 978 ia .. sgG8 j —O'ols o107 GC +19°0 
20559 95 Cet 5°7 +104 +0" Kr IV +0253 —o'oso GC,Y +278 
20618 59 Ari ‘ ree ee sg Gs : o’021 oo71 GC,Y ol 
20855 43 +105 ieee Ko lV , +0°066 +0102 GC 95 
21467 66 Ari . aa .. sgG6 j 0000 o100 GC,Y +49°0 
22849 +29°592 ; ied Ki IV o'o4! o062 Y + Ba 
23249 4Eri * + 0°92 Ko IV MK 0093 +0°746 GC, N30 6°4 
23817 B Ret . + 1°16 .. KolV ; +0312 +0°077 GC, N30 + 50°5 
26766 +-29°678 be. sas = Kr 1V +0120 o170 GC,Y + go 
29613 HR 1487 5°43 +1°04 oon Os j +O'123 0124 GC,Y +560 
31782 +25°766 3. ad aia KolV 0080 0074 Z,Y -67°8 
32963 +26°789 is +066 +0 Gs lV 0078 0054 Y 62°9 
34538 HR 1737 tog! ae sg Go ; —O°005 oogo GC,Y +75°3 
34642 oCol 0°93 sg Ko ; +0080 0°336 GC, N30 +212 
37160 40 Ori : +O'9O5 sg G6 $+ O°095 o305 GC,Y +98°7 
37601 24 Cam . “ad ee sg Go ; +0020 +0°030 GC, N30 28°9 
37981 +14°991 . ae wo ae —-0'038 oo1o GC,Y +62°7 
38455 —12°1258 ; ai -. Ga +0°009 o12a1 Y -33°0 
39169 —1°1038 : sae «. sg Ko +0012 0044 GC,Z + 30 
39280 — 44°2271 , +O'94 vou ED. ; o'000 +0°030 GC, CP +29°5 
39364 SLep 3 + 1°00 <- @wGF / +0231 o645 GC +99°3 
39861 +81°201 ° ve ws «= WG : —0°O77 0375 GC, GAZ 27°0 
40409 HR 2102 + 1°07 = BKB ; +0135 +0°s40 GC +2g‘1 
40959 +27°963 se we BGS )  —--er00s —orogo GC, Y 2°8 
44019 —O°1255 eee aes K2Iv ] O'013 0024 Y +39°0 
45410 6Lyn +094 +0°66 sgGB8 j —0°020 0°338 GC, N30 +36°0 
46480 8 Lyn aa sg G7 j ~~ ©°200 o279 GC 46°4 


eorrerroraeocrrrsanacrrrsrsg® 


47205 v»CMa 110 as oe +0064 0076 GC,Y + a's 
50806 HR 2576 +0°72 we - ne . +0°273 —0'437 GC,Y +715 


§2071 +27°1296 one os eo 0000 oooo GC, Z +94°7 
52711 HR 2643 . 0°59 ooo  <aeae +0'162 0833 GC,Y +28 
§3590 —0°1587 ven ae ar I +0°013 oo3r Y +15°0 
55280 18 Lyn ; +1°09—Cti«w“a“ kg Kz j —O'Og! e260 GC +23°8 
58895 —58°1741 + 0°70 a re —oo82 +0136 GC +19°7 
59468 —51°2507 7 +070 ... Gs IV-V » —0'291 ooro GC, CP + 3°4 
62644 HR 2998 5 +0°77 ao, - 0072 0563 GC +27°0 
67767 # Cnc . ape sa sg G6 ; 0065 0350 GC, N30 —43'0 
68978 —31°5719 +o°61 .. GsIV-V ; 0336 +0143 GC,CP + 48'5 
70642 --39°4247 +o7o ... Gg IV-V ; 0200 +o240 GC+ +481 
71952 HR 3351 = -. sg Ko F toro1s —oro82 GC +43°6 
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TasLe I—continued 
Subgiants of spectral type later than Gs 
HD Name ‘—e B-V U-B_ Sp Ss Me Hy 


m m ra 

73593 34Lyn 5" AES ~~»  aeG@e Ww +o'029 +0°084 GC 
74543 73°38 , + 1'O4 ne kKolV Cc to'o2z5 +o 152 GC 
76291 HR 3545 5° og +106 Ki IV R +0'126 o'042 GC 
77818 +59°1218 ‘65 +101 +082 sg Ko Ww +o'o12 on Y 
78249 +59°1221 ’ o99 +082 sg K2 Ww —0'163 0033 GC, Y 
79452 HR 3664 5° +086 +037 sgG3 Ww 0147 +0°046 GC 
81265 +31°1978 ° as as | «QE MP +0°075 0183 GC 
82210 24 UMa z o'70 | «6 ~(Gaeary e061 +0075 GC, N30 
84406 +63°861 , me vee sg Ko 1 o'116 0143 GC 
84453 14 LMi , o9os +074 sg Ko j +0°056 0138 GC 
86006 — 45°5556 , +or71 a ee , +0034 +0140 GC 
86359 26 Leo . ois « Gy j O04! oo1g GC,Z 
87783 —46°5806 5° +088 ws eae > o’o10 0056 GC+ 
88416 +27°1860 “ at vid Ko lV 0086 --o'o70 Y 
88532 +29°2012 3 nae we KolV +0038 ~oo016 Y 
90572 +4°2333 ° 7 ww Bae 1 o'116 0028 GC,Y 
giorr 31 Sex . oak ie sg Ko F +0°063 o045 GC,Z 
92588 33 Sex ° ve ara sg K1 | O14! 0124 GC,Y 
93636 43 LMi “« _ «sae j -0°085 0048 GC 
94132 HR 4243 ; han ; sg Go j 0394 ~—o077 GC, GAZ 
94178 +33°2049 . 0 a sg G7 F +oo18 0053 G 
94264 46 LMi 3° hae rm sg K2 ; +0090 0285 GC,N30 
94834 +24 2291 , Bre! mS Ki IV ’ o'058 +0°042 Y 
96074 +66°697 ’ ave - sg G8 j +O°O10 oo1ro GC, GAZ 
96436 65 Leo ° a .. sgG7 4 0°384 oo88 GC,Z 
97033 +66°704 ‘ ak iio ea j 0°336 0132 GC, GAZ 
97561AB + 20°2572 6° +o'75 +033 G7IV 0° 382 0134 GC+,Y 
98818 —60°3502 *57 1°07 ia KolV > o'1ol 0037 GC 
100030 +48 1952 ‘ as we G8 IV —o'230 —o 079 GC 
100901 ‘64 +1° as Thee ; 0148 +0°022 GC 
101408 —45°7180 . +4" sc eee > 0026 oogo GC,CP 
102326 +77°440 "75 + to4s G8SIV-V +o'181 oogo GC, GAZ 
102928 —4°3152 ° 7% sg Ko ’ +0°002 o'oos GC,Y 
108309 —48°7426 Gs IV-V ; 0635 oro8o GC, N30, CP 
111028 33 Vir : sg Ki +0°275 o452 CC 
111417 ~-45 8040 K3 IV ; 0'053 0043 CP 
113817 14°3644 7 wo. e2Gs f -0'042 +001: GC,Y 
114174 +6 2697 oe 2a j +0°o84 0°684 GC, N30 
115202 57 Vir ° en sg K1 1 +0°306 o12z1 GC,Y 
115577 —27°9145 ooo. | eae ; 0074 0084 GC,Y 
117176 70 Vir ow en f 0°238 0582 GC, N30 
121146 HR 5227 ; hae os “oe ; ~O18s oo60 GC + 
121416 HR 5239 fans meee ; —O'159 oo80 GC, CP 
123409 +29 2943 ‘ a) sg G6 j —0'066 o'oog GC, Y 
123797 — 48°8820 t oe) ; —0°078 0078 GC,CP 
125968 ~— 27°9806 «. GsIV-V ; o'264 o252 GC, Y,CR 
128095 +29°2553 veal eas Kr IV +0°056 o11r Y 
129580A +-58°1523 . ae oo. eee +O°131 0196 GC 

B eta -. 8g@G7 F +O'131 0'196 
136136 +44 2444 en ia sg G8 ; 0004 o'e0g GC 
138716 37 Lib ‘ oe §«=6 d +0°305 0240 GC, N30 
140301 HR 5847 ci vee sg Ko J 0005 e100 GC, N30 
140687 —24°12275 7°5 .. sg Ko f —0' 089 e069 GC, Z 
140690 —42°10755 , - Case . -O'116 0158 GC, NP 
142091 «CrB ; oo ~o'006 -—o 355 GC, N30 
142898 +27°2555 : oss K1 IV +0°034 —o'005 Y 
142980 ¢ Ser +1" Ki IV -0'123 —o'087 GC, N30 
143586 +30°2735 one sg G9 F 
144087A —11°4057 . sg G6 ; 0'070 0030 GC,Y 
145148 HR 6014 i KolIV +0230 —0'736 GC, N30 
150275 HR 6191 } sg Ko : —O'101 0276 GC, GAZ 
1s10g0 41 Her , ae .. 8g G6 § o'213 e269 GC 
1§1216 +-33°2775 , ee ve ee hee \e oe 
152781 HR 6284 s oi os - eh ; +0076 +0°034 GC 
153226 HR 6301 “ KolV 00386 «=6+0'°068 GC, N30 
153344 +62°1520 . Gs IV —0'330 —ooss GC 
154088 28°12769 *e -. G8IV-V ; +o°o82 0262 GC,Z 
1sSar1 rs t eee G8 IV > +o 100 ~—0o 006 GC 
168614 HR 6516 .. G8IV-V —o123 ~0172 GC,Y 
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TasLe I—continued 

Subgiants of spectral type later than 

B—V U-B_ Sp Ss H, 

m m oe 
+ 1°02 ee —0'072 
6s ae Gs Ww —0°030 
+0O°75 Gs IV —o308 
seh as: as W +0020 
+1°23 ieee * i —oo16 
+0°65 ote ee —0°037 
at en ae +O°135 
+10! ae +0016 
aay we ee +0°018 
= ia: ene —0°007 
+0°94 ‘72 sg Ko +0113 
+0°94 ‘65 G8IV —O'555 
ive en Sore +0°006 
eas i. aa +0012 
+0°89 me) +0°066 
ans sao.) ee — 0°007 
+075 +0 Ko IV 
+0°97 one: ae 
+o°70 ae) are 
+1°09 rn KolV 
ae KrIv 


HD Name 


gE N 
& 


159809 —45°11742 
160042 —21°4682 
161797 » Her 
162076 HR 6638 
162619 —47°11895 
165271 —46°12117 
165438 HR 6756 
165753 —44°12373 
166070 +27°2960 
166867 +29°3198 
167042 HR 6817 
168723 7 Ser 
170829 HR 6950 
171994 HR 6995 
173791 HR 7065 
173949 HR 7075 
175518 —5°4811 
176578 —47°12692 
177565 HR 7232 
178395 —42°13943 
179094 HR 7275 
181391 26 Aql a eo ow 
182156 5 +1°17 sos. ee 
182572 d . +o°79 +0 sg G7 
183877 ; +0°68 icc’ , ee 
184010 ai “i «-- Gs 
wh ies .. KolV 
185351 5° +0°97 .. sgKo 
186975 Z +1°10 aww eee 
188326 ° +078 +035 G8IV 
188376 ° ais a ees 
188512 . +086 , sg G8 
190248 . +0°76 ube Gs IV 
190360 HR 7670 ° +o72 +0 sg G8 
191026 27 Cyg , +o85 +0 sg G6 
IQIIQO —47 13275 ‘ +116 pa < 4 
191615 +25°4124 : +0°94 
194433 HR 7808 , +098 
196755 « Del " +0°74 
196925 74 Dra 5" saa 
197900 — 44° 14143 , +118 
198149 7 Cep . +0o'91 
198896 +42°3746 8 ica 
199190 —70°1800 , +o°61 
199223 HR 8010 , ‘sii 
199580 +42°3915 : +0°97 
199976 —8°5544 ose 
200026B — 43°14327 6° +0°97 
200553 : +0°93 
206341 —28°17407 . + 1°08 
208627 —44°14717 6 +089 
210441 —44°14805 , +1°00 
ae +15°4622 ees 
212943 35 Peg ' + 1°05 
215549 +29°4753 ‘37. +0°92 
218101 HR 8784 ° +0°83 
218935 60 Peg ies 
219834 94 Aqr 
220167 +59°2703 
221148 HR 8924 
221585 +62°2244 
221639 +59°2744 
222107 AAnd 
222404 yCep 
222480 —32°17623 
222803 HR 8993 
223524 HR 9027 
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Notes to Table I 
HD HD 
28 Kr III (4550) R, Sp. B. P=739, V (41-7) 142980 gK4(—0o-1) W, WB (+2:0) 
G8 III R, V (+1°4) 144087 AD 10 
ADS755, WB (+40) 145148 a itso W, WB (+3°7), Kr II-IV 
G8 ITI-IV (wk lines), R, V (+0°4) 
Ko IVR 151090 omy K3, 164” common motion 
dG8 (+5:0) W 150275 Kr i R, WB (+2:0) 
11™ 6” 153344 dG4W 
dG2w 158614 dG6(+4°8) W, ADS10598 
ADS2459, gKo (+0°9) W 159809 SDS 
dKo (+4:6) W, WB (+5:1) 161797. dG4(+42) W 
G2 V DD 162076 WB (+2:1) 
Ko III R, G6 III-IV KK 165438 WB (+3°6) 
ADS4442 167042 WB(+4'1), Kr IIR 
G8 III KK 168723 Ko III-IV MKK 
Optical comp. to ADS4629 170829 = Sp. B. 
Ko IVR 173949 ADS11661 
; 175518 dG8(+5-3) W 
179094 Sp. B. P=294 
181391 Sp. B.P= 2664, G8 III-IV R, WB (+1°3) 
182572 G8IVR 
+19°4061 Vel. by Fehrenbach 
185351 Ko III R, WB (+2°3) 
ty! C+ 3-7) Ww 188326 dG4(+46) W 
Ko IVR 188512 G8 IV MKK, ADS13110 
The values from CP are 4,= +0057, 190360 G6 IVR, 14™-8 M6 2’-5 
ig= +0°'030; these are supported by 191026 KoIVR 
relative values determined from parallax 191615 Ko III DD 
plates. 194433 K2IV-VC,SDS 
SDS, o” 3 196755 dG2(+3°8)W 
ADS777 196925 9™@ Go 3’, Sp. B range 16 km/sec 
K1 Ill R WB (+2:1) 198149 KoIV MKK, WB (+2:9) 
ADS8060, Am~6™ 189223 ADS14430 
dG4 W; ‘ADS8094, Am~o™ 198896 Sp. B with range of 36 km/sec, optical 
dG7 W companion 
Sp. B., range 28 km/sec. 199580 KrIVR 
V=10"-4, B—V=+1™-28, common 199976 ADS14528 
proper motion 200026 SDS 
K1 IV R, Ki, HI-IV KK, WB (+3°7) 212943 Ko III R, WB (+3:2) 
K1 IV (4550) R 215549 Ko III-IV R, ADS16248 
Gs IV-V MKK 219834 ADS16672, dG4 (+36) W 
ADS9346 220167 ~+Vel. by Fehrenbach 
ADS9573, optical 221148 K3ITIIR 
WB (+3°0), Kir III R 221585 dG3 (+44) W 
SDS 073 222107. Sp. B. period=214, G8 III-IV MK 
Ko III-IV R 222404 K1 IV MKK, WB (+2°6) 


U, V, W: The vectors of the space motion directed, respectively, away from 
the galactic centre in the direction of galactic rotation and toward the north 
galactic pole. The trigonometric parallaxes together with the proper motions and 
radial velocities in Table I were used in computing these motions. 

The colours and absolute visual magnitudes from the trigonometric parallaxes 
are plotted in Fig. 1 where stars with (1) 7 >0"-05 are indicated by filled circles 
(2) 7<0”-05 as determined at more than 2 observatories are indicated by open 
circles and (3) 7<0”-o5 from observations at only two observatories are indicated 
by crosses. The standard main sequence shown in the figure is that of the Hyades 
cluster (Sandage and Eggen 1959). 

Also shown in Fig. 1 are the subgiant sequences of a globular cluster (M13), 
two high velocity groups (y Leonis and 61 Cygni) and one galactic cluster (M67). 
The only available photometric results for M67 are those of Johnson and Sandage 
(1955). From photoelectric (U, B, V)-photometry of about 50 stars in the cluster 
Johnson and Sandage found an ultra-violet excess which they interpreted as due 
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Stars in Table I with trigonometric parallaxes at least as large as 0"-02 from determinations 
at 2 or more observatories 


HD m, n B-V My My(W) U Ww 
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3546 
5395 
8512 
11937 
16141 
18907 
23249 
23817 
37160 
39304 
40409 
45410 
47205 
50806 
62644 
67767 
94132 
96436 
108309 
111028 
114174 
115202 
117176 
138716 
142091 
144087 
145148 
151090 
#158614 
161797 
167042 
168723 
170829 
175518 
177565 
181391 
182572 
183877 
188512 
190248 
190360 
191026 
194433 
196755 
198149 
199190 
212943 
215549 
219834 
222107 
222404 


31 
40 
34 
52 
33 
35 
109 
42 
a4 
22 
40 
23 
52 
35 
50 
29 
24 
33 
36 
21 
42 
4! 
41 
24 
36 
39 
26 
27 


4 
3 
4 
2 
2 
2 
4 
2 
4 
3 
2 
2 
2 
2 
2 
2 
2 
3 
3 
4 
2 
2 
2 
3 
2 
3 
5 
3 
6 
2 
2 
4 
2 
3 
2 
2 
3 
2 
5 
2 
2 
4 
2 
4 
4 
2 
3 
4 
2 
3 
4 


m 
+088 
+ 1°03 
+ 1°00 
+o°85 
+0°66 
+o°79 
+0°92 
+1°16 
+0°95 
+ 1°00 
+ 1°07 
+0'94 
+ 1°00 
+0°72 
+0°77 
+0°7: 
+o°8: 
+0°7: 
+ 0°67 
+0'99 
+ 0°66 
+10: 
+o°71 
+101 
+ 1°03 
+0°74 
+ 1°02 
+oO'7: 
+0°72 
+0°75 
+ 0°94 
+0°94 
+o'8: 
+0°75 
+0°70 
+0'9: 
+079 
+0°68 
+0°86 
+0°76 
+072 
+085 
+098 
+0°74 
+o°91 
+0o°61 
+ 1°05 
+0°92 
+0o'7: 
+1°02 
+ 1°03 


m 
+1°8 
+2°6 
+12 
+2'1 
+4°4 
+36 
+3°7 
+1°9 
+10 
+04 
+2°6 
+2°7 
+2°4 
3°7 
+3°6 
+3°2 
+3°0 
+3°3 
+4°0 
+2°3 
+5°0 
+3°2 
+31 
+1°5 
+2'4 
+5°5 
+3°1 
+3°9 
+4°5 
+3°6 
+24 
+1°9 
+38 
+4°9 
+5°2 
+2°3 
+4°0 
+4°0 
+28 
+48 
+4°0 
+2°9 
+3°9 
+ i497 
+2°7 
+4°3 
+1°4 
+3°4 
+2'8 
+2°0 
+2°2 


m 
+2°6 
+2'8 
+ 2°! 


+4°1 
+4°6 


+2°5 
+31 


+2°5 
+0o°7 
+4°5 
+4°2 
+4°4 
+26 


+2°4 
+4°2 
+2°1 
+4°4 
+2°1 

+2°0 
+43 
+5°6 
+4°5 
+48 
+4°2 


+39 
+5°3 
+2°3 
+4'5 
+4°7 
+4°0 
+4°6 
+4°5 
+2°5 
+38 
+26 


+2°3 
+4'8 
+3°6 
+2°3 
+2°! 


+ 


Heh tee i tei titi 


+++ 


' + 


15 

3 
28 
11 

9 
20 
13 
16 
33 


28 
11 
2 
6 
35 
45 
22 
27 
II 
18 
10 
9 
8 


35 
18 
13 
88 

6 
19 

5 
22 
15 
32 

8 
17 
12 
21 
46 

2 
89 

2 
10 

4 
12 
19 
82 
31 
10 


9 
I 


*158614 Luminosity corrected for equal components before plotting in Fig. 1. 
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to a reddening of E, y= +0™06. After correcting for this reddening, they 

fitted the main sequence of M67 to that of nearby stars and found a corrected 

modulus of V,—M=9™.5. From a more refined analysis of the data, Sandage 
(1957) later altered the modulus to 9™-2. 

My 

-2r 








L i 


+14 








Fic. 1.—The (My, B—V)-diagram for the stars in Table II. The subgiant sequences of 


the disk-population groups, y Leo and 61 Cyg, and cluster M67 and the halo population globular 
cluster M13, are shown schematically. 


The unevolved main sequence of M67 begins at apparent magnitude V, = 14™. 
Photoelectric observations are available for the following 11 cluster members 
fainter than V, = 14™ (Johnson and Sandage 1955, Table 1): 


Star Vez (B—Vy,) (B-V) Mod(1) (B-—V)* Mod(2) 


m m m m 
14°95 +0°69 +0°75 9°35 +0°63 10°03 


1498 14°91 0°78 0°84 : o"72 9°48 
i 15°23 0°84 0-90 ; 0-78 9°48 
14°83 0-69 0-75 , 0°63 9°91 
16°19 0°66 o"72 , o'60 (11°41) 

14°05 0°64 o’70 , o's58 9°39 

15°78 0°84 090 F 0-78 10°03 

15°64 089 0°95 ‘ 0°83 9°64 

16°05 1°07 1°13 A 1'o1 9°25 

132a 15°83 1°08 I'14 1°02 8-94 
149b 14°76 +0°48 +0°54 (10°28) +0°42 (11°20) 


got+o™3 (AD) 9°6 + 0°3 (AD) 


The photoelectric (U, B, V)-photometry of 34 cluster stars with 
(B—V),4,> +0™-4 indicates a mean ultra-violet excess of 5(U—B)= +0™-06 
from the standard (U—B, B—V)-relation for the Hyades main sequence stars 


32" 
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(Sandage and Eggen 1959). The A(B~V)« F[8(U — B)] relation previously 
derived then indicates that if the excess is duc to line-blanketing, the observed 
values of B~ V must be corrected by + 0-06 before the main sequence stars can be 
fitted to the Hyades main sequence. The values of (B— V)! listed above result 
from this correction and the values of modulus (1) represent V,— My (Hyades) 
for the corrected colours, However, if we interpret the observed excess as a 
reddening effect, following Johnson and Sandage, we obtain the corrected values 
of (B~ V)* and of modulus (2), ‘There is no strong evidence to support either 
of these two extreme interpretations of the observed ultra-violet excess, ‘T'wo or 
three early-type cluster stars have been observed but these are cither (1) field 
stars, (2) unevolved main-sequence cluster members, in which case their presence 
is unexplained by current theory unless they were formed much later than the 
majority of the cluster stars, or (3) cluster members which originated as F-G 
dwarfs and whose evolutionary development has carried them to their present 
position in the (M,, B—V)-diagram. Although either alternative (2) or (3) is 
probably the true situation, neither carries any confidence that the colours of these 
stars can be directly compared with unevolved A-B stars, even if we knew the 
appropriate (U — B, B= V)-relation to which they should be referred. 





A 











gaye 40 
Pia, a=The (UB, BV «relation for the few stars in Table 1 with available (U, B, V)« 
photometry, The full curve represents the standard relation (Hyades) and the broken line that 
for M13. 

The sub-giant sequence of M67, obtained from modulus (1) is indicated in 
Fig. 1 together with the subgiant sequence for the high velocity groups y Leonis 
(Eggen t9sga, d) and 61 Cygni (Eggen 1959 b) and the globular cluster M14 
(Eggen and Sandage 1959). Lf we used modulus (2) for M67 the tip of the sub- 
giant sequence in Fig, 1 would join the 61 Cygni sequence. It should be 
emphasized that all the subgiant sequences shown in Fig. 1 are highly schematic. 

The (U ~ B, B= V)-diagram for the stars in Table | for which three-colour 
photoelectric observations are available is shown in Fig. 2 where the standard 
relation for the Hyades (Sandage and Eggen 1959) is shown as the continuous 
curve and that for M13 and the very high velocity group, Groombridge 1830, asa 
broken line, If we accept the interpretation that the observed excesses of stars 
are caused by line-blanketing and if the stars in Fig. 2 are representative of those in 
Table I, then the subgiants represent (1) a few stars with ultra-violet deficiency 
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and therefore stronger metal lines than the Hyades standards, (2) a large number 
(more than 50 per cent of the total) with an observed excess between o and o™1, 
similar to the excesses observed in the y Leonis and 61 Cygni groups and in M67, 
and (3) a few stars with excesses as large (~ +02) as those observed in some 
globular clusters. Because the number of stars with available (U, B, V)-photo- 
metry is small, and because there exists the possibility of contamination by 
unknown companions, it can not now be determined whether there is a continuous 
variation from 6(U ~ B) «0 to the observed extreme of +025, That is to say, if 
we adopt the nomenclature proposed by the Vatican conference on stellar popula- 
tions (O'Connell 1958), it can not now be decided whether or not the subgiants in 
the solar neighbourhood represent a continuum from ‘‘ disk '’-population stars 
with high or moderately high metal abundances, represented by M67 and the high- 
velocity groups, to ‘‘halo’'-population stars with very low metal abundances, 
represented by globular clusters such as M13 and the very high velocity groups 
(Groombridge 1830). This is a question of some importance that might be 
answerable when (U, B, V)-photometry is available for more of the stars in Table I. 
However, because of the very high dwarf/subgiant ratio in globular clusters and 
the relatively few globular cluster dwarfs in the vicinity of the Sun (Eggen and 
Sandage 1959) we would expect very few globular cluster (halo) subgiants among 
the objects in Table I, 

‘The lower limit to the luminosities of the subgiants seems to be well defined 
in Fig. 1 and corresponds to the luminosity of the M67 subgiants if modulus (1) is 
correct. This lower boundary to the luminosity of the subgiants indicates an 
upper limit of ~ 15 x 10" years as the age of the disk-population near the Sun, 
using Hoyie’s (1959) dating procedure. 

‘The marked variation in slope shown by the subgiant sequences illustrated in 
Fig. 1 seems to be unrelated to the ultra-violet excesses which range from +025 
for M13 through + 0-04 for y Leonis, 40-07 for 61 Cygni group and — 0-02 for 
the very similar { Hercules group, (Section 4, Table I11) too or + 0-06, depending 
upon the amount of reddening, for M67. At present no group of disk-population 
stars with a mean ultra-violet excess between o1 and o™2 is known, 

4. Motion.The components of the space motions listed in ‘Table II are 
plotted in Figs. 3(U, V) and 4(V, W). The ellipse shown in Fig. 3 is that found 
previously (Eggen 1959 c) to contain all the known Ap and Am stars brighter than 
V,~6™ and probably represents the (U, V)-limits for the ‘‘ older population 1”’ 
stars (O'Connell 1958). The lack of subgiants falling within this ellipse in Fig. 3 
is at least partly caused by observational selection because parallax programmes 
mainly contain stars with appreciable proper motion, Space motions have also 
been computed for an additional 63 stars not in ‘Table I] but for which values of 
M,(W) determined at Mt Wilson (Adams et al, 1935) are available. Six, or ten 
per cent, of these stars have values of (U, V) that fall within the ellipse drawn in 
Fig. 3. ‘The (1) median and (2) mean values of U, V and W for the §1 stars in 
Table I] plus the 63 additional stars for which values of M,(W) determined at 
Mt Wilson are available, together with the (3) median and (4) mean values 
previously found for the Ap and Am stars are as follows : 

U V Ww 
(1) +8 ~ a’ —8 km/sec 
(a) +9°3 ~33'9 “94 
(3) +7 - 6 8 ” 
(4) +7'8 ee ii . TS @ 
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The median and mean values for the subgiants give the same apex for the solar 
motion (a= 19", = + 44°) but the total solar motion from the median values is 
30 whereas the mean values give 36km/sec. The average motions of the subgiants, 
which are largely a disk population, differ from that of the older population 1 stars 
only in the value of V. This merely means that the disk-population contains stars 
whose galactic orbits have a larger range in semi-major axes (a) than the stars 
contained in the older population 1, but both populations have the same upper limit, 
near a= 1°1 RO, where RO is the Sun’s distance from the galactic centre. 
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Fic. 3.—The (U, V)-vectors for stars in Table II. The ellipse encloses the (U, V)-vectors of 
all the known Ap and Am stars (older population I) brighter than Vg=6™. 
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Fic. 4.—The (V, W)-vectors for the stars in Fig. 3. 


Several groups of high velocity stars, which contain some of the subgiants in 
Table I, have previously been discussed. ‘The largest are the 61 Cygni group 
(Eggen 1959b) the y Leonis group (Eggen 1959 a, d) and the ¢ Herculis group 
(Eggen 1958). The stars now considered to be probable members of the later 
group are listed in Table II1-and the motions of the stars are in Table IV. The 


subgiant sequence of the { Herculis groupis nearly identical to that of the 61 Cygni 
group in Fig. r. 
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There are undoubtedly additional groups among the stars in Table I. How- 
ever, if the very large ratio of dwarfs to subgiants that is observed in M67 also holds 
in the groups, it is more profitable to search for moving groups among the dwarfs, 
and observational programmes are now being carried out to increase the material 
available for suchasearch. An example of a possible group of subgiants in Table I 
which, if the dwarf/subgiant ratio is as high as that in M67, should contain very 
many dwarfs, is the following: 


HD My U V W Remarks 


m 

+1°48 —62 —11 £B Ret; 7=32Y(16), 67C(7) 
+3°2: — 63 774 = 12A(20) 

+ 5°05 —64 - 6 7, = §6A(16), 61M(7) 
+2°3: — 67 m™ = 22G(8) 

+284 —68 — 9 


Tase III 


Probable members of the { Herculis group 
HD Name Veg B-V U-B XU-B) Sp. 
m m m m 
2151 B Hyi 2°77 +0°61 a toc. aa 
4075 +75° 36 7°4: +07: Ae ves ie 
9166 +67° 133 6°76 +123. +1°39 ae | ae 
19058 p Per 3°5 var +1°70 te .. M4II-Ill 
25728 « Ret 4°92 +1°45 bis in - Oe 
43899 HR 2263 5°4: +11: ay gKo 
68788 +73°407 8-38 +084 +0°53 Ki V 
70523 —17° 2464 5°77 +105 +089 Kr III 
71377 —12° 2524 5°54 +119 +1°14 . Kell 
83425 +5°2207 4°66 +130 «64+1°46 K3 Ill 
89668 —0° 2726 9°41 +109 4+1°00 K3 V 
90250 29L Mi 6°49 +109 §=6++0°99 Ki Ill 
98824 71 Leo 7°04 +105 +0°96 K1 Ill 
100733 HR 4463 57 var +1°66 oie oss eee 
119425 84 Vir 5°37 +110 6©+1°05 —0o'06 Kr Ill 
129245 +80°448 6:26 +1300 «0 +1°46 .- Ko HI 
136422 ¢ Lup 3°59 +1°55 “5 ... gMo 
*150680 ( Her 2°82 +064 +021 -—0'03 GolIV 
166207 +50° 2525 6:1: +10: a a ay 
210967 « Oct 48var +1°6: M6 
217166 +8°4973 6°64 +065 ai vee GE 
*219829 +4°49904 8:03 +080 +0°43 ooo Ko V 
221354 +56° 2605 6°75 +083 +0°50 ooo Ko V 


mean= —0-02 
*150680 ADS 10157 
219829 ADS 16645, Am=1™, P=108 yrs. 


The individual values of the absolute trigonometric parallax listed under 
‘* Remarks ’’ are in units of 0”-oo1 and are followed by the observatory abbrevia- 
tion and, in parenthesis, the weight. If these stars belong to a moving group they 
are among the faintest, and therefore oldest, subgiants known since they indicate a 
subgiant sequence similar to that obtained for M67 from modulus (1) in Section 3 
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Taste IV 


Motions of the probable members of the [ Herculis group 

p Ha Me S U V Ww 

km/sec x ” km/sec 
2151 +23°0a +2°228 +0324 GC, N30 +55 — 44 —28 
40°75 — 84b +0°380 —oog5 GC+ +59 —42 —20 
9166 —13'9b +0°104 —oo52 GC+ +53 — 48 —30 
19058 +28-2b +0°131 —o'106 GC, N30 +57 —47 — 16 
25728 +60°5 a +0°067 +o'092 GC, N30 +51 — 52 — 30 
43899 +70 d — 0007 +0085 GC, CP +54 — 50 —20 
68788 — 34b — 0°300 —0'483 GC+ +56 — 50 —22 
70523 +67°7b —0o'100 -—oo1o GC,Y +56 —53 —25 
71377 +65 c — 0085 —o'024 GC,Y +54 — 53 —25 
83425 +45'24a —o'164 —o056 GC, N30 +60 —48 —20 
89668 +35 Cc —0°655 —o165 Y, CR, Ci +59 —49 —18 
go250 +11'4b —O'r110 -—o071 GC +50 —50 —21 
98824 + 5 3b —0O*120 -—oo20 GC+ +59 —43 —20 
100733 +184b — 0°087 -—oors GC, CP +54 — 54 —16 
119425 —42°4b — 0293 -—o077. GC +59 —44 — 28 
129245 —22°6b —o'100 +o°095 GC+ +59 —4!1 —25 
136422 —29°'4a — 0°093 —o'086 GC, N30 +53 —48 —18 
150680 —69°9 a —o'471 +0390 GC, N30 +54 —47 —24 
166207 —57'1b —0°006 +o:100 GC, N30 +56 —44 —17 
210967 +117 4 +0°073 —o-o040 GC, N30 +56 —47 —19 
oo —26°5b +0°392 —o142 GC+ +50 —54 —23 
19829 ~—14'3b +0°466 —o108 GC+ +57 —47 —21 
@21354 —25‘1b + 1°095 +o111 GC +59 —51 —12 


HD 


and shown in Fig. 1. Also, it is of interest to note that the short period (04-06) 
variable CY Aqr may also be a member: 


CY Aqr; #g=0"'059, pg= —0"'065, p= —32 km/sec. 
7 My U V W 

. m 
0°005 4°5 +26 —81 — 30 km/sec 
0006 48 +24 —71 —21 
0°007 51 +21 — 64 —14 
0°008 5°4 +20 — 54 —- 9 


” 
” 


There is no direct determination of the distance of CY Aqr but the value M, = 5™ 
derived from the assumption that it shares the space motion of 8 Ret is about a 
magnitude fainter than that derived from the trigonometric parallaxes of two 
similar stars, SX Phe (04-06) and AI Vel (04-11), for which M, = + 4™. 

5. Masses.—The incidence of binaries among the stars in Table I is not 
high. Only 5 per cent of the stars are known spectroscopic binaries and an addi- 
tional 10 per cent have known visual, mostly widely separated, companions. 
Orbits have been published for only two of the visual binaries, ADS 755 (36 And) 
and 10598 (HD 158614). Preliminary orbits have now been computed for an 
additional three stars, 95 Cet (ADS 2459), HD 87783 (SDS), and HD 194433 
(SDS). The five binaries are discussed below. 

(i) 36 And(ADS 755). The period is 166-5 years (Baize 1946) and the semi- 
major axis, adjusted to the ‘‘standard observers’’ (Eggen 1956) is a=0"-965. 
The magnitudes are nearly equal with the mean component having V,=6-2, 
B—V=+1™03. The spectral type is sgK1. The system was discussed 
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previously (Eggen 1957) when it was shown that the parallax is probably in the 
following range: 


My(AB) m(AB) 
° m 
0°020 +2°7 20 © 
0°023 +3°0 13 
0°025 +3°2 1°o 


where AB indicates the mean component. There is no trigonometric parallax 
available but the spectroscopic absolute magnitude is M, = + 2™-3 (Adams et al. 
1935). Wilson and Bappu (1957) found M, = +4™-o from the widths of the 
emission H- and K-lines; this value would give 7=0"-036 and the improbably 


low mean mass of 0°35. The most likely values are 7 = 0-023 and m(AB) = 1-30. 
(ii) 95 Cet (ADS 2459). I have derived the following preliminary orbit: 


P= 200 yrs. A= +0°393 
a=1°030 B= +0'410 
e=0'48 F=+0'162 
T= 1841.1 G= +0887 


The comparison with the more reliable observations is as follows: 


t PA Sep n Obs o-C 
1854.81 7 i 0-70 + 2 Da z +0°03 + 
1888.77 118 0°45 + 0°02 
1898.15 137 O'51 +O°14 
1899.81 157 0°35 — 0°04 
1900.78 st aa 
1907.14 172 0°40 
1916.78 198 "50 
1921.69 211 o'52 
1929.79 218 o’72 
1935-34 220 0°76 
1942.15 224 o'80 
1948.99 233 0-92 
1956.11 236 0°96 
1957-73 237 1°08 


pee 


+0°02 
+o°o! 
—O°'0s§ 
+ 0°04 
—o’ol 
— 006 
+o°o! 
—0'07 
+0°04 


< 
oo 


<w 
je) 


wears On NPP eH DN 
N 


Dew err ry 


The total magnitude and colour are V, = 5™-37, B— V = + 1™-04 and the spectral 
type is K1 IV. The components differ visually by 2™-o. Values of the absolute 
magnitudes of the individual components, the total mass of the system, Xm, and 
the mass of the secondary component computed from the mass-luminosity relation 
(Eggen 1956) on the assumption that it is a main sequence star, m(B),, are listed 
for four values of the parallax; the difference Xm — m(B), is taken to be the mass 
of the subgiant, m(A). 


M,(A) My(B) =m m(B), m(A) 
° m m 

O-Ors 1°3 3°3 810 160 650 

0-018 1°7 3°7 4°7 1*4 3°3 

0°020 1°9 3°9 3°4 13 21 

0°022 2°1 4°1 26 12 14 
If the subgiant has a mass of 1-40, the second component is ~ F6 and the parallax 
is 0"-022; these seem to be the most probable values. The trigonometric parallax 
determinations, in units of 0’-o01, are 16A(20), 2M(10) and — 17Y(10). 
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(iii) HD 87783 (SDS). 
by D. Jones (1959): 


Olin 7. Eggen 
The following preliminary orbit has been computed 


Vol. 120 


P=232 years 
a=0'64 
e=0°72 
T= 1932.8 


A= —0"607 
=—O211 
F= o165 
= —0°495 


The comparison with the observations is as follows: 


t 


1897.10 
1901.70 
11.27 
14.24 
26.40 
39.34 
30.43 
35-97 
45.29 
47.33 
48.54 
51.30 
56.25 


PA Sep 
o6+ 
0°53 
ost 
0°20 
o15? 
O14 
O14 
0°30 
0°30 
o°31 
o-41 
0°47 


n Obs 


seNPNK SNW HR SS Me 


wawwwawdenmonmeen 


N 
wes OW KH OCNWH OKO 


Lat + 


+ 


os 


a 


Oo 


+ 


wv 


O-C 
oot 
— 0°05 
+o1rt 


— 0°02 
—0°03? 
— 0°04 
— 006 
o"00 
— 0°03 
—0°04 
+0°02 
+o°o!I 


The total magnitude and colour are Vp = 5™-07, B— V = +0™-88 and the spectral 
type is Ko IV. The components differ visually by 1™-5. Three assumed values 
of the parallax give the following results : 


My(A) My(B) 
’ m m 
o'o10 o°3 18 
ool! o's 2°0 
o-o12 0-7 22 


im m(B), m(A) 
490 2°50 2°40 
3°7 2°4 1°3 
28 2°2 o"4 


If the preliminary elements are nearly correct, the parallax is almost certainly in 
the range shown above and the luminosity of the bright component is rather high 
for asubgiant. It is possible that the parallax is o”-o12 and both components are 
subgiants with equal mass, 1-40. 

(iv) HD 158614 (ADS 10598). The period is 46-1 years (Ashbrook and 
Duncombe 1951) and the semi-major axis, adjusted to the standard observers, is 
a=0"-g2. The mean magnitude and colour of the nearly equal components are 
V,=6™0, B—V=+0™-72. The absolute magnitudes and masses of the mean 
component for three values of the parallax are as follows: 

My(AB) 


m(AB) m( AB), 


oO 
0-050 
0'055 
0060 


m 
+4°5 1°40 10 
+4°7 I'l ere) 
+49 0°85 0°95 


where m(AB), is computed from the mass-luminosity relation for main sequence 
stars. The trigonometric parallax determinations are 55A(20), 31M(8), 58Y(10), 


67Yk(4), 495(7), 39C(8), mean = 50. 


The stars are probably only ~o™-5 above 
the main sequence with m= 0"-05. 
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(v) HD 194433(SDS). The following orbital elements have been derived by 
D. Jones (1959); we are indebted to W. H. van den Bos for the 1959 observations : 


P=135 years A=+0°214 
a=o'79 B= —0°665 
e€=0'59 F=-—0-'750 
T= 1949.8 G=—-—o'125 


The comparison with the observations is as follows: 


t : Sep n Obs o-C 
” ° 
1880.86 97 

9°.73 96 1-08 
96.67 100 0-78 
98.68 93 1°44 
1902.53 83 1°31 
13.62 71 0°93 
16.69 71 1°02 
25-53 54 1°00 
27.22 50 0°66 
28.68 43 o'72 
29.42 44 o'71 
30.62 40 0°66 
33-70 32 0°67 
33-88 33 061 
36.40 23 0°50 
36.69 23 0°50 
41.07 3 0°42 
42.36 358 0°44 
43-48 338 0°37 
44-75 340 0°34 
45.46 330 0°36 
45-75 330 0-38 
47-70 317 0°34 
49.79 292 028 
52.70 245 0-25 
59-57 188 o°52 


— 0°03 

o'0°0 
—0°23 
+0°31 
+0°35 
+0°04 

o'17 
+0°26 
—O°05 
+0°04 
+ 0°04 
+o'o! 
+0°06 

0°00 
— 0°05 
a O04 
—O'0s 

0"00 
—0O°O§ 
— 0°05 
+ 0°03 
+0°04 
+0°02 
—o'o! 
—0°05 
+0°02 


The total magnitude and colour are V, =6™-23, B— V = + o™-98 and the spectral 


type is sg K1, Mt Wilson or Kr IV-V, Cape. The components differ visually 
by 1™17, Four assumed values of the parallax give the following results : 


My(A) My(B) Lm m(B)e m(A) 
» m m 

0°020 +3°1 +4°2 120 220 

o’o2I +3°2 +43 : 1°2 "7 

0°022 +33 +44 : 1‘t 14 

0°023 +3°4 +4°5 jf I'l I'l 


IHD NKWNFWHWSEU DHE OBE NNWNNN HWM 
See eee lls oe 


If the orbital elements are correct the mass of the brighter component can be 
quite accurately determined. If M,(A)< +3™2 the derived mass is greater 
than that of a main sequence star of the same luminosity. Although this result 
depends on the assumption that the fainter star is a main sequence object, the 
fact that the magnitude difference between the components is 1™-2 makes the 
assumption almost certainly true. On the other hand, if the luminosity of the 
subgiant is as faint as M,(A)= + 3°6 the mass is rather lower than expected for 
an evolved main sequence star. We will adopt 7=0-023, m(A)=11©. The 
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trigonometric parallax determinations are 25C(10) and 43Y(8). This important 
system deserves further study since it may contain the oldest stars for which a 
mass is now obtainable. 


It seems well established from the binaries discussed above as well as from a few 
others outside the range of luminosity and spectral type discussed here ({ Her, 
y Leo) that the masses of the late-type subgiants are in the range 1 to 1°5©. 

In addition to the 5 pairs discussed above there are only two others in Table I 


that may offer an opportunity for a mass determination in the near future. These 
stars are the following : 


HD 90572 (ADS 7779) HD 97561 (ADS 8094) 
PA Sep n PA Sep 
58 082 
63 1°21 
66 1-06 
71 Ils 
75 1°07 
80 o’78 
82 0°72 
9! 
100 


288 1°05 
287 0°87 
284 0°68 
277 0°54 
276 0°42 
275 0°40 
266 0°39 
262 0°35 
254 0°34 
249 O25 
241 

203 

189 


ww AOU AN 


Royal Greenwich Observatory, 
Herstmonceux Castle, 
Hailsham, Sussex: 
1959 September. 
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STELLAR GROUPS, VI. 
SPACE MOTIONS OF THE DWARF A-TYPE 
AND GIANT M-TYPE STARS IN THE SOLAR NEIGHBOURHOOD 


Olin J. Eggen 
(Communicated by the Astronomer Royal) 
(Received 1959 October 9) 


Summary 


The A-type dwarfs and M-type giants fall in restricted regions of the 
(My, B—V)-diagram permitting the derivation of homogeneously accurate 
parallaxes and, therefore, space motions for objects in these two classes that 
are brighter than visual magnitude 5 and for which accurate radial velocities 
are available. ‘The space motion vectors, (U, V), of the A-type stars, which 
are a pure sample of the older population 1, are not randomly distributed and 
indicate the possibility that the majority of this population may exist in 
isoperiodic groups similar to the Hyades and Sirius groups. The restricted 
region of the (U, V’)-plane defined by the A-type dwarfs is used to isolate the 
disk population among the M-type giants. The M-type giants are almost 
equally divided between disk population and older population 1 stars. 





1. A-type dwarfs.—In a discussion of the luminosities of the stars brighter 
than visual magnitude 5 (Eggen 1957) the following two features of the 
distribution of the A-type stars in the (M,, B— V)-diagram were noted: 

(a) The so-called peculiar (Ap) and metallic-line (Am) stars in galactic 
clusters or with well-determined trigonometric parallaxes lie in a well-defined 
region and their luminosities show an average deviation of about o™-1 from 
the expression (Eggen 1959c) 

M, =7:0(B—V)+0™-4. (1) 

(b) The stars placed in luminosity class V of the MK-system occupy a 
region of the (M,, B—V)-diagram bounded by the Ap and Am stars (eqn. (1)) 
and the Pleiades main sequence. 

The Pleiades main sequence and equation (1) are represented in Fig. 1. 
Since the region in the (M,,B—V)-diagram occupied by the Ap, Am and 
AV stars is quite restricted, relatively accurate luminosities, and therefore 
distances, of these stars can be obtained from the spectroscopic classification 
alone. This is important since these stars are defining members of what has 
been called ‘‘ older population 1’’ (Oort: see O’Connell 1958) and, at least for 
the brighter examples, the distance is all that is lacking of the data needed for 
computing space motions. 

All the stars brighter than visual (HR) magnitude 5™-o and which satisfy 
the following conditions are listed in Table I: 

(i) Classified (MK) as Ap, Am or Bg-FoV. The authorities for the Ap 
and Am classifications are listed by Bertaud (1959) and those for the dwarfs are 
mainly by Slettebak (1954, 1955), Johnson and Morgan (1953) and de Vaucouleurs 
(1957); one or two are by Woods (1955). In a few cases stars were omitted 
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B-V ° ° +02 +03 


A. 








Fic. 1.—Schematic representation of the Pleiades main seq e and equation (1). The 
luminosities of the A-type dwarfs are bracketed by these features. 


because of disagreement among the classifiers as to the luminosity class and 
there are undoubtedly one or two mis-classified stars, which therefore fall outside 
the region defined in Fig. 1, remaining in Table I. 

(ii) Colours between B—V= —o0™-o5 and +0™-35. The main sources of 
the photometric data are Eggen (1955), Johnson and Morgan (1953), Hogg 
(1958), Stoy (1953), Johnson and Knuckles (1957), Naur (1955) and Osawa 
(1959). A few colours transformed from the values of C,/T listed by Hertzsprung 
(1940), are identified by colons. 

(iii) Proper motion listed in at least one of the fundamental catalogues given 
in the column labelled ‘‘S”’ as 1= GC, 2= N30 and 3 = FK3. 

(iv) Radial velocity of quality ‘‘a’’ or ‘“‘b’’ in the general catalogue (Wilson 
1953). A few recent determinations have also been incorporated into the means. 


TaBLe I 
A-type dwarfs brighter than visual magnitude 5™-o for which accurate 
radial velocities are available 
Name Vz B-V U-B_ Sp p Q He Me 
km/sec. 


14b —o-0s5 
8-ob —0-067 
8-8b +0102 
+O°1s5t 
+0'228 


6 And 4: +0°06 A2V 

o And +0°05 A2V 
«x Phe +018 yak A3V 

pw And 3°88: +012: AsV 
6 Cas . +0°17 A7V 
5 Cas . +0°14 AsV +0°297 

B Ari . +0°13 AsV +0°100 
50 Cas "00: +0°05: Be A1rV y —0°044 
58 And > +0°12: ive AsV . +0°157 
t Cas ° +0°13 Ap , —0o'013 

o Cet : +o°10 A2V , —o°14! 
mw Per 58: 0°00: AoV : +0°031 
+ 0°05 A1V . — 0°043 

38 Tau 3° +0°05 A1V : +0°003 
w Tau 4° +0°27 ies Am ; —o'o4! 
y Dor ; +0°34 FoV “ +o°101 
64 Tau : +016 A7V +o°114 
x Tau 4° +013 A7V , +0°099 

+0°17 A7V ° +0°10§ 

p Tau ‘ +0°24 . FoV ° +o°1o!1 

88 Tau +018 =e Am . +0°056 


l+++41+ 
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TaBLe I—continued 


A-type dwarfs brighter than visual magnitude 5™-o for which accurate 
radial velocities are available 


Veg B-V U-B_ Sp pO Me 


4°68 AsV + 36-9b 
Cam 4°47 AlV — 79a 
Aur 4°95: : AoV + 6-9a 
Tau 4°64 A7V +42°'2a 
491 oo BoV +24°2b 
4°70 ree j +29°6b 
Lep 3°72 . f — 1°6b 
Lyn 4°38: ‘00: 5 4 — 3°6a 
CMa —1-46 , — 76a 
Pic 3°28 +20°6b 
Gem 4°60: : EX: +25-9b 
4°46 ‘ — o8a 
Gem 3°57 g°2b 
Gem 4°16 5°7b 
Gem 18: d 30a 
Gem 4°95 8-ob 
3°89 . +10°0b 
Cnc 4°73 . AlV + 28-7b 
Vel 1-99 * AoV + 2-2b 
Hya 4°38: : +32°8b 
UMa 3:12 +12°2a 
Cnc 4:22 d — 13°8b 
UMa 4°46 — ola 
UMa 4°06 — goa 
UMa 4°82 + O° —15*1b 
Lyn 3°82 + 16b 
UMa 4°60: : ig: +23:1b 
4°51 ’ +24°4b 
LMi 4°48 / —17°8b 
LMi 4°74 +13°4a 
UMa 4'80: : AlV —17°4a 
Leo 4°42 f —10°2b 
Leo 2°55 : — 20°6b 
Leo 3°41 78a 
Leo 4°00: 53b 
Leo 2:18 4 orb 
UMa 2°44 r 4 12°9a 
UMa 3°30 j 12-9b 
Vir 3°92 aah ] 2°3b 
Vir 2°75 . F 19°7a 
Vir 4°95 j 1°6b 
Vir 4°38 ee 4 2°9b 
Cen 2°68 y 
UMa 4:01 AsV 75a 
Vir 4°04 : Ap 8-0a 
Vir 3°37 
Vir = 4°55 
Cir 3°20 
Vir 3°74 
Boo 4:30 
Ser 3°55 


~J 


ete Seawed- cove or B&B Dw 


= 





5892 
5980 
6031 
6117 
6129 
6153 
6168 
6254 
6324 
6436 
6446 
6554/5 
6581 
6771 
6789 
7OOI 
7056 
7069 
714! 
7235 
7420 
759° 
7653 
7724 
7871 
795° 
7990 
8075 
8140 
8151 
8322 
8576 
8585 
8709 
8717 
8728 
8830 
8949 
8984 


Name 


e Ser 
8 Nor 
% Sco 
w Her 
3 Oph 
w Oph 
o Her 
52 Her 
e Her 
69 Her 
53 Ser 
v Dra 
o Ser 
72 Oph 
6 UMi 
a Lyr 
¢ LyrA 
111 Her 
6 SerA 
{ Aql 
« Cyg 
e Pav 
5 Vul 
p Aql 
f Del 
e Aqr 
pw Aaqr 
6 Cap 
6 Ind 
6 Mic 
5 Cap 
PsA 
Lac 
Aqr 
Peg 
PsA 
And 
t Phe 
A Psc 


r 


Ve 


4°52 


B-V 


+017 
+0°24 
+o°10 
+o°o! 
+o18 
+0°13 
+0°02: 
+0°09 
— 0°03 
+0°04 
+ 0°05 
+0°265 
+008 
+015 
+0°o! 
0°00 
+018 
+O°11: 
+O°1§ 
+o°o! 
+0o°16 
— 0°02 
+016 
+0°06: 
+015: 
+ 0°03 
+0°32 
— 0°03 
+0°17 
+0°02 
+0°29 
+0°02 
+ 0°03: 
+0°05 
+o°o! 
+0°09 
+018 
+ 0°09 
+o°19 


TaBLe I—continued 


A-type dwarfs brighter than visual magnitude 5™-o for which accurate 
radial velocities are available 


U-B 


— 0°03 


+ 0°04 
~0-09 
—o'ol 


+008 


—o-'o1 
+0°17 


+0°09 


—o'o2 


+oO°o! 


+0°06 
0°00 
+0°08 


+0°07 


Sp 


Am 
Am 
Am 
Ap 

Am 


Ap 
BoV 
Ap 
AoV 
A2V 
Al1V 
Am 
A2V 
A4Vv 
A1V 
AoV 


AoV 
AsV 
Ap 
Am 
AoV 
A2V 
A3V 
AlV 
A3V 
FoV 
Ap 
A7V 


p Q 
km/sec. 


— 94a 
—15°5b 
— 5-7b 
— 66a 
—30°6b 
+ 2°§a 
—10°9b 

o’oa 
—25‘1b 
— gob 
+ 4°8b 
—15*6a 
—30°0b 
—23°9a 
— 7-6b 
~—13°9a 
— 26°0a 
—44°6b 
—44°0b 
— 26°3b 
—19°5b 
+ ob 
—20°9b 
—23°0b 
—25-0a 
— 16-0b 
— gb 
—10'9b 
~14°5b 
+ 2°3b 
— 63a 
+ 6°3a 
— 4:0b 
+18-0b 
—i0°1b 
+ 6°5a 
+12°5b 
+19°4b 
+12°4b 


Fe 


+0°125 
+0°002 
—o’o10 
+0°042 
—0°079 
+0°013 
- o°014 
rT 0°023 
—0°050 
—o°o4! 
+0°041 
+0°143 
—0°074 
—_ 0°064 
+0o°o10 
+0°199 
+0°023 
+0°069 
+0°042 
—o°o10 
+0°018 
+0°083 
+0°054 
+0°052 
+0°041 
+0°030 
+0°041 
+0082 
+0°105 
+0°071 
+0°261 
+0°064 
+0°134 
—o°o4I 
+0°073 
+0°232 
+0°089 
+0°039 
— 0°133 


He 


+0061 
+0°028 
—0°020 
—o'061 
+0°004 
+0°030 
+0°042 
— 0062 
+0°025 
+0°058 
+0°oo! 
+0°053 
—o'058 
+0°081 
+0°054 
+0282 
+0°021 
+O°117 
+0°032 
— 0°097 
+0°130 
— 0°133 
+0°006 
+0°056 
+Oor°or! 
—0°032 
—0°028 
—0°058 
—0°065 
+0°'001 
— 0°293 
—o°ol2 
+0°020 
—0°022 
+0oO°o17 
—o'162 
+0°098 
+0°004 
—0'148 


12 


12 


123 
123 
12 
123 
123 
123 
123 
123 
I 
123 
I 
13 
123 


Parallaxes of the 23 Ap and Am stars were obtained by the use of equation (1). 
The luminosity of each of the 88 dwarfs was assumed to be half-way between 
that given by equation (1) and that of a star on the Pleiades main sequence with 
the appropriate colour. The resulting parallaxes are listed in Table II together 
with the computed space motion vectors. ‘These vectors are directed as follows: 
U (away from the galactic centre), V (in the direction of galactic rotation) and 
W (toward the north galactic pole). Ninety per cent of the 111 stars in Table I 
are within 50 parsecs of the Sun. The (U, V)- and (V, W)-vectors in Table II 
are plotted in Figs. 2 and 3, respectively. 
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Tase Il 
Space motion vectors for the stars in Table I 
HR J V Ww HR 7 
km/ 
63 - 2 4356 23 
68 ng 4534 76 
100 —10 4554 5° 
269 + 4660 44 
343 46389 
403 4825/6 
$53 4828 
580 4963 
620 5028 
7°7 5062 
8o4 5105 
879 §107 
1046 5359 
125! 5463 
1329 $511 
1338 5733 
1380 5881 
1387 5892 
1427 5980 
1444 6031 
1458 6117 
1479 6129 
1568 6153 
1§92 6168 
1620 6254 
Ip2g 6324 
1762 6436 
2685 6446 
2238 6554/ 
2491 6581 
2550 6771 
2564 i 6789 
2740 5 - : 7001 
2763 - . 7056 
2852 
2890 
3067 
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The distribution of the (U,V)-vectors in Fig. 2 confirms two previous 
conclusions obtained from a preliminary discussion of the known Ap and Am 
stars brighter than V,=6™: 

(i) The (U, V)-vectors lie in a very restricted region of the (U, V)-plane— 
the limiting ellipse shown in the figure is that found previously for the Ap and 
Am stars (Eggen 1959b). 

(ii) The vectors are not randomly distributed within this ellipse. 














Fic. 














i i 
Vo -20 o +10 


Fic. 3.—The distribution of the (V, W)-vectors of the dwarf A-type stars. 


The (U, V)-limits in Fig. 2 are indicated by an ellipse because of the analogy 
with the velocity (dispersion) ellipsoid usually derived from statistical discussions 
of stellar motions. In fact, however, the distribution of the (U, V)-vectors of 
the dwarf A-type stars appears to be only roughly elliptical and the scarcity of 
stars with negative values of both U and V—that is, stars inbound toward the 
galactic centre on orbits with semi-major axes less than the Sun’s distance from the 
galactic centre—seems to be real. Moreover, the tilt, or deviation of the vertex, 
of the velocity (dispersion) ellipse from parallelism with the U-axis, usually 
found for A-type stars (cf. Alexander 1958), is clearly reflected in the tilt of the 
limiting ellipse and seems to arise partly from the concentration of vectors near 
(U, V)=(— 16,0) which represent the Sirius group (Eggen 1958, Paper I). 

It is also apparent from Fig. 2 why the Sirius and the Hyades, 
(U, V)=(+ 40, —17), groups could be detected by the ‘‘convergent-point’’ 
type of analysis described in Paper I since they lie at the extremes of the 
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(U, V)-distribution of the older population 1 stars. ‘The objects that represent 


these two groups in Fig. 2 are the following (Hyades cluster members are 
identified by an asterisk) : 


Hyades Sirius 
V WwW Y HR J V W 
—17 68 = — 3 And 
—13 r 804 -11 Cet 
—17 fer 2491 CMa 
14 3572 Cnc 
~13 ey 3662 UMa 
—13 : 3974 LMi 
~18 4357 Leo 
—19 4554 UMa 
—19 ‘ 4660 UMa 
5062 UMa 
§105 Vir 
6524 Her 


~NNOONN DONROW WN 


The non-randomness of the distribution of the rest of the (U, V)-vectors in 
Fig. 2 is undoubtedly real and is confirmed by the space motions of the early-type 
visual binaries (Eggen 1959¢). Since lines of constant Vx —7 and — 19 km/sec 
represent, over the small range of U in the figure, lines of constant semi-major 
axis (a) for galactic orbits with azo-93R© and o-g8RGO, where RO is the 
Sun-galactic centre distance, the concentrations of (U,V)-vectors in Fig. 2 
would seem to indicate that these values of V represent two preferred values 


of a for the older population 1 other than the Hyades and Sirius groups. Older 


® | 
@- 
© 





Fic. 4.—~Suggested interpretation of the distribution of the (U, V)-vectors of dwarf A-type stars 
shown in Fig. 3. The circles marked “‘ S"’ and ‘‘ H” contain the members of the Sirius and 
Hyades groups, respectively. 


population 1 therefore consists mainly of stars outbound from the galactic centre 
on orbits with these values of a and with eccentricities varying uniformly between 
about o and 0-05. However, by analogy with the known Hyades and Sirius 
groups, the interpretation of Fig. 2 that is shown in Fig. 4 seems more likely. 
Here seventy-five per cent of the stars in Table I are assigned to seven groups, 
other than the Hyades and Sirius groups. Some of the groups have so nearly 
a common orbital period and eccentricity that the slight blurring of the 
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(U, V)-vectors caused by the ~o™-5 uncertainty in the assumed luminosities of 
the stars taken from Fig. 1, and the uncertainties in the motions of even these 
well-observed bright objects, is enough to merge one group with another. It is 
obvious from Fig. 2 that the convergent-point method of assigning stars to groups 
is impracticable for these low velocity groups. In fact, it was previously found 
(Eggen 1959 a) that the Hyades and Sirius groups are at the limit of detectability 
by this method. The matter is also further complicated by the fact that the 
motion in W, although generally small, may be independent of that in U and V 
for these low velocity stars and the ‘‘convergent-point’’ of the group motion 
is consequently difficult to assign with accuracy. 

The possible groups of older population 1 stars, other than the Hyades and 
Sirius groups, which are indicated in Fig. 4 will be more fully discussed in a 
later paper of this series concerning the motions of the G-K giants. However, 
whatever interpretation is placed on Fig. 2, the obvious non-randomness in the 
(U, V)-distribution makes the classical, statistical treatment of the motions of 
these stars, for such constants as those of the solar motion, of dubious significance 
and it will not, therefore, be applied to the material in Table I. : 

The obvious absence of a ‘‘ solar group’’ with (U, V) = (0,0) in Fig. 2 should 
be noted. Since the stars in Table I are probably all younger than about ~ 
5 x 10° years and the Sun is nearly 5 x 10° years old, members of a solar group 
would not be expected in this material. 

2. M-type giants.—On the basis of current evolutionary theory we can expect 
the M-type stars that are classified as giants on the Mount Wilson system or as 
of luminosity class III on the MK system to be a mixture of three populations ; 
the older population 1, typified by the A-type dwarfs discussed above, the disk 
population, typified by the subgiants (Paper V, Eggen 1960) and a relatively 
small number of representatives of the halo population. Furthermore these stars 
probably fall in a rather restricted region of the (M,, B—V)-diagram. The older 
population 1, with ages up to about 5x 10* years, will mainly lie between 
M,=o™ and —1™ (cf. Hyades and Sirius groups, Paper 1) and the disk 
population, with ages up to about 10” years, will cover the same luminosity range 
(cf. 61 Cygni and y Leonis groups, Paper V). Because of the very large 
dwarf/giant ratio in the halo population globular clusters, and the relatively 
few halo dwarfs in the solar neighbourhood (Eggen and Sandage 1959, Paper IV), 
very few halo population giants, which have My, = —2™, will be found among 
the bright stars. Therefore, if we confine ourselves to the stars for which 
accurate proper motions and radial velocity material are available and that are 
redder than about B—V=1™-5 (that is, spectral type later than ~Ks), the 
true space motions can be reliably bracketed between those computed for 
assumed luminosities of M,=o™ and —1™, 

The stars brighter than visual (HR) magnitude 5™-o and with B—V >1™-5 
are listed in Table III. The arrangement and source ot the material in Table ITI 
are similar to that in Table I. Known supergiants of Class I have been omitted 
although a few stars recognized as of luminosity Class II are included and there 
are undoubtedly several others of that class not recognized. The spectral types 
are mainly by the Mount Wilson observers (Wilson 1953) and by Hoffleit (1942); 
very few M-giants have been classified on the MK system. Most of the stars 
whose magnitudes are followed by a ‘‘v’’ in the table were discovered to be 
small-amplitude variables by the Cape observers (Stoy 1953). 





Olin F. Eggen 
Taste III 
M-type giants brighter than visual magnitude 5™-o 


Vy  B-V  & 


p Q Hg 
km/sec 


m m ” 
4°80 +158 Moiil —45'8b +0°091 
+166 gM —22'5a —0'028 
4°85 +166 gM4 +12:2b —o’022 
4°27 +155 KslIll +32°3a §=6+ 0083 
4°78 +158 gMo +15°8a +0°006 
2°06 +1°58 Molll + oja +0178 
+160 gMs5 + I ¢ —0°'090 
4°05 +158 gM +180a +0°130 
5°00 +1°50 gMo —30°6a —0'032 
+1°65 M2 Ill —25°'9a —0'009 
+ 1°69 M4 Ill +28:2b +0131 
+1°58 gM3 +41'7a +0°053 
+1°65 gM2 +45°7a +0°047 
+164 Moiil +15°5a +0°051 
+ 1°60 Mo III +61-7a =+0°063 
+164 gM — 14b +0009 
+ 1°63 2M6 -— 7ob —o-002 
+150 gKs +24:1a +0°058 
+ 1°52 Ks III +54:1a +0°067 
+160 gM4 —33°6a +0°021 
+r5s0 KslIIl + roa +0°025 
+1°52 gKs +36°0a §=+0°062 
+160 gMr +37°7a +0°036 
+ 1°65 M3 II + o'9a = =6+0°003 
+ 1°66 M3 III +19g;0a —0°065 
+159 g—M3 +3454 —0'024 
+ 1°67 M3 III +54°8a +0°059 
+ 1°69 gMr +22'la —0o-'002 
+16 aMs +§3°0a §=6 +0104 
+1°51 Ks Ill +88:1a —0°064 
+ 1°60 Mo Ill —20°6a —0°033 
+ 1°62 gKs +32°7a —0°'017 
+1°74 gM2 +17'la 0009 
+1°51 K4 Ill +26°7a +0060 
+1°52 gzMo +25°3a —Oo°or! 
+165 gMo +15"9a —0-008 
+ 1°60 Ks Ill +24'4a —0'017 
+1°56 gM3 + 46a -—o-021 
+162 gKs —44°7b +0°039 
+ 1°55 gzMo +37°6a —0'220 
+168 gMo +20°0b —o-or2 
+160 KslIlIl +26-7b —0-026 
+158 gKs —13'9a8 —0°033 
+167 Moiil +23'4a —0°032 
+155 K4lll +41'1Ia —o-o81 
+1°57, gMo + 86a -—0°023 
+163 gMo +12°9a —0'009 
+155 Moll —20°5b —o-081 
+1°51 Ks III +346a —o0-'128 
+169 gMo +11-'2a —o-or! 
+166 gMo + 99a —o°0I9 
+160 gK6 —22°4a —0'037 


>Re 


Ser ox 4 





Stellar groups: V1 

TasLe 11 ]—continued 
M-type giants brighter than visual magnitude 5™-o 
y FF & 


e QO & 
km/sec 


m « 
+1°73.  gKs + 93a —0-028 
+163 KslIlIl —13'5a +0015 
+166 gM2 +15°6a —o-021 
+168 Molll + 7'2a —o-040 
+151 Molll +50°7a  —0'020 
+156 gM2 +37'4b +0°023 
+150 gMo — 24b —o-071 
+156 g@Ms + 71a 0234 
+1°61 gMo — 71b —0-078 
+160 M@3lll +21°3a ©+0°028 
+160 gM3 +17°6a —0°025 
+157, M3lll —178a —0-469 
+166 gMo — 16a —0'038 
+1°53 gKs —10°3a8 —OoI! 
+1°71 gM2 —26°9a —0'008 
+166 gM3 +18:2a  —o-r102 
+170 03=—« M2 III —171a —0'024 
+148 MriIll +40°78 —0°'047 
+160 KglIlIl — 56a —o'095 
+156 gM2 —43°6a -—0'024 
+150 gM3 —10°78 —o'0OI 
+1°57, gM3 + 56a —o°'0I7 
+ 1°65 aM4 — 4°38 —0°'073 
+155 gKs —29°44 —0°'004 
+162 KglIIl — 38b —0-039 
+168 gM: —38°7a —o-050 
+165 Mg4lIil —12°08 —0'004 
+160 gMo —19°98 ~—0'046 
+16 2M6 + 3°4a +0°027 
+150 KslIII + 3:14 —0o'180 
+160 gM2 — 21a +0°020 
+161 gMo + goa +0°040 
+163 gMr — 60a —o-014 
+150 Kglll — 67a —0°047 
+163 ~MBlll +0°019 
+167. gMr —o'olo 
+152 KsIIlI —oroll 
+152 gM3 — o3b —oo1s 
+158 M3lIlI —0°130 
+1°52 gMo +0°013 
+148 gMr 0-000 
+160 MgslIlIlI +o°o19 
+151 Moll —0°130 
+152 KslIll +0060 
+164 gM4 ; +0°034 
+166 gK6 —0°009 
+165 gM3 — 0°005 
+155 gMr +0°043 
+165 KslIlIl — 0-006 
+161 gMr —0°026 
+165  MriIill +0°018 
+161 gM4 +0°041 
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Taste I11—continued 
M-type giants brighter than visual magnitude 5™-o 


Veg B-V Sp Mg 


m * 
+1°61 M3 Il + 1° +0°135 
+1:59 Mo —~O'O1§ 
+164 gMa2 , + 0°006 
+ 1°66 Ma Ill ’ +0187 
+ 1°56 M2lll + 0°007 
+158 gM2 . +0'036 
+160 gMa2 , +0034 
+158 gM3 — 0°040 
+164 gM2 +0°050 


Tastz IV 


Space motion vectors for the stars in Table III 


U V W 
km/sec 
- 48 + 28 
20 +19 
13 —10 
10 “a 
° ~16 
18 —10 
I ~ I 
41 - 6 
8 
9 
14 
15 
7 
17 
45 


3 
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18 
15 

& 
12 
45 
20 
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° 
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20 
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76 
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24 
16 
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16 
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34 
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32 
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Stellar groups: VI 


TasBLe [V—continued 
Space motion vectors for the stars in Table III 


U V W ™ U 
km/sec 
eee 
+ 43 
— 2 
— 21 
— 21 

13 
23 
44 
8 
13 
2 10 
52 33 
10 , ° 
11 4 
18 24 
11 10 
3 22 
10 14 
- 6 18 
68 31 
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19 

94 

12 

18 

20 
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HR 


6815 
6832 
6855 
6868 
7187 
7495 
7635 
7650 
7673 
7951 
7952 
7980 
8080 
8225 
8582 
8636 
8679 
8608 
8775 
8795 
8834 
8860 
9064 
go89 


It 
23 


1 
16 
13 
16 
13 
10 
12°5 
10°5 
15 
12 
13 
I! 
36 
15 
17°5 
33 
12°5 
14 
II's 
11's 
13 
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TaBLe 1V—continued 
Space motion vectors for the stars in Table III 


U V 

km/sec 

+ I 
42 

11 

34 

21 

103 

20 

11 

7 

25 

14 

I 

8 

8 

10 


+34 


—24 
+11 
—2! 
+ 8 — 24 


Ww 


N J 


+ I 


9 
3 
° 
2 
I 
7 
° 
4 
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The space motion vectors of the stars in Table III are given in Table IV 
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where the first set of values for (U,V, W) were computed on the basis of an 
assumed luminosity of o™ and the second set on the basis of —1™, The values 
of (U,V) and (V, W) based on the assumption that M, = —1™ are plotted in 
Figs. 5 and 6, respectively. If the true values of the luminosities lie between o™ 
and — 1™, then the (U, V) and (V, W) distributions seen in Figs. 5 and 6 represent 
the most ‘‘expanded’’ form of these distributions for the stars in Table III; 
the use of the assumption M, =o™ makes all the space motions smaller and the 
distributions in the (U, V)- and (V, W)-planes somewhat more compact. The 
ellipse drawn in Fig. 5 is the limiting ellipse of the (U, V)-distribution of the 
dwarf A-type stars in Fig. 2. The M-type giants that fall within this ellipse are 
indicated by open circles; the crossed circles indicate the additional stars that 
would fall within the ellipse if the assumption that M,=o™ was used. The six 
stars falling near (U, V) = (+ 40, — 17) and indicated by filled circles are probably 
members of the Hyades group. Such members would be expected to have 
My —1™. 

The remaining stars of Table III are represented by crosses in Fig. 5. The 
majority of these stars are undoubtedly members of the disk population. It is 
of interest to note that the 113 stars of Table III are evenly divided between 
stars falling within and those falling outside the limiting ellipse for the A-type 
dwarfs. Therefore, since we can identify the stars marked as crosses in Fig. 5 
with the disk population, it appears that at least 50 per cent of the M-type giants 


in ‘Table III originated as main sequence stars of spectral type F2—Go and with 
masses ~ 1'5(). 
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Fic. 5.—The distribution of the (U, V)-vectors of the M-type giants computed on the basis of 
the assumption that the absolute magnitudes are —1™. The six stars shown as filled circles are probably 
members of the Hyades group and the crosses indicat’ probable members of the disk population. 
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Fic. 6.—The distribution of the (V, W)-vectors of the M-type giants computed on the basis of 
the assumption that the absolute magnitudes are —1™, The symbols are as in Fig. 5. 


Although the assumption that M, = —1™ will only be correct for a fraction 
of the stars involved, the distribution of the (U, V)-vectors within the limiting 
ellipse in Fig. 5 bears some resemblance to that of the same vectors for the dwarf 
A-type stars in Fig. 2. This is especially striking in the marked, apparent scarcity 
of stars with negative values of both U and V. There are also, apparently, very 
few M-type giant members of the Sirius group. Although we would expect 





Olin 7. Eggen Vol. 120 


such members to have M, ~o™ (Paper I), inspection of the vectors in Table V 
of that paper, computed on the basis of this luminosity, indicates only one member : 
61 Leo, HR 4299. 


Finally, a comparison of the (U, V)-distribution in Figs. 3 and 6 indicates 
that both the M-type giant and A-type dwarf members of population 1 show 
a small dispersion about a mean value of Wx —8 km/sec whereas, as expected, 


the disk population M-type giants show a larger dispersion about approximately 
the same mean value. 


Royal Greenwich Observatory, 
Herstmonceux Castle: 
1959 October. 
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47 TUCANAE: RADIAL VELOCITIES AND SPECTRAL TYPES OF 
INDIVIDUAL STARS 


M. W. Feast and A. D. Thackeray 


(Received 1959 October 9) 


Summary 


A discussion is given of the radial velocities and spectral types of 32 stars 
in 47 Tuc determined from 64 plates. The spectral types of 23 red giants 
lie in the range of MK types G8-—M2 and luminosity classes II-III. There isa 
slight CN discrepancy, but the spectra are much more normal than the 
extreme population 11 globular cluster stars. The metal/hydrogen ratio is 
unlikely to be less than 1/4 normal. The spectroscopic evidence suggests 
that the top of the giant branch lies below and to the red of that of other 
globular clusters in the HR diagram. One red giant, not known to be a 
variable star shows a variable spectrum. A mean velocity of ~ 14°2 km/sec for 
the cluster is deduced. A dependence of radial velocity on spectral type is 
found. It is unlikely that this is due to the wavelength system employed. 
From the velocity dispersion of the red giants it is found that under almost 
any assumption the mass of 47 Tuc is unlikely to be greater than 6-0 x 10° 
solar masses and taking into account the dependence of velocity on type may 
well be less than 2°5 x 10° solar masses. ‘The mass/light ratio (o°3 solar units) 
is much less than in intermediate spirals and elliptical galaxies. "The member- 
ship and spectrum of a bright B8 III star is discussed. ‘Two irregular and one 
semi-regular red variables have spectral types in the range M2-3 II-III; 
three regular (2004) red variables have types in the range M2-—3 II-Ib, and 
are normal Me variables. The velocities of these stars are discussed. The 
distance modulus of 47 Tuc is deduced from the RR Lyrae member 
(HV 814). It is shown that a distance of 4:4 kpc, which follows if HV 814 
has a median absolute visual magnitude of +0°5, fits best with the spectro- 
scopic absolute magnitudes of the red giants, the red variables and the B-type 
member. 





Introduction and observations.—Between 1951 and 1959 a considerable number 
of spectra of individual stars in the globular cluster 47 Tuc (NGC 104) were 
obtained. They were taken with the two-prism spectrograph at the Cassegrain focus 
of the 74-inch Radcliffe reflector. Altogether we now have for discussion 64 
acceptable spectra of 32 stars (three spectra of two RR Lyrae variables are not 
included in this number). The details are listed in Table I. Exposure times 
were usually between 3 and 5 hours, the photographic magnitudes being mostly 
near 13™-o. 

In the present paper we discuss the radial velocities and spectral types derived 
from these spectra. The programme was initially planned in conjunction with a 
photometric programme being carried out by Dr A. J. Wesselink which will be 
published separately. A preliminary condensed account of the spectral types and 
the photometry has already been given (1) and another paper (2) has discussed the 
membership of two RR Lyrae variables in the cluster and their spectra. This 
investigation has been restricted to the brightest stars in the cluster. These 
include a selection of the red giants, the long period and irregular variables and a 
bright blue cluster member. 





M. W. Feast and A. D. Thackeray 
Tasie | 


Spectra of individual stars in 47 Tuc 
Dispersion No. of Usual Projected Width of 
at spectra slit-width slit-width spectrum 
Hy (A/mm) (microns) (km/s) (mm) 


b 29 I 60 or1o 
e 49 39 go 0°20: 
d 86 24 175-200 0°23 


The stars studied are identified by number or letter in Plates 10 and 11(i)a 
(which are reproduced from a 74 inch Newtonian plate of the cluster in blue light). 
Three ‘stars’ are noted in Table II as blends of more than one star. All others 
observed were well resolved on the slit. 

The adopted velocities and spectral types are listed in Table II. 


Tasie II 
Radial velocities and spectral types in 47 Tuc 


Member Stars 


Velocity 
Star Sp. Date (U.T.) km/sec 


No.1 BS III: 1952 VII 29.077 —14°7 
1954 VI 29.118 —18'5 
VII 26.103 — 9°5 
—142 
K3 III (—?) VIII 17.028 —13°4 
IX 23°93! —18+1 
—15'8 
G8 III: (11?) VIII 1.077 —10°2 
G8 III IX 22.996 + 2°4 
XI 16.810 —17°3 
—- 74 
Ko III < 22.800 —13°2 
Ko II or III < 15.947 —17°2 
25.091 + 0O4 
— 84 
K3 II 6.046 —14'0 
27.074 —22°3 
29.926 —24°2 
29.983 —22°4 
—20°7 
K4 II-III 18.889 —25'8 
26.922 —25°2 
28.085 —14'8 
4.036 — 20°4 
—21'6 

Mz II (or 
II-III) 9.894 —12°8 
31.045 gst 
—- o4 





No. 5, 1960 47 Tuc: radial velocities and spectral types 


TasLe II—continued 
Member Stars—continued Velocity 
Sp. Date (U.T.) km/sec 


.13 G8: (III:) 1954 VIII17.036 — 45 
Yo. 14 Mr Il 1952 X 7.927 —23'3 
Jo.15 G8III 1952 VIII 12.042 —- 60 

G8 II-III 1954 IX 21.000 —10°8 

X 11.936 + 68 

— 20 

M IX 16.996 - 58 

K 23-774 —3rE 

Ko (III:) 25.879 —15°4 

31.090 73 

— 164 

Mr-2 III 2.911 —14°0 

M1 II 27.933 —29°0 
M2 II (or 

II-III) 19.999 —I13'1 

9.068 —29°6 

—214 

Yo. 36 M2 (no velocity) 

«gt. Kall 1957 19.093 — 18:1 
Yo. 52 Kolll 1957 VI 22.127 — 18 
vo. 58 Kolll 1957 VII 13.095 —1Iro 
Yo.79 Kail 1957 VII _ 16.095 —20°4 

Red Variables 


Velocity km/sec 
Star Sp. Date (U.T.) Emission Absorption Camera Obs, 
Vi M2-3 Ib 
(II-Ib?) 1952 VII 2.118 —19°0 
VII 14.090 —17'8 
1953 IX 3.119 —17°3 
IX 4.047 —21°§ 
IX 12.944 — 184 
xX 24.840 —12°7 
X 30.859 —15°7 
XI 20.853 
1954 XI 197290 —1I7'1 
1955 VII 12.118 —14°'0 


oR BeBe Re BeBe Be Be Ber 


enaanacananna 


—17'0 
M2-3 II-Ib 1952 VII 7.124 — 50 

VII 15.115 —12°2 

VIII 3.106 — 14°5 

1957 VIII 2.075 — 48 

— gt 

V3 M2-3 II-Ib 1952 XI 9.809 —40°1 
XII «5.799 — 46-0 

1953 XI 28.840 —44°1 

XII —s.1.797 —41'°8 

1959 VIII 15.052 — 34°9 


—41'°4 


hj try on 


he Rene 
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Taste I]—continued 
Red Variables—continued 
Velocity km/sec 
Star Sp. Date (U.T.) Emission Absorption Camera Obs. 
Vs M2 II or Ill 


CHD 1953. XII 13.785 ~14°6 d = 
1955 VI 25.125 —15°3 d T 
~—I5'0 
V7 M2: 1953 IX - 77 
v8 M2-3 II: 1953 X —378 —231 
Non-Members 
Velocity 
Star Sp. Date (U.T.) km/sec 
Comp b G8 V: 1953 IX 28.061 + 13 
XI —s.23.919 +27°3 
+143 
No. 22 F8V 1953 X 23.915 —13°5 


Notes on Table II1.—The table gives the star number, the spectral type, the time of mid- 
exposure (U.T.), the velocity (km/sec), the camera used (cf. Table I) and the initial of the 


observer. (B) in the last column indicates that the ‘ star’ observed is a blend of more than 
one star. 


i 
These results will be discussed in the following sections; it is convenient, how- 
ever, at this stage to draw attention to stars No. 22 and 6 (a Harvard comparison 


star). These two stars are classified as late type dwarfs and their types leave no 
doubt that they are foreground stars. The magnitude of No. 22 is about the same 
as that of the red giant members and that of 5 is rather brighter than this. The 
velocity of 5 is definitely discordant from the adopted cluster velocity of 
— 14:2 km/sec (see below). Though the velocity of No. 22 agrees rather well with 
the cluster velocity this agreement must be fortuitous. These twostars are omitted 
from further discussion. 

Spectral types of the red giants.—Spectral types of the red giants were obtained 
onthe MK system by comparison with standard stars. The 23 stars inthe category 
are a selection of the brightest non-variable stars, distributed fairly uniformly about 
the centre of the cluster and fairly close to the centre compared with the distance to 
which the cluster can be traced (26'-8). 

It should be emphasized that the photograph reproduced in Plates 10 and 
11(i)a is a light exposure to show the stars studied and does not give a correct im- 
pression of the size of the cluster. The stars range from the brightest non-variable 
stars to stars about a magnitude fainter. ‘The spectral types in Table II show a 
continuous range of types from G8 to M2 with luminosity class between II and 
Ill. There are 5 stars of type G, 11 of type K and 7 of type M. It is well known 
(cf. (3)) that some globular clusters show marked spectral peculiarities. Weak- 
ening of the metal lines with respect to hydrogen has been particularly noted (e.g. 
in Mgz2) and this is usually interpreted as a real metal deficiency in these stars. The 
giants in 47 Tuc were therefore examined for any sign of this effect. In particular, 
following Deutsch (3) we have studied line ratios involving hydrogen and the 
metals (e.g. CH or Fer 4325/Hy) and ratios involving metals alone (e.g. Cri 
4254/Fe1 4250). Whilst it is found that in a cluster such as Mg2 much earlier 
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74-inch Newtonian photograph of 47 Tuc (blue light) identifying stars discussed in text. Some 
stars not discussed are given Radcliffe numbers. The Bailey variables are marked“ V'"’ and HV 810 
(RR Lyrae type) is also marked. 


M. W. Feast and A D. Thackeray, 47 Tucanae. 
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(i) 

(a) Enlargement of central region of 47 Tuc (blue light); (b)-+f) Spectra in the violet region 
(approx. 4100-4350A) of stars in 47 Tuc and comparison stars. All spectra are enlargements from 
¢ camera plates; (b) 47 Tuc V.1 (2124 Me variable) with Iron Arc comparison; (c) RV Pup 1904 
galactic Me variable; (d) 47 Tuc No. 6 G8 III (with Iron Arc comparison); (e) 8 Boo standard 
G8 III+CN-1; (f) 47 Tuc No. 14 M1 II. 
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(ii) 

Spectra in the blue region (approx. 4380—4980A) of stars in 47 Tuc and comparison stars. All 
spectra are enlargements from c camera plates. 

(a) 47 Tuc V1 (212d Me variable) with Iron Arc comparison; (b) RV Pup 1904 galactic Me 
variable; (c) 47 Tuc V3 (1924 Me variable); (d) 47 Tuc No. 14 M1 II; (e) x Ser standard Mi III; 
(f) 47 Tuc No. 26 M2 II (or I-IIl); (g) 47 Tuc No. 3 K3 III (—); (hk) & Boo standard 

G8 IIT+CN-1; (i) 47 Tuc No. 6 G8 III with iron are comparison. 


M. W. Feast and A. D. Thackeray, 47 Tucanee. 
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types are indicated by the former ratigs than the latter, in the case of 47 Tuc the two 
methods agree well. Thus 47 Tuc is much closer to the ‘‘ normal ’’ population 1 
giants than to the extreme population 11 giants of some globular clusters. Exam- 
ination of the 47 Tuc G—K spectra in the light of the criteria established by 
Keenan and Keller (4) for high velocity stars shows that whilst the 47 Tuc spectra 
are much more like the normal giants of the MK system than the metal weak 
globular cluster giants, there is a slight weakening of CN in the 47 Tuc stars of a 
comparable amount to that found in the high velocity stars. In the M type stars 
the types were determined in the usual way from the strength of the TiO bands. 
These types were found to agree well with those determined from atomic line 
ratios (line criteria for M types being taken from Keenan’s work (§)). In Plate 
11 (i) and (ii) some typical spectra of stars in 47 Tuc are reproduced. The spectra 
shown were all taken at an original dispersion of 49 A/mm at Hy. Some standard 
stars are also shown for comparison including 6 Boo which is a high velocity 
star with a slight CN weakening (G8 III + CN — 1) (see Keenan (6)). 

The fact that the giants in 47 Tuc appear from their spectra to be nearer to 
‘*normal’’ population I giants than to the extreme population 11 stars lends great 
interest to the HR diagram for this cluster. Whilst a full discussion must natur- 
ally await the publication of the photometric results certain facts can be ascertained 
from the spectra themselves : 

(1) 47 Tuc contains non-variable M giants, that is, the stars at the tip of the giant 
branch are later in type than those found in any other globular cluster so far ; 

(2) the luminosity classes assigned to the giants in 47 Tuc (II-III) are fainter 
than those generally assigned in other clusters (I]—Ib) (cf. (7)). 

We may deduce from these spectroscopic results that : 

(3) the giant branch in the HR diagram in 47 Tuc lies below and to the red of 
those of other globular clusters thus far studied. The preliminary photometric 
results (1) bear out these deductions. 

These observations appear to fit well with presently held ideas as to the inter- 
pretation of the HR diagram of globular clusters. ‘The spectral peculiarities of 
globular cluster giants is attributed to a low metal/hydrogen ratio, and Hoyle and 
Schwarzschild (8) showed that the position of the giant branch in the HR diagram 
is sensitive to this ratio. ‘The giant branch is depressed and extends farther to the 
red for the clusters with a more normal (population 1) ratio. The exact position 
of the giant branch computed theoretically is critically dependent on the mixing 
length in the stellar model (9) and has not yet been fixed completely satisfactorily. 
Kinman (7) has investigated in a semi-quantitative manner the spectral types of 
stars in globular clusters and their relation to the metal/hydrogen abundance 
ratio. He finds no peculiarities in the integrated spectrum of 47 Tuc and deduces 
that the metal/hydrogen ratio is unlikely to be less than 1/10 normal. Our 
observed CN weakening in the giants is near the limit of detection on our plates and 
places the 47 Tuc stars in the same category as the high velocity stars showing this 
abnormality. The work of Schwarzschild et al. (10) and the discussion by Keenan 
(11) indicate that in the giants showing extreme high velocity characteristics the 
metal/hydrogen ratio does not drop below that of ordinary giants by a factor of 
more than 3 or 4. Until a proper curve of growth analysis of 47 Tuc can be carried 
out, we may tentatively adopt a metal/hydrogen ratio of not less than about 1/4 
normal, since the CN weakening in the 47 Tuc stars is not particularly marked. 
It should be noted that the above remarks apply only to the G and early K stars in 
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47 Tuc. For stars of later type additional sources of opacity other than hydrogen 
are undoubtedly present and an analysis such as that given by Kinman breaks 
down. It may be recalled that the spectra of the RR Lyrae stars in 47 Tuc (2) show 
evidence of a moderate abundance anomaly. We may conclude that the cluster 
shows some metal deficiency compared with normal population 1 stars but to a 
less extreme degree than any other globular cluster so far studied. 

It should be mentioned that star No. 10 appears on at least one spectrogram 
to have exceedingly weak hydrogen lines whilst on other spectrograms these appear 
to be normal for the type. The most probable explanation is that there is weak H 
emission occuring from time totime. Perhaps the star is variable and it would be 
interesting to examine it photometrically. 

The distance modulus of 47 Tuc.—In later sections we shall require a value for 
the distance of 47’ Tuc. Shapley (12) deduced a modulus of 14-7. This distance 
modulus would make 47 Tuc extremely bright with an integrated photographic 
absolute magnitude of — 10-0 (13) compared with an average value of about 
— 8-2 (6). It would also make the Me variable members much brighter than simi- 
lar stars in the general field (see below). Shapley’s estimate was based on the 
magnitude of the brightest stars, assuming this to be the same in all clusters. The 
dependence of the position of the giant branch on chemical composition discussed 
above renders this assumption invalid. In particular the giants in 47 Tuc will be 
absolutely fainter than those in other clusters so that Shapley’s modulus will have 

|! been overestimated. Preliminary consideration (1) of the luminosity classes of 

, the 47 Tuc giants and of the colour magnitude diagram led to a modulus in the 
range 13™-5 to14™-o. Since that time it has been shown (2) that two RR Lyrae 
stars are cluster members. The c-type variable which has been studied photo- 
electrically has a median magnitude (V) of 13™-9. Provided we adopt a value for 
the absolute magnitude of the RR Lyrae stars (M_,) this star then gives a distance 
modulus. Unfortunately the absolute magnitude of the RR Lyrae stars at median 
light remains in some doubt. A value of o™-o was adopted for many years but 
recently values near + 0™-5 have found favour ((14), (15); see also (16) and (7)). 
In the following we shall derive quantities based on both assumptions and it 
will appear that in general the results are best interpreted if the absolute magnitude 
of the RR Lyrae stars in 47 Tuc is +0™-5, 

It is not known exactly how much correction for absorption is necessary. 
However studies of foreground stars in the direction of the Small Magellanic 
Cloud (i.e. in the same region of the sky) indicate (17) a probable absorption 


correction of about o™-2. The distances deduced from these data are given in 
Table III. 


Tasie III 


Distance of 47 Tuc and absolute magnitudes 

Assumed M, of RR Lyrae Stars oM-o 
Uncorrected distance modulus 13-9 
Corrected distance modulus 13™-4 
Distance 5°5 kpc 
Mean m, for red giants studied 11™-3 
Mean M, of red giants — 20-6 
Mean M, of red giants from spectra (II-III) —1M-4 
Integrated photographic apparent magnitude of cluster (13) 4™-7 
Integrated photographic absolute magnitude 

of cluster —g@-2 
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Since the preliminary photometric results (Wesselink, unpublished) indicate 
that the red giants studied spectroscopically have an average visual magnitude 
(m,) of near 11™-3, the absolute magnitudes for these stars on the two hypotheses 
are shown in Table III. An absolute magnitude of — 2™-6 would place the stars 
in luminosity class II according to the calibration of the MK system by Keenan 
and Morgan (18). But an absolute magnitude of — 2™-1 would be in good agree- 
ment with the assigned luminosity class (II I—I1) as discussed above. 

Arp (19) finds the integrated absolute magnitudes of 7 clusters ranging from 
—7°6 to —8-9. Since they depend on an absolute magnitude for the RR Lyrae 
stars of o™-o these figures must be compared with the value derived for 47 Tuc on 
our first assumption, viz. —g™-2 (assuming the absolute magnitude of the RR 
Lyrae stars does not change from cluster to cluster, as is possible). The inte- 
grated magnitude of 47 Tuc is, then, greater than that of the clusters in Arp’s list 
but not surprisingly so if account is taken of the considerable spread amongst 
these 7 clusters. 


The radial velocities of the red giants 


(a) ¢ camera velocities (49A/mm).—All the plates (19 plates of 15 stars) were 
measured on the Observatory’s projection measuring machine (20) by MWF and 
16 of these plates on a Hilger machine by ADT. ‘The measurements were reduced 
in both cases using wavelengths deduced froma Boo(21). Inmany cases weighted 
means of lines seen multiple at high dispersion had to be used. Using the 16 
spectra in common between the two measurers we deduce the following means 
from these spectra: 

ADT — 17°6 km/sec 
MWF — 13°9 km/sec. 


When, however, the measurements were reduced, using for each spectrum only 
lines measured by both measurers, we find for the means : 


ADT — 153 km/sec 
MWF — 15°5 km/sec. 


This improvement in the agreement of the two sets of measures is probably due to 
the rejection of spurious lines, lines at the margin of visibility, misidentified lines, 
etc. Ithas, therefore, appeared best for the c camera results to adopt the velocities 
derived in this second manner and to take a straight mean of the results of the two 
measurers. A list of the wavelengths used is given in Table IV. 

Where possible velocities were also derived from wavelengths for G and K stars 
at this dispersion used at Victoria (22). Unfortunately, sufficient Victoria lines 
to give a significant velocity could not be measured in a good many of the spectra, 
so that the results obtained in this way were not used in deriving the mean velocities 
of the stars. However the ‘ Victoria’’ velocities will be needed in the later 
discussion. 


(b) d camera velocities (86A/mm).—The 16 plates of 15 stars were measured on 
the projection micrometer. The wavelengths adopted for the reductions were 
those used at this dispersion by Kinman (23) in his globular cluster work. ‘These 
wavelengths are provided with corrections according to spectral type to Ko but 
not later. The Ko corrections have therefore been used for all the later type 
stars. 
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Wavelengths used for the reduction of the c camera plates (with principal 
contributors in the case of « Boo) 


4667°84 Fe+Ti 4430°62 Fe 
46°80 Fe+Cr 29°83 Call 
05°77 Ni+V 22°01 Fe+V 
4586-14 Ca+V+O 18-04 Titi+Fe 
71°32 Mg+Ti 1511 Fe 
64°67 Cr+ Fe+ Tit. 08-47 Fe+V 
49°35 Tiu+Ti 04°74 Ti+ Fe 
45°32 Tin+Cr+V 00°49 Sc 11+ SiF 
22°78 Ti 4395715 Tiu+V 
18-04 Ti 83°50 Fe 
12°75 Ti 7598 Fe 
or-29 Ti 67-86 Several 
449653 Ti+Cr 51°49 Cr 
94°56 Fe 40°48 Hy 
88-61 Several 37°50 Fe+ Ti 
68-52 Till 25°55 Fe+Scir 
62°05 Fe+V 20°85 Scu+Tim 
54°76 Ca 4314°78 Ti 
43°82 Ti 4254°36 Cr 
35°64 Several 50°46 Fe 


(c) Discussion.—The mean velocities of the cluster derived from c and d plates 
appear to be significantly different. We find a mean difference between the results 
of the two cameras of c—d= + 10°0 + 3°5 (s.e.) km/sec. A difference of about 


10 km/sec between c and d camera plates of the same star has also been found in the 
case of early type stars in the Magellanic Clouds (i.e. at the same zenith distance) 
(unpublished material) whereas the standard B star + Sco which passes near the 
zenith in Pretoria shows very good agreement between c and d camera velocities. 
We may infer from this that the effect depends on zenith distance. (‘The con- 
siderable amount of earlier work with the c camera shows that the effect is not in the 
velocities from that camera.) It may be that atmospheric dispersion affects the 
guiding at the wide slit used with the d camera. Atmospheric dispersion might be 
expected to give an effect of the right sign and magnitude. However Kinman (23) 
has argued from the agreement of the velocities derived from integrated spectra of 
clusters and from the spectra of individual members that any effect due to atmo- 
spheric dispersion is unimportant. If such is the case the effect must have some 
other cause. In view of these findings we have applied a correction of + 10 km/sec 
to all the d camera velocities. The final means (‘Table IT) are derived giving equal 
weight to the c and d camera velocities. 

The mean velocity of 47 Tuc derived from 22 red giants (35 spectra) giving 
unit weight to each star is —14:2km/sec. Kinman (23) gives a velocity of 
—24km/sec from 5 accordant integrated spectra of 47 Tuc obtained with the d 
camera. It will be noted that if our correction of + 10km/sec is applied to this 
mean it agrees exactly with ours. 

Velocity dispersion of the red giants.—One of the aims of the present investiga- 
tion was to attempt to obtain an estimate of the mass of 47 Tuc from the velocity 
dispersion in the cluster. Before proceeding to this attempt we investigate the 
velocities for a dependence on spectral type. 
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Plots of the velocities of the individual red giants against spectral type are shown 
in Fig. 1. This figure appears to show, and statistical tests indicate, a significant 
variation of the mean velocity with spectral type. As will be seen from Figs. 1 to 4, 
this variation is shown not only by the adopted velocities but also (i) by the 
c camera velocities alone derived as discussed above from the « Boo wavelengths 
(ii) by the d camera velocities alone derived from Kinman’s globular cluster 
wavelengths; and (iii) by the c camera velocities reduced according to the Victoria 
wavelengths. The fact that the different dispersions and different wavelength 
systems all show the same effect indicates that it is unlikely to be due to erroneous 
wavelengths. The c camera line residuals were examined to see if any residuals 
showed large dependence on spectral type but that was not found to be the case. 
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Fic. 1.—Plot of the adopted velocities of individual red giants in 47 Tuc against spectral type. 
The mean velocity of the stars is indicated by the horizontal line and a line joins the means for the 
groups G8. Ko. K2—K4 and M1-Mz2. 
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Fic. 2.—Same as Fig. 1 but using only the c camera velocities from Table II (« Boo wavelengths). 


In Mgz (24), the only other cluster for which radial velocities of an appreciable 
number of stars have been determined, the spectral types of these stars have not 
been given so a test for this effect cannot be made. In Fig. 5 the velocities are 
shown on the plane of the sky. There appears to be no significant variation of 
velocity with position in the cluster. In particular, the cluster does not exhibit 
significant rotation. 

Although an instrumental origin for the dependence of radial velocity on spec- 
tral type cannot be completely disproved, it is of interest to note that a physical 
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origin can be suggested which is in line with modern ideas concerning atmospheres 
of red giants. It is known that late type giants of comparable absolute magnitudes 
to those in 47 Tuc are suffering mass loss with velocities of the order of 10 km/sec 
((z5) etc). In the case of the globular cluster giants studied, these lie at the tip of 
the giant branch where helium burning accompanied by a considerable change in 
stellar structure is believed to occur (8). The details of stellar evolution at this 


point are only qualitatively defined but the possibility of considerable mass loss 
here has been discussed (cf. (26)). 
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Fic. 3.—Same as Fig. 1 but using only the d camera velocities from Table I. 
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4.— Same as Fig. 1 but using only those c camera velocities determined with the Victoria 
system of wavelengths. 


It appears therefore feasible that an outward flow of matter in the atmosphere 
of the late K and M giants in 47 ‘Tuc would explain the observations. _ If this is the 
case we might expect to find a relation between the velocity derived from any given 
line and the line’s excitation potential (Deutsch (25) lists examples of such a relation 
for Me variables and some other stars). Unfortunately the present data, in 
particular the fact that so many lines are blends, is insufficient to test for this 
relationship. 

If we accept the correlation of velocity and spectral type as real, whether instru- 
mental or physical in origin, attention must be paid to it ; for it will have a profound 
effect on the velocity dispersion and the mass deduced for the cluster. Bearing this 
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Fic. 5.—Plot of the stars with measured radial velocity on the plane of the sky. The radial 
velocities to the nearest km|sec are indicated. The red variables are marked V and the two field 
stars F. The line represents the major axis of the cluster (H. Shapley and H. B. Sawyer, H.B. 852 
1927). 


in mind we now proceed to a discussion of the velocity dispersion in 47 Tuc. We 
have 35 spectra of 22 red giants not known to be variable. For 2 stars there are 
4 spectra, for 7 stars 2 spectra, and for 13 stars only 1 spectrum each. From 
these data we calculate the dispersions listed in Table V. _ In this table a, is the 
observed velocity dispersion obtained by giving equal weighttoeachstar. » is the 
dispersion per spectrum of unit weight due to observational error (and possible 
real variation in the measured stellar radial velocity). This quantity is obtained 
from the g stars with more than one spectrum each. If o is the true velocity 
dispersion then we have (27) 


i 
B=a2— 2 yo9°/t (1) 
where i= the number of plates of unit weight per star 


v,;= the fraction of stars with i plates of unit weight o,, given in Table V is 
defined by 


oy = > 44,06°/i. 
The dispersions given in Table V have been calculated in the usual manner, that is, 
from the sums of the squares of the deviations from the means (suitably weighted) 
divided by the number of degrees of freedom in each case. The errors quoted are 
standard errors. 

We have derived the true dispersion (a) both ignoring and taking into account 
the variation of velocity with spectral type discussed above. In the case where we 
ignore the variation of velocity with type (i.e. where we regard it as a chance 
fluctuation in the velocities sampled) we give for comparison the results calculated 
giving equal weights to the c and d camera plates and those obtained giving doubled 
weight to the c camera plates. The results do not differ significantly. In taking 
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TaBLe V 


o,* oa Corresponding 

mass (cf. 
Table VI). 

No correction applied for 

dependence of velocity on 

spectral type. 

(1) Equal weight to c and 

d plates. 69°6+22 62:9+25 48°6+19 21'0+29 
« (2) Double weight to 


c plates. 65°2+20 88°9+35 50°7+20 14°§+28 
After correction for dependence 
of velocity on type. Equal 
weight to c and d plates. 38's +13 62:9+25 48°6+19 


Moz reduced by the same 


method 41°2+16 21-0o+10 15°4+8 


Units are (km/sec)* 


into account the observed variation of velocity with type we have divided the stars 
into four groups (G8, Ko, K2-4 and M1-2) and obtained a joint estimate of ¢, 
f It will be noticed that in this case since o, < c,, the value of 
6 (the true dispersion) derived is negative though of course not significantly less 

The considerable difference between the values of o? obtained by 


ignoring or taking into account the variation of velocity with spectral type will be 


om these four groups. 


than zero. 


noted. It is therefore quite clear that any mass determination must be treated 
with reserve. However, it will be shown that a good upper limit to the mass may 
still be determined. 

The only other work so far published on the velocity dispersion in globular 
clusters is the investigation of Mg2 by Wilson and Coffeen (24) whose results have 
also been re-discussed by Schwarzschild and Bernstein (28). Forcomparison 
with our results we have reduced the Mgz results according to the procedure 
adopted for 47 Tuc ind outlined above. From 23 spectra of 15 stars we deduce 
the dispersions given in Table V for Mg2. Our velocity dispersion and its error 
for Mgz2 are somewhat larger than that given by (24) and (28) (o?=19°5 + 
(km/sec)*) for the following reasons. In the case of the dispersion itself we have 
derived o, and o, from the squares of the deviations from the mean whereas 
Wilson and Coffeen (24) preferred to employ the modulus of these deviations. 
Further our result is obtained by combining ¢, and o, according to equation (1). 
The error quoted for o? by Schwarzschild and Bernstein (28) (viz. 9 (km/sec)*) is a 
probable error and therefore needs multiplying by 1-48 to be comparable with our 
quoted standard error. Apart from this trivial difference there is another more 
important difference. The variance, o”, is obtained as the difference of two 
variances and its standard error is thus determined from the standard errors of the 
two variances differenced. The manner in which these two standard errors are 
combined to form the standard error of o* depends on whether or not the two 
variances differenced are correlated or not. Schwarzschild and Bernstein (28) 
assume a fairly high degree of correlation in deducing their uncertainty for o*. 
If they had taken the variances differenced as uncorrelated then the uncertainty 
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of o® would have been increased by a factor of about 4/2. We have calculated the 
standard error of o? assuming that the two variances differenced in equation (1) 
are uncorrelated, as appears to be indicated by statistical theory (29). Should 
correlation exist, our standard error is then an over-estimate. The matter is of 
some importance since the value of the standard error of o* enters into the deter- 
mination of an upper limit of the mass of the clusters. Schwarzschild and 
Bernstein’s procedure would lead to a somewhat lower upper limit. 

In determining the mass we have followed Schwarszchild and Bernstein (28) 
and used the virial theorem in the form 


V?=GM/# (2) 
where G =the gravitational constant ; 
M =the mass of the cluster 


{Praot and F=a({" ay)" /[" 5 @) 


S(q) being the number of stars in a strip of unit width at perpendicular distance q¢ 
from the centre of the cluster. Thus 


S(q)=2 i) * n(d) d(I— q*)2 (4) 


q 
where n(J/) is the number of stars per unit area at distance / from the centre of the 
cluster. 

For Mg2 Schwarzschild and Bernstein obtained n(J) from published star 
counts and deduced # by numerical integration. For 47 Tuc we have used the 
surface photometry of Gascoigne and Burr (13) who tabulate the equivalent number 
of tenth magnitude stars per unit area at different distances from the cluster 
centre. We therefore assume a constant mass to light ratio throughout the 
cluster. Whilst this may not be strictly true, in view of the radial variation of 
colour in the cluster (13) the error thus introduced will probably be negligible 
for our purpose. 

Then by numerical integration of (4) and (3) we find 

7 = 409". 

‘The conversion of this to linear measure depends on the adopted distance modulus 
discussed above. For the two alternative distances we obtain the results listed in 
Table VI. This table also contains the mass deduced from the dispersions of 
Table V according to equation (2). Since the value of o* derived if allowance is 
made for the dependence of velocity on type is negative we in principle derive a 
negative mass. The result is not significantly different from zero and has not 
been entered in Table VI. The result still however allows a determination of the 
upper limit of the mass and this is included in Table VI. The upper limits to the 
mass are determined from the quoted standard errors. If the distribution is 
normal there is 0-84 probability that the mass is not greater than the deduced mass 
by more than its standard error and 0-98 probability that it is not greater by more 
than twice its standard error. 

Despite the problem of the variation of the velocity of the 47 Tuc stars with 
spectral type, Table VI shows that we can still deduce useful upper limits to the 
mass. ‘Thus it would appear very unlikely under any assumption that the mass is 
greater than 6-0 x 10° solar masses and on our most favoured distance there appears 
a good probability that it is less than about 4-5 x 10° solar masses. If the variation 





M. W. Feast and A. D. Thackeray 
Taste VI 


Limits to masses of 47 Tuc and Moz 


47 Tucanae 
Adopted M, of RR Lyrae stars oM-o +o™-5 


7 10°9 psc 8-6 psc 


Mass (uncorrected for variation of velocity 
with spectral type and equal weight to 


c and d plates). Solar masses (M,) (1°6+2°2) x 10° (1-3 41°7) x 10° 


Mass (uncorrected for variation of velocity 
with spectral type and double weight to 


c plates). Solar masses (M,) (1*1 +2°2) x 108 (o-9 +1°7) x 10° 


Upper limits to mass: 


With 0°84 probability M, <3°8 x 10° 


<3°0x 105 
With 0-98 probability M, < 6-0 x 10° 


<4°7x 10° 
Mass (M,) after correction for variation of 


velocity with spectral type and with equal 
weights to c and d plates. 


Upper limits to mass: 
With 0°84 probability M, 
With 0-98 probability M, 


Mass (M,,) to light ratio (solar units) 


Moz re-reduced 
Mass of Moz (solar masses) ‘8+1°2)x 105 


Upper limit to mass of Moz: 
With 0°84 probability <3°0x 10° 
With 0-98 probability <4'2x 10° 


Mass/light ratio (Mg2) 1'°0+0°7 


of velocity with type is considered established then we reduce the upper limit even 
more; the mass being unlikely to be more than about 2-5 x 10° solar masses. As 
shown in Table VI results of the same order of magnitude are found for Mg2 and 
the results may also be compared with the mass of M3 (1-4 x 10° solar masses) 
derived from the luminosity function (30). 

The mass to light ratio in 47 Tuc is low (see Table VI). This confirms the 
previous estimates (31) that this ratio is much less in globular clusters than in 
elliptical galaxies. A recent analysis by Minch (32) suggests that mass/luminosity 
ratio may be the same in intermediate spirals and giant ellipticals and about 
300r 40. Hence it would appear that the ratio is significantly different in globular 
clusters from either of these two types of systems. 

The radial velocity and spectral type of the bright blue star.—Star No. 1 which 
is easily the brightest star photographically has a visual magnitude of 10°26 
(Wesselink, unpublished) and is quite blue. As seen in Table I, its velocity is 
~—14'2km/sec. This is in good agreement with the mean velocity of the red 
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giants (— 14:2 km/sec), the exact agreement of course being fortuitous. There 
is therefore good reason to believe the star to be a cluster member especially as 
to the evidence of the radial velocity we may add the facts that the star lies quite 
close to the cluster centre and that blue stars of this magnitude are relatively rare 
at this galactic latitude. We have classified the star as B8 III (H, Het, Mgu, 
Catt and Sitt all being found) and no peculiarities were noted in the spectrum. 
Two blue stars such as this have been proved as members of M13 and M3 
respectively from radial velocity considerations and more stars of this sort are 
believed to be members of other clusters (cf. (33)). The star in M13 is classified 
as B2V and that in M3 as O8. In neither case have peculiarities been noted 
either in the spectra or the colours. The presence of such stars has been generally 
regarded as a problem as it seems they must be considerably younger than the 
maiority of the stars in the clusters. A similar situation is found to exist in young 
galactic clusters (34) where some stars are also found which are apparently younger 
than those of the main body of the cluster. Perhaps these bright blue stars in 
globular clusters represent a second generation of stars formed out of matter 
ejected by stars that have passed through the red giant phase. Further 
investigation of these stars is clearly desirable. 

From the measured visual magnitude of Star No. 1, we deduce its absolute 
magnitude on our two alternative distances (Table VII). 


Taste VII 
The absolute magnitude of the bright blue ® No. 1 (B8 III) 


Adopted M, of RR Lyrae stars oM-6 +oM-5 
Absolute visual magnitude of No. 1 — 3™-6 —3™1 
Absolute visual magnitude from spectrum —3™-0 


The spectral type gives the absolute magnitude shown according to the 
calibration of Morgan and Keenan (18). The result may be taken as another 
indication that the observations are better satisfied with the smaller distance. 

Radial velocities and spectral types of the red variables.—47 Tuc contains 
several bright red variable stars and a considerable amount of time has been spent 
studying them spectroscopically. The red variables of greatest interest are 
Bailey Nos. 1, 2, 3 (35) which appear to be regular long period variables with 
periods near 200%. ‘The irregular variables Bailey Nos. 5 and 7 and the semi- 
regular variable Bailey No. 8 have also been observed. Altogether we have 22 
spectra of the 6 stars, 16 of them on the c camera and the rest on the d camera 
with all but 4 of the spectra referring to variables 1, 2 and 3. The variables are 
identified on Plates 10 and 11 (i) and Table VIII is reproduced for convenience 
from data in Mrs Hogg’s catalogue (36). 

Consider first the three regular variables. Their nearly identical magnitudes 
and periods as well as their relative closeness to the centre of the cluster (distances 
of 2'-0, 3':4 and 5’°5 respectively compared with an estimated radius of 47 Tuc 
of 26'-8 (37)) strongly suggests cluster membership. These points were all 
stressed by Shapley (12). We can make membership even more likely by 
computing the expected number of foreground long period variables. At this 
galactic latitude Shapley (38) finds 19 long period variables down to 14™ in 400 
square degrees or 0-048 variables per square degree. Within a radius of 6’ we 
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Taste VIII 


Red variables in 47 Tucanae 
Mag (pg) 
Max. Min. 


113 < 16°0 
11°55 15°3 


11°35 16°1 
130 14°0 
13°3 13°8 
12°7 14°7 1504 + (cycles) 


therefore expect only 0-0015 long variables and furthermore this refers to stars of 
all periods and of all magnitudes down to 14™. Similar considerations apply to 
the irregular variables. 

The radial velocities for the ¢ camera spectra were obtained from measures 
with the projection micrometer using the absorption line wavelengths adopted 
by Merrill (39) in his extensive investigations of Me variables. Also following 
Merrill, laboratory wavelengths were adopted for the hydrogen emission lines. 
In the case of the few d camera spectra these were reduced using Kinman’s 
globular cluster wavelengths (23), as discussed above, for the absorption lines 
and laboratory wavelengths for the emission lines. It may be mentioned that 
on account of their brightness the spectra of the regular variables at maximum 
are of better quality and the velocities derived are somewhat more accurate than 
tHose of the red giants. The velocities adopted are listed in Table II. As in 
the case of the red giants a correction of +10km/sec was applied to the 
d camera velocities and the c and d camera results were given equal weight in 
taking the mean. 

The velocities of the variables 5, 7 and 8 (viz. — 15-0 — 7-7 and — 23:1 km/sec) 
agree with the cluster velocity of — 14-2 km/sec within the limits to be anticipated 
from the velocity dispersion as discussed above. The measured radial velocities 
of the three long period variables are rather puzzling and require some discussion. 
Whilst the mean of the three radial velocities deduced from the absorption lines 
(—7°8, —5:2 and —37-7km/sec) is — 16-9 km/sec and clearly agrees as well 
as can be expected with the mean cluster velocity (— 14:2 km/sec) yet the three 
velocities scatter rather widely about the mean. Examination of the velocities 
in Table II shows that whilst several red giant members have velocities as high 
or higher than variables 1 and 2, the velocity of variable 3 is lower algebraically 
than that of any other star measured in the cluster. We may note in this 
connection that the mean square space velocity in the cluster is 3 x 21 = 63 (km/sec)? 
taking the larger estimate of the velocity dispersion given above. According 
to Chandrasekhar ((40) eq. 5-306, p. 206) the mean square velocity of escape 
is four times this. Thus the root mean square velocity of escape is 
V 63x 4+=16km/sec. The radial velocity of variable 3 differs from that of the 
cluster (— 14:2 km/sec) by 23-5 km/sec which taken at its face value means that 
the variable could be (but is, of course, not necessarily) escaping from the cluster 
(though it is worth noting that Kurth (41) considers the above equation for the 
velocity of escape as “‘ of little significance ’’). 
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The question naturally arises, does this place the membership of variable 3 
and perhaps also of variables 1 and 2 in doubt? The predicted numbers of 
field variables given above is so small that even allowing for a fairly large margin 
of error in the estimates we still find the probability of the stars being foreground 
variables as practically negligible. Further, besides the nearly identical magni- 
tudes and periods of these three stars, we have the fact that so far as can be judged 
from the available data the median visual magnitudes of the variables are about 
the same as those of the stars at the tip of the giant branch. This result would 
probably be expected if the variables are really cluster members. Furthermore 
it is shown below that the variables are normal Me variables. Merrill’s work 
(42) shows that the Me variables of about 200 day period have large peculiar 
motions so that even an approximate agreement of radial velocities for the three 
variables with that of the cluster is unlikely unless they are really members. It 
appears rather safe therefore to adopt the view that these variables are indeed 
cluster members. ‘The rather high velocity dispersion indicated for these three 
stars and in particular the low velocity of variable 3 remains an outstanding 
problem. We have already discussed a possible dependence of velocity on spectral 
type for the red giants and perhaps we may see in the results for the variables an 
extension of this phenomenon. In this case the measured radial velocities may 
reflect mass motions in the atmospheres of the variables connected with the 
evolution of the stars from the top of the giant branch to the horizontal branch, 
an evolution which may be quite rapid (cf. (26)). On the other hand it is possible 
that variable 3 and perhaps variables 1 and 2 also may be members of binary 
systems showing orbital motion and it is intended to test this hypothesis by 
observing these three stars for radial velocity from time to time in the next few 
years. 

We now turn to the spectral types of the variables. An earlier report (43) 
of the three regular variables was based on low dispersion Newtonian spectra. 
We find all six of the variables studied to be of M type with TiO bands somewhat 
stronger than in the latest type non-variable giants and we classify them as M2 
or 3. The irregular variables 5 and 7 do not show any emission lines on our 
spectra. H emission has been found in all three long period variables and in 
the semi-regular variable 8. The long period variables have luminosity classes 
in the range II—Ib whilst the other variables appear to be somewhat fainter 
(11—II1). 

The long period variables have not been found to show any other emission 
lines except H although these latter were often very strong. This is consistent 
with the known behaviour of the Me variables of moderately short period. We 
find no spectroscopic difference between these three stars and normal Me variables 
of comparable periods outside clusters. One slight peculiarity was, however, 
noted on some of the spectra. That is, the ratio Cri 4254/Fe1 4250 indicated a 
slightly earlier type than the Ti O bands (i.e. Cr 1 4254 is too weak for the assigned 
type). This effect was not noticed in any of the non-variable M stars in the 
cluster but Radcliffe spectra of several galactic Me variables also show this weak- 
ening of 4254 and it may possibly be a general phenomenon. Incipient emission 
is one possible cause. In Plate 11(i) and (ii) spectra of variables 1 and 3 are 
reproduced from 49A/mm(Hy) plates. The normal galactic Me variable RV 
Pup of period 1904 is shown for comparison. It is known (44) that in Me 
variables there is a dependence on period of the difference between the absorption 
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and emission velocities. In the case of the three Me variables in 47 Tuc we find 
an average difference A—E=-+5-6km/sec whereas Merrill gives a value of 
+9km/sec, with a considerable scatter for this period. The agreement is 
satisfactory. 

Since most of the spectra of the three Me variables were taken at or near 
maximum light there is little evidence of any variation in the spectra with phase. 
However the spectrum of variable 1 taken on 1953 November 20, when the star 
was judged to be near minimum light shows a weak apparently normal absorption 
spectrum and Hé only in emission, the other H lines being absent. 

A definitive determination of the absolute magnitudes must await accurate 
magnitudes and colours for the variables. However adopting the magnitudes at 
maximum given in Table VIII and a colour index of + 1-7 as is indicated both 
by the preliminary colour magnitude diagram and by what is known of galactic 
stars of this type, we obtain the results of Table IX for our two assumptions as 
to the distance of 47 Tuc. 


TaBLe IX 
Absolute magnitudes of the Me variables (V1, V2, V3) 


Assumed M, of the RR Lyrae stars of. 
Mean apparent photographic magnitude of Me variables at 
maximum r1M-g 
dopted colour index + 1m.7 
dopted visual magnitude at maximum g™-8 
Deduced absolute visual magnitude at maximum —4m-; 
Absolute visual magnitude from the spectrum (luminosity 
class II-Ib) —3zM-5 
Absolute visual magnitude from galactic statistics (45) —2m.2 


Irregular and semi-regular variables (Vs, V7, V8) 
Mean apparent photographic magnitude at maximum 13™-o 
Adopted colour + m-7 
Mean visual magnitude at maximum 11%-3 
Deduced absolute visual magnitude at maximum 
Absolute visual magnitude from spectrum (luminosity 
class II-III) —2M-o 


In the case of the irregular variables it is difficult to compare the deduced 
absolute magnitude with that of galactic irregulars as these have a rather large 
spread in absolute magnitude (cf. (46)). ‘The assigned spectral types place the 
stars between luminosity classes II and III but probably near II. The estimated 
absolute magnitude from the spectrum was obtained from these luminosity classes 
and the calibrations of Keenan and Morgan (18). ‘The agreement is satisfactory 
especially with the value derived from the smaller distance. Table IX also shows 
that in the case of the three Me variables the agreement between the absolute 
magnitude based on the distance of the cluster and that derived from the spectrum 
is best if the smaller distance is used. For Me variables of period about 200 days, 
Merrill and Wilson (45) derive an absolute magnitude of —2™-2 from space 
motions. In view of the uncertainties involved the difference from —3™-6 
cannot be considered as significant. Studies in the Magellanic Clouds will 
probably eventually show whether all 200 day long-period variables are as bright 
as those in 47 Tuc. 
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The three 200-day Me variables in 47 Tuc appear to be the only known Me 
variables in globular clusters although the motions of the variables of this period 
in the general field suggest that these latter stars are related kinematically to 
the globular clusters. A number of variables of periods near 1004 exist in 
globular clusters and Arp (33, p. 124) has suggested that these stars represent the 
short period end of the Me sequence. Whilst it is possible that this is the case, 
the spectroscopic work on these 100% variables in clusters (47) throws some 
doubt on such an interpretation. The observed spectra are quite different from 
normal Me spectra, being considerably earlier in type (early G). It is, of course, 
true that there are few Me variables with periods quite as short as 100 days in 
the general field. However SS Her(1074, M3e) and T Cen (g14, Moe) are two 
Me variables that lie in this range and yet they have much later types than the 
globular cluster variables. 
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Summary 


The distribution of neutral hydrogen away from the galactic plane has 
been measured between Dec. = + 90° and Dec.= —32°. It shows a number 
of neutral hydrogen clouds which coincide in position with dust complexes. 
The major feature of the distribution is the excess of neutral hydrogen in the 
position of the local system (Gould’s Belt) of early-type stars and dust. 
The system (~ 10° MO) appears to represent a later formation in the local 
spiral arm. An extended region in Cepheus is probably similar. 





1. Introduction.—Gould (1) noted that the brighter naked-eye stars of the 
northern and southern Milky Way formed a belt inclined at an angle of 20° to the 
galactic plane. Shapley (2) found this belt of stars to be a sub-system of the 
Galaxy localized within several hundred parsecs of the Sun. The B stars within 


200 parsecs of the Sun have the spatial distribution shown in Fig. 1 with a semi- 
amplitude of 12°. These stars are mainly responsible for the naked-eye Gould 
Belt. Bright Be stars of magnitude less than 6-5 show the same distribution. 
Some constituents of the interstellar medium, notably reflection nebulae which 
are plotted in Fig. 1, also indicate a system inclined at about 20° to the galactic 
plane (3). These facts suggest that there is a galactic sub-system within which 
the Sun lies having the properties summarized in the following table (4, 5) 


Tasie I 


Properties of the Local System 
N. Pole l'=170° b= +74° 
Centre lI’ =240° b= -—3° 
Radius ~ 500 parsecs 
Mass ~10° MO ? 


The reality of this proposed sub-system has been questioned (e.g. 4, 3); 
it may be the observational effect of a dust density which falls at increasing 
distances or it may be the effect of a fortuitous irregular distribution of stars. 

Strong new arguments in favour of the reality of a local system can be derived 
from the neutral hydrogen distribution away from the galactic plane. Lilley and 
Heeschen (6) have noted a correlation between neutral hydrogen gas and dust in 
the southern parts of the sub-system. A preliminary investigation of Jodrell 
Bank hydrogen-line records showed a much more extensive correlation between 
the neutral hydrogen distribution and the local system (7). 
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Reflection Nebulae 
Relative Sizes ~ Hubble 
+ 


i+ 30 














B-Type Stars~- Shepley 





























Fic. 1.—The full-line represents in galactic coordinates the best fitting plane drawn through 
the bright B stars (after Shapley). The positions of large reflection nebulae (after Hubble) are 
also shown. 


2. Programme of observations.—During the course of the hydrogen-line 
programme of the 30ft paraboloid (beamwidth 1°-5 x 1°-7) in 1955 and 1956 a 
large number of drift curves at constant frequency were taken at various decli- 
nations. ‘These employed bandwidths of either 18 Kc/s or 40 Kc/s and showed a 
clear separation from the galactic plane of a broad feature centred about 20° 


away from the plane. 

The preliminary analysis of this data and of Muller and Westerhout’s (8) 
spectra taken near the galactic plane pointed to the desirability of a complete 
spectral survey of the sky visible from Jodrell Bank. This complementary 
survey consisted of setting the 30ft paraboloid at various declination settings 5° 
apart and allowing the rotation of the Earth to sweep the beam through 24 hours 
in right ascension. During each drift the hydrogen comparison radiometer was 
constantly sweeping a frequency band 1-2 Mc/s in width so that spectra were 
obtained at least every 5° in the sky. The observing bandwidth was 100 Kc/s; 
this was adequate to measure the line integral of neutral hydrogen at each point. 
Each declination run was repeated. The survey covered declinations — 32° to 
+go0°. In all, 2500 spectra were taken and analysed. 

The velocity characteristics of the neutral hydrogen in the local sub-system 
well away from the galactic plane were investigated by taking narrow-band 
spectra in sample regions. Bandwidths of 100, 25 and 3°5 Kc/s were used 
simultaneously ; the 25 Kc/s channel gave the spectra required while the 3-5 Kc/s 
channel was used as a check to show that the 25 Kc/s channel produced no 
broadening in the observed profiles. 


3. Analysis and results 


(i) Drift curves at a single frequency.—In certain regions of the sky a drift 
curve across the galactic plane did not show a single broad peak but a broad peak 
centred on the radio galactic plane with a second fainter peak sometimes as much 
as 20° in galactic latitude from the first. Such constant declination drift curves 
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Fic. 2.—Constant frequency drift curves across the galactic plane and the local system at (a) 
Dec.= 40°, (6) Dec.=26°, (c) Dec.=—10°. Galactic coordinates and the frequency difference 
Av from the Ls.r. are given. 
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Fic. 3.—Drift curves across the galactic plane and the local system at Dec.=39° for various 
frequencies relative to the Ls.r. 
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Fic. 4.—The distribution of N (the integrated number of neutral hydrogen atoms per cm* in the line 
of sight) for the northern sky in units of 10® cm~*. 


taken at frequencies near the local standard of rest (l.s.r.) are shown in Fig. 2. 
The ordinates are brightness temperature while the abscissae show, as well as the 
right ascension, the galactic coordinates (/’, 6’) and the changing frequency Av 
(although constant relative to the Earth) relative to the local standard of rest in 
the solar neighbourhood of the Galaxy. Fig. 2(a) is a drift curve across the 
Cygnus region and shows a strong asymmetry south of the plane lying between 
galactic latitudes —5° and —15°. Fig. 2(6) in the Taurus region gives a clear 
separation between the main galactic peak and the neutral hydrogen belt at 
latitudes —15° to —30°. A similar separation is visible above the plane in the 
Ophiuchus—Scorpio region shown in Fig. 2 (c). 

Hydrogen at different distances can be distinguished due to galactic rotation 
by taking drift curves at different frequencies. This criterion still applies in 
the anticentre region where galactic rotation is nearly perpendicular to the line of 
sight. Such a series of drift curves at different frequencies was taken at Dec. = 39° 
and the results are shown in Fig. 3. The frequency relative to the local standard 
of rest at any right ascension for a given curve is obtained by adding the frequency 
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Fic. 4.—Continued. 


correction Av at that right ascension to the frequency assigned to that curve. The 
H1 spectrum in this region of the sky is a single line 160 Kc/s in width at the 
half temperature points. The drift curve for higher frequencies is narrow 
and consequently refers to more distant neutral hydrogen. The same pheno- 
menon is observed in the absorption spectrum of Taurus-A (9) where absorption 
which occurs in nearby hydrogen is observed at the lowest frequencies. Hydrogen 
at negative frequencies (recession) relative to the local standard of rest has one 
broad peak on the galactic plane and a secondary maximum 5° to 10° below the 
plane. This local hydrogen some distance from the plane at zero or negative 
frequencies relative to the l.s.r. is visible over a range of galactic longitudes. It 
occurs above the plane around I’ = 330° and below the plane around I’ = 150°. 
(ii) The survey of neutral hydrogen spectra in the Northern Sky.—This survey, 
taken with a bandwidth of 100 Kc/s was designed to measure N, the in 
number of neutral hydrogen atoms in a 1 cm* column in the line of sight at each 


point. 


N= | Ngdl = 2-3 x 108 | 1(0)dv 
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where N,, is the number of atoms per cm', / is the depth, r(v) is the optical depth 
of neutral hydrogen with a velocity v cm/sec. For this evaluation the kinetic 
temperature of the neutral hydrogen was taken to be 125°K. 


The group of spectra obtained at a given declination was analysed in the 
following way. The integrated quantity {7(v)dv was measured with a plani- 
meter for each spectrum. 7(v) is the observed temperature on an arbitrary scale. 
These quantities are standardized by comparing spectra taken near the plane in 
the survey with similar spectra in Muller and Westerhout’s (8) catalogue of 
profiles. All brightness temperature spectra were then converted into optical 
depth spectra and the value of N was derived for each. The value of N for the 
declination range — 32° to +90° were plotted in galactic coordinates. Smooth 
contours were drawn through these points and the resultant map is shown in 
Figs. 4 (a), (6), (c) and (d). No attempt has been made to represent detail less 
than 5° in diameter. The units of N are 10” neutral hydrogen atoms per cm? 
column in the line of sight. 


This map may be compared with a map prepared by Erickson et al. (10) on 
a 10° grid for latitudes more than 20° from the plane. The agreement with this 
less detailed map, allowing for the difference in grid spacing, is good. 

(iii) Narrow band spectra.—Spectra obtained with half-power bandwidths 
of 25 and 3-5 Kc/s were taken at the positions of enhanced N found away from 
the plane in Fig. 4; regions in the Gould Belt area were included. The central 
fgequency, peak temperature and half-width of these profiles were compared 
with the same properties of profiles obtained nearby on the galactic plane which 
were checked against the Westerhout and Muller profiles. 

Averages of three profiles taken with the 25Kc/s bandwith channel in each of 
eight regions away from the plane are reproduced in Fig. 5. The profiles 
obtained with the 3:5 Kc/s channel were not significantly different from those 
taken with the 25 Kc/s bandwidth. Peak brightness temperatures of up to 
55°K were measured 20° off the plane along both axes of the Gould Belt system. 


A number of the brightness temperature spectra show half widths to half 
temperature points of 7km/sec. If the kinetic temperature is taken as 125°K 
the optical depth spectra indicate a value of halfwidth (y) of 5-8km/sec. This 
may be compared with values of 12 to 10 km/sec for the galactic plane 
profiles at distances of 4 to 8-2 kiloparsecs from the galactic centre (11). 
The latter profiles probably include many neutral hydrogen clouds of the Gould 
Belt type in relative motion. 

Other profiles in the present survey were broader and showed more complex 
structure. Those at /'=330°, b’'=+20° and l'=330°, b’=+40° showed a 
broad base with a narrow bright peak. Some of these regions more than 20° 
away from the plane have values of » as large as 12 km/sec. 


The centre frequency referred to the local standard of rest of the profiles 
plotted in Fig. 5 was determined to an accuracy of +5 Kc/s(+1km/sec) by 
referring spectra taken near the galactic plane to spectra in the catalogue of 
Muller and Westerhout. There is a significant departure of the centre frequency 
of most spectra from the zero of the |.s.r. The spectra representative of the 
Gould Belt region, i.e. near /’ = 330°, b’= +20° and l'=150°, b’=20° show a 
centre frequency at — 15 Kc/s( +3 km/sec) which represents small but significant 
recession of this neutral hydrogen. 
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Fic. 5.—Narrow band (25 Ke/s) profiles taken in regions away from the galactic plane. 


Tase II 


Properties of spectra taken in some regions of large N observed away from 
the galactic plane 

R.A. Peak N WidthofT 7 centre freq. 
(1958) b’ temp ounitsof spectrum, Km/ _ relative to 
h m °K =s.:10"cm-* Ke/s Ls.r., Ke/s 
17 43 93 95 170 —13 
16 40 53 20 70 , —13 
15 30 27 9 7° ‘ + 3 


43 106 80 105 . — § 
28 37 19 127 . —15 
20 55 16 67 : —20 
14 48 28 125 ° —10 
5° 42 16 80 . — 3 


30 18 9 100 r +13 
10 10 x . +10 
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Erickson et al. have plotted the departure of the centre frequency from the 
Ls.r...0f spectra for which |b’|>20°. These show a velocity contraction of 
5 km/sec for |b’|> 40° and at their low latitude limit of b’= + 20° they check the 
present results. Thus there appears to be an inward movement of neutral 
hydrogen near the Sun at high galactic latitudes as distinct from an outward 
movement within the plane of the Gould Belt system. The properties of each 
of the spectra taken in the present 25 Kc/s survey are summarized in Table II. 

4. Features of the spatial contours of N.—The main features of the distribution 
of N, the integrated number of neutral hydrogen atoms per cm? in the line of 
sight, are to be found plotted in Fig. 4 (a), (0), (c) and (d) and may be described 
as follows: 

(i) The greatest values of N lie along the ‘‘radio’’ galactic plane (b” = 0) which 
is inclined at 1°-5 to the Lund plane (6’=0). ‘The maximum value of N on the 
plane is about 150 x 10” atoms per cm? and the minimum is 70 x 10” atoms per cm? 














Fic. 6.—The distribution of N in longitude between latitudes 10° and 60°. The plot for the 
northern hemisphere is a full line above, and that for the southern a broken line below a fiducial 
level for each latitude. 
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assuming a kinetic temperature of 125°K. ‘There are significant amounts 
of neutral hydrogen in both polar regions where values of N in the range 2 to 
4 10” atoms per cm? are found at |b’|> 70°. 

(ii) The main departures from symmetry of N about the galactic plane and about 
the galactic centre are summarized in Fig. 6 which shows the distribution of N with 
galactic longitude plotted at 10° intervals of latitude. The values of N for the 
northern hemisphere are plotted above, and those for the southern hemisphere 
are plotted below a fiducial baseline. The major asymmetries are 
(a) a bulge about the galactid plane in the longitude range 60° to 120° 
extending from b’= +15° to b’=+5° and (6) the local sub-system of Gould’s 
Belt seen above the plane from /’ = 220° to 350° and below the plane from I’ = 110° 
to 190°; it can be traced out to latitudes of 30° and 40° above and below the plane. 

-(iii) Concentrations of neutral hydrogen of limited extent were seen away from 
the galactic plane in a number of regions. Fourteen which were separated from 
the background distribution are listed in Table III. For each is assigned a 
position, extent and a value of AN, the increment of N above the surroundings. 
In addition the last three columns of the table show optical data which will be 
discussed below. 


Taste III 
The properties of 14 regions of enhanced neutral hydrogen emission 


H1 region Extent N units Obscuring Bright Extragalactic 
designation b’ lxb of 10” — clouds nebulae absorption 
cm~* (mags) 
o —20 35X10 Present ? 1-2 
60 +35 10x 7 1-2 
65 +50 15x 8 1-2 
60 + § 
to 120 to+15 60x10 , Avoidance 
70 —28 10x § ° 
70 —§0 I§X § I-2 
110 ° 
190 to—40 80x40 Present Avoidance 
120 —30 I0X1I§ 8 Present 1-2 
140 —40 I§XI5 15 Present C (Pleiades) 2 
155 —55 I0X10 5 ° 
K 175 —20 10X25 18 Present C+E(Orion) Avoidance 
L 220 +5 Present C (Scorpio etc) Avoidance 
350 to +30 130X25 10 
M 320 +45. 30x20 6 Present 1-2 
N 340 —40 30X10 5 o? 

5. The distribution of N in latitude.—A representative picture of the amount 
of neutral hydrogen at high galactic latitudes can be gained from a study of the 
latitude distribution of N. Accordingly the average value of N observed at 
each latitude was plotted against latitude for the northern and southern hemisphere 
in Fig. 7. To avoid any effect of the biasing of these plots by the sum of localized 
regions a second set of plots for each hemisphere was derived from the average of 
the ten lowest values of N (found in any 10° longitude interval) for each latitude. 
These are of similar shape to the former plots but are reduced by about 30 per cent 
in magnitude ; they probably represent the broad distribution of neutral hydrogen 
at high galactic latitudes. 





492 R. D. Davies, A study of neutral hydrogen Vol. 120 


The plots were fitted to a uniform thickness model of the galactic neutral 
hydrogen using a secant-law analysis. The minimal distributions were fitted by 
layers above and below the Sun of N=2-3x10"cm~* and 3:1 x 10%cm~* 
respectively. The average distributions were equivalent to layers above and 








Fic. 7.—The distribution of N in latitude averaged over all values of longitude. Minimal 


values are also plotted. The full line is for the northern hemisphere, the broken line for the 
southern. 


below the Sun of N=3°5 x 10%cm~* and 4:3 x 10%cm-~ respectively. The 
difference between the minimal and average distributions may be ascribed to the 
effect of the non-uniform neutral hydrogen distribution near the Sun already 
discussed. Both distributions suggest that the Sun lies above the neutral 
hydrogen galactic plane with some 20 per cent more hydrogen below the plane 
than above it. If the distribution of neutral hydrogen in the solar neighbour- 
hood is taken as uniform with a thickness of 220 parsecs (12) then the above 20 
per cent difference would be accounted for if the Sun were about 22 parsecs above 
the neutral hydrogen plane. Analysis of stellar distributions indicates that the 
Sun is 13 parsecs above the stellar galactic plane (13). 

The observed amount of neutral hydrogen at high galactic latitudes can be 
compared with that expected from the observed value of N within the galactic 
plane and the known dimensions of the galaxy as mapped in neutral hydrogen. 
For example, the value of N within the galactic plane in the Cygnus direction is 
150 x 10"%cm~* which is integrated over a depth of about 8kpc. The value 
expected at high latitudes with a depth of about 100 parsecs is accordingly about 
2x 10” cm~* if we assume the H1 density in the solar neighbourhood to be that 
of the average in the line of sight in the direction of Cygnus. In view of this 
assumption, the agreement is good with the value of 2—3 x 10” cm~* observed 
for the minimal distributions. 





No. 5, 1960 in the solar neighbourhood of the Milky Way 493 


6. A comparison of the neutral hydrogen distribution with optical features.—The 
correlation between the density of neutral hydrogen and the existence of dust 
clouds has already been established for a number of regions within the galaxy 
(14, 15, 16). The present study has included an examination of types of 
Population I objects which are known to co-exist with neutral hydrogen in the 
spiral arms of the galaxy. 


(i) Obscuring clouds 

A map prepared by Lundmark (17) shows the positions and sizes of those 
dust regions in which the star density is less than 5 or 6 times the surroundings. 
In consequence the higher latitude regions may be in error, being only statistical 
fluctuations in star density. All those dust clouds greater than 1-o degree in 
diameter have been plotted in the maps of N in Fig. 4 (a), (6), (c) and (d). 

A correlation is visible between many features of the distribution of obscuring 
matter and the distribution of neutral hydrogen away from the galactic plane. 
There is an excess of obscuring matter in the plane of the local system (regions 
G and L) in the well-known dust regions of Taurus, Orion and Auriga which lie 
below the galactic anticentre and the regions of Ophiuchus, Scorpio and Scutum 
above the galactic centre. The bulge seen in the neutral hydrogen distribution 
above the plane between longitudes 60° and 120° in Cepheus (region D) also 
has a counterpart in Lundmark’s map; this region is a dust complex about 
500 parsecs distant (4). There also appears to be good evidence for an excess of 
obscuring matter in the small regions A, H, I, K and M. These correlations 
are indicated in column 6 of Table III. 

Hubble’s (18) measurements of the effect of galactic obscuring matter on 
counts of extragalactic nebulae give another indication of the distribution of local 
obscuring matter. The agreement with the general features of Lundmark’s 
catalogue is good and in all cases where there was a correlation between the 21 cm 
measures and Lundmark’s obscuring clouds there was a clear deficiency in the 
counts of extragalactic nebulae. Hubble’s observed number deficiencies have 
been converted into magnitudes of absorption in column 8 of Table III. The 
zone of avoidance of extragalactic nebulae coincides with the shape of the local 
system and also with the Cepheus bulge. 

(ii) Bright nebulae 

Cederblad’s (19) catalogue of continuous reflection and of emission nebulae 
gives the distance as well as the size of clearly distinguishable galactic nebulosities. 
Those nebulae with distances less than 1000 parsecs have been divided among 
the quadrants centred on the axis of the local system (/’=150°, /’=330°) in 
Fig. 8(a) and the perpendicular axis in Fig. 8(6). Their positions above or 
below the plane and their diameters have been plotted. The local system is 
clearly seen separated from the galactic distribution. It lies inclined at 20° to 
the radio galactic plane (6”=0). These clouds of dust and gas extend for 
200 parsecs in the direction of l’ = 330° and 300 parsecs in the direction of I’ = 150°. 
The more extensive bright complexes are noted in column 7 of Table III where 
their spectral character is given (E for gaseous emission and C for continuum 
reflection from dust). These regions are interwoven with the obscuring matter 
plotted on Lundmark’s map and are at the same distance. The variation of N 
with galactic latitude integrated over each quadrant is given as a polar plot by the 
full-line in Fig. 8(a) and (b). The correlation between N and the existence of 
interstellar dust is again marked. 








R. D. Davies, A study of neutral hydrogen Vol, 120 


1» 205° yo0"=15° 
N40 


1's 105" ~ 195° 
N40 











b's ==2$ 
oY coors eee ee : 7 ho 




















Ne 
4 





9 
” 








bd rr} 
Se Bh | ee 9 BQyeo 5 . 











Uaioe’ 

@ <tpee 

© ipwead <s psc 
@ *Spm 


F1a, 8.—The distribution in height (#) above the galactic plane of bright nebulae (after Cederblad) 
averaged for each quadrant of longitude, The corresponding distribution of N is given by the full 
line polar plot. 


(ifi) The distribution of stars 

The spatial distribution of stars in the solar neighbourhood has already been 
discussed in the introduction. ‘The nearby early-type stars lie within the local 
system inclined at about 20° to the galactic plane. A large number of these are 
of spectral type B, to B, which indicates (20) an age of 20 x 10° years for the local 
system. ‘This may be compared with 100-200 x 10° years for the bulk of the 
stars in the solar spiral arm. 

7. A comparison with radio continuum surveys.--Recent publications have 
indicated that there may be some correlation between the distribution of neutral 
hydrogen and the brightness distribution of the continuum radiation from the 
Galaxy. Mills (21) has found 81 Me/s radiation from spiral structure within 
the Galaxy and Blythe (22) has shown that some 38 Mc/s continuum details seen 
in the Taurus region are also visible in the neutral hydrogen contours, More 
evidence relating to this correlation was derived from a study of the present 
neutral hydrogen survey and from continuum surveys which have been made 
with beamwidths of less than about 10° i.e, comparable with the resolution of the 
present survey, ‘The results are plotted in Fig. g(a) and (4) for frequencies of 
48 (22), 200 (23), and 400 (24) Mc/s. ‘The values are the average brightness 
temperatures within + 5° of the latitudes indicated. 

An examination of these plots shows that there is no close correlation between 
the continuum emission at any frequency and the neutral hydrogen distribution 
at galactic latitudes near + 20°. ‘The only feature which produces an effect in 
both line and continuum radiation is the Orion complex of gas and dust. This 
is not unexpected since the continuum radiation comes from an extended region 
showing the Ha emission centred on the Orion Nebula. Moreover the continuum 
spur at /' #0” reaching to high positive latitudes has no counterpart in the neutral 
hydrogen distribution with the resolution of the present survey, 
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Fic. 6.—The distribution in galactic longitude of N and continuum radiation for (a) b= + 1§° 
and (b) b' = aoe", 


8. The mass of gas, dust and stars in the local system.—-The mass of neutral 
hydrogen in the local system (Gould’s Belt) can only be estimated if the 
dimensions of the system can be derived, Since the neutral hydrogen distribution 
is 80 closely correlated in position with the dust complexes of the local system, it 
will be assumed here that the hydrogen is distribute! . depth the same way as 
the dust. ‘Thus the hydrogen will be taken to hay: . uniform density towards 
I’ = 330° out to 200 parsecs (the Ophiuchus and pio clouds) and towards 
I'= 150° out to 300 parsecs (the Taurus, Orion clouds). ‘The equivalent depth 
of the hydrogen perpendicular to the galactic pla: is 30 parsecs as derived from 
the width to half-intensity in latitude of ~15°. ‘he corresponding breadth in 
the direction /' = 60° to I’ « 240° is estimated at 150 parsecs (width ~ 70°) although 
it is difficult to determine because the local system here crosses the galactic plane 
and might just as well have been taken as 500 parsecs i.e, making the region a 
circular disk. ‘The neutral hydrogen in the former case is 10° solar masses (MQ) 
and in the latter 26x10°M©. The hydrogen gas density in this idealized 
representation of the local system is 2°0 atoms per cm in excess of the gas density 
within the local spiral arm, 
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The mass of the dust in the system was calculated using the same dimensions 
and taking the obscuration along the major axis to be 2 magnitudes as read from 
Hubble’s counts of extragalactic nebulae. The line of sight dust mass, m’, was 
derived using the expression (16) m’ = 3 x 10-5 mgr cm? where m is the absorption 
in magnitudes. The dust masses of the two models used in the neutral hydrogen 
calculations were 2°5 x 10° M© and 6x10* M© respectively. The neutral 
hydrogen gas mass is 40 times that of the dust. 

The mass of stars in the local system is not readily calculable. Schmidt (5) 
gives an approximate estimate of 10° M© for a system 500 parsecs in radius. 
The major uncertainty lies in determining the number of low luminosity stars in 
the system. An approach to the problem can be made using the fact that the A 
type stars with apparent magnitudes less than 6 lie in a plane inclined to the 
galactic plane at 4° compared with 16° for the B stars within 200 parsecs (2). 
This indicates that A stars within the local system amount to only about 30 per 
cent of the A star density in the field stars of the local spiral arm. Assuming the 
mass of the local system consists mainly of A and later type stars with a space 
density 30 per cent of that in the local spiral arm, a useful estimate can be made 
from the known density of field stars (25). The stellar mass of the local system 
calculated in this way lies in the range of 10* to 10° solar masses, and would 
appear to be less than the mass of hydrogen in the system. This result may be 
a with the fraction of neutral hydrogen in the spiral arm near the Sun 

thich amounts to 10 per cent (26). 

‘‘ g. Conclusion.—The present study of the distribution of neutral hydrogen 
away from the galactic plane shows a close similarity to the distribution of young 
stars in the local system proposed by Shapley. Neither distribution is smooth 
although the major features of one appear on the other. The reality of this 
system is further strengthened because it consistently exhibits the characteristics 
of a young Population I object—young stars and a high percentage of gas and dust. 
It then would appear to be a ‘‘recent’’ phenomenon in the local spiral arm if 
this data can be interpreted in an evolutionary sense. The ‘‘bulge’’ seen in 
the dust and neutral hydrogen distributions in the direction of Cepheus may be 
a similar type of system. This can only be substantiated if data on the type and 
distance of the stars embedded within it become available. 

Acknowledgments.—The author wishes to thank Dr R. C. Jennison, who 
played a major part in the development and operation of the equipment used in 
the spectral survey described in Section 3(ii). 


Jodrell Bank Experimental Station, 
Lower Withington, 
Cheshire: 
1959 September. 


References 
(1) Gould, B. A., Uranometria Argentina, 335, 1879. 
(2) Shapley, H., Ap. 7., 49, 311, 1919. 
(3) Hubble, E., Ap. 7., 56, 162, 1922. 
(4) Bok, B. J., The Distribution of Stars in Space, Chicago, 1937. 
(5) Schmidt, H., Bonn. Univ. Veroff, No. 35, 1949. 
(6) Lilley, A. E., and Heeschen, D. S., Pub. Nat. Acad. Sci. 40, 1095, 1954. 


(7) Davies, R. D., Paris Symposium on Radio Astronomy, Ed. R. N. Bracewell, Stanford 
. Univ. Press, Stanford, 1959. 





No. 5, 1960 in the solar neighbourhood of the Milky Way 497 


(8) Muller, C. A., and Westerhout, G., B.A.N., 13, 151, 1957. 

(9) Muller, C. A., Ap. F., 125, 830, 1957. 
(ro) Erickson, W. C., Helfer, H. L., and Tatel, H. E., Paris Symposium on Radio 

Astronomy, Ed. R. N. Bracewell, Stanford Univ. Press, Stanford, 1959. 
(rz) Kwee, K. K., Muller, C. A. and Westerhout, G., B.A.N., 12, 211, 1954. 
(12) Schmidt, M., B.A.N., 13, 247, 1957. 
(13) van Tulder, J. J. M., B.A.N., 9, 315, 1942. 
(x4) Lilley, A. E., Ap. 7., 121, 559, 1955. 
(15) Heeschen, D. S., Ap. 7., 121, 569, 1955. 
(16) Davies, R. D., M.N., 116, 443, 1956. 
(17) Lundmark, K., Med. fr. Astron. Obs. Upsala, No. 30, 1927. 
(18) Hubble, E., Ap. 7., 79, 8, 1934. 
(19) Cederblad, S., Med. fr. Astron. Obs. Lund, Ser. II, 119, 1946. 
(20) Sandage, A., Stellar Populations, Vatican Observatory, p. 41, 1958. 
(2x) Mills, B. Y., Paris Symposium on Radio Astronomy, Ed. R. N. Bracewell, Stanford 
Univ. Press, Stanford, 1959. 

(22) Blythe, J. H., M.N., 117, 652, 1957. 
(23) Droge, F., and Priester, W., Z.f. Ap., 40, 236, 1956. 
(24) Seeger, C. L., Westerhout, G., and Conway, R. G. (in preparation). 
(25) Allen, C. W., Astrophysical Quantities, London, 1955. 
(26) Kerr, F. J., and Hindman, J. V., P.A.S.P., 69, 558, 1957. 





CONVERGENCE OF CHANDRASEKHAR’S METHOD FOR 
THE PROBLEM OF DIFFUSE REFLECTION* 


Philip M. Anselone 
(Communicated by Z. Kopal) 


(Received 1959 October 8) 


Summary 


Let I(r, »), J(7) and B(r) denote the intensity, average intensity and 
source function for radiation in a semi-infinite, plane-parallel, isotropically 
scattering atmosphere, with albedo w,S1 and the only external source due 
to an incident parallel beam. Let J,,(7, »), J,,(7) and §,,(7) denote the 
corresponding Chandrasekhar approximations, derived with the aid of either 
the Gauss or double Gauss quadrature formula. It is proved that J,,—J, 
JJ and %,,-+3 uniformly as m--0o. Error bounds are obtained in the 
non-conservative case (m,< 1). 





1. Introduction.—The general problem of the diffuse reflection of a parallel 
beam of radiation by a plane-parallel atmosphere is formulated by Chandrasekhar 
ip, ((1), pp. 20-23). In this paper we are concerned with the important special 
de of a semi-infinite atmosphere which scatters radiation isotropically. Our 
notation is standard. Let IJ(r,u), J(r) and S(r) denote, respectively, the 
intensity, average intensity and source function of diffuse radiation. Let w, 
(o<@,<1) denote the albedo for single scattering, 7F (F>o) the net flux of 
the beam per unit area normal to itself, and —» (j49>0) the direction cosine 
of the beam relative to the outward normal to the atmosphere. The functions /, 
J and S are positive, bounded, and satisfy the equations 


ol 4 
pA) - 1¢7,n)- 87), (1) 
B(7) = we) (7) + fooF e-", (2) 

1 

I)=4[ Urn) di (3) 
I(o,4)=0 for p<o. (4) 
The method of discrete ordinates used by Chandrasekhar in ((1), pp. 80-87), 
involves replacing the integral in (3) by the Gauss quadrature formula sum 
of order 2m, m=1, 2, 3,..., and solving the resulting problem. For each m 


the Chandrasekhar approximations /,,, J,, and 3,, are the unique bounded 
functions which satisfy the equations 


p Slate) =f(s,4)—-G._(7), (5) 


S,,(7) = wel (7) + LargF e~", (6) 
Jn(1)= 1 > Gyula 7: by) (7) 


I,,(0,#)=0 for p<o. (8) 
* Presented to the American Mathematical Society, 1959 September 3. 
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We make a slight generalization. The constants a,,; and y,,; in (7) pertain either 
to the Gauss formula or to the ‘‘double Gauss”’ formula used by Sykes (2). 
In either case, 

O< bin < +++ SPmm Sts) bm, -i= ~Hmir Gm, -i= 4mi > O- (9) 


For a general discussion of numerical integration see Kopal (3), especially 
Chapter VII. 
This is a statement of our principal results. 


Theorem.—As m-> oo, 


In(7 HL (7, 4), Im(t}-I(7), Sy(7)}>S8(7), 
uniformly for oS7< o and —1 Spt. 

In (4), (5) and (6) this writer derived analogous results for the ‘‘ classical 
transfer problem ’’—no incident radiation, pure isotropic scattering, and constant 
net flux. The same general approach that was adopted in those papers is also 
used below, but with important differences and extensions. Most of the analysis 
is expressed in terms of 3 and §,,. Estimates of the deviation of 5,, from 5 are 
obtained in the non-conservative case (ap < 1). 


2. Integral equations for 3 and 3,,.—It follows from equations (1) and 
(4) that 


I(0, x) = B(0); 
Kru)= [ema paces 
0 —- 
lessen I erupt) nro. 
T ad 
Equations (10), (3) and (2) yield 
8(r)= boy | ~ Ey(lr—t)8(¢) dt + berg e-"™, (11) 


where E, is the exponential integral function of order one; more generally, 





EAx)= [. e~4yr—? dy, a6 £21. (12) 
It follows from equations (5) and (8) that 

In(O,#)=5,,(0) ; 

re ee [es eee 





Litule | ¥ dH Bq(t) ell 
Equations (13), (7) and (6) yield 
Bn(7)= ber |” Eqa(lr—#1)8q(t) dt + borg em, 
where E,,, is a special case of : 


Enr(*)= 2 Bini ec 7m (mil, #20,r2 1. 
1 
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3. Neumann series representations of 3 and %,,.—Let S denote the set of 
bounded measurable functions f defined on the interval oS r< «©. With the 
norm of uniform convergence, || f||=sup|/f(7)|, S is a Banach space. Let A 
and A,,, m21, be the integral operators on S such that if feS then Af and A,, f 
are the functions defined for o < r< 0 by 


(Af)(r)=4 i)  Ex(le—afe)ae, (16) 


(An f(r) =4 | : Exus(\7—t]) f(t)dt (17) 


The following statements are easily verified. If feS then AfeS and A,, feS. 
For each feS, |[Af||< |] f]| and |JA,,f[|<|[ fl], so that A and A,, are bounded 
operators. More specifically, ||A{|=||A,,|]=1, where the operator norms 
have the usual definitions, ||A||=sup {||A/]|:||/||=1} and similarly for ||A,,]]. 
If feS and f20 then Af2o and A,,f2o0. For further properties, see (4) and the 
monograph (7) by Hopf, where classes of integral operators which include A 
and A,, are considered. 

The integral equations (11) and (14) are expressed equivalently as 


(I— w,A)3=A, (I—a,A,, )5,, =A, (18) 
where I is the identity operator on S and 


1 h(7)= fuk e-"™, (19) 
Note that heS and 
[Vl = toro (20) 


Suppose that o<wy<1. Then |{wAj|=|}~A,,||=a.<1, so that the 
operators in parentheses in (18) have the unique bounded inverses given by 
the uniformly convergent Neumann operator series, 


(wh) t= ¥ ay"A",  L-arpy) t= F yA" (21) 


n=0 
It follows from a cursory examination of the effect of these inverse operators 
on constant functions and from the Neumann series that 


A — aA) I1= [1 — aA.) = 1/1 — 0). (22) 
These results imply that the equations in (18) have the unique solutions in S 
given by the Neumann series 


B= SF arA"h, 8,= F wr"A,"h, (23) 


n=0 n=0 
at least if o< w,<1, and in this case equations (18), (20) and (22) yield 


Isls = —, Balls =. (24) 


4(1— mp) 


Now let w,=1. Hopf (7) proved that the equation (I-A)3=A has a 
unique bounded solution in the form of the Neumann series. His reasoning 
applies, mutatis mutandis, to the equation (I—A,,)5,,=A. Therefore, for each 
value of the albedo wy, 0< wy <1, the integral equations in (18) have the unique 
bounded solutions given by the Neumann series in (23). It follows from (23) 
and (19) that 5 and 56,,, are positive functions, which was to be expected. For 
a quite different representation of 5,,, see Chandrasekhar (p. 83, formula (98)(1)). 
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4. Convergence of S,, to & in the non-conservative case.—Throughout this 
section assume that 0< wy, <1. 
It follows from (18) that 


8 —G,, = a AS — w,A,,5,, = (A —A,,)5 + wA,, (5 — G,,), 
(I— a,A,,)(5 — 5,,) = ao(A — A,,)5, 
3 —,, =a (I— w,A,,)“*(A — A, )5. 
By (25), (22) and (24), 
IIS-S,.1s = — A—A,,I1 ISH, 


2F 
IS — 3,11 s ——2—] ||A-A, II. 


4(1—a)? 
The definitions of A and A,, imply that 


(Af—AnfMr)=4 ["  CEa(tsl)- Emi(lsD IPC — 2), 
WAF—An SIS SILA { ~ 1Ev(ts1)— Ena( leis, 


A—A, I< | * [Ea(#)~ Ema(s)las (28) 
The inequalities (27) and (28) yield 
2F a 
IS Sulls me | [E,(x) — Ema(x)|dx. (29) 


1—w)* Jo 
We prove now that the sequence of integrals appearing in (28) and 
(29) converges to zero as m tends to infinity. According to equations (24) and 


(26) of (5), 
E,(y)= in E,(x)dx, Eno(y)= | ” Emi(x)dx, y2o, (30) 


and E,(0)=E£,,,(0)=1. Hence, if o<t<7< o then 
i) * [Ei(#)— Ena(2)ide= { i) a | ey } 1E,(*) — Eqa(x)|d 
< 2—E,(t)— Ems(t)+ i) " [E\(*) — Ey (x)|dx + Eg(r) + Equa(1)- 


According to Theorem 1 of (5), E,,;(x)—>£,(x) uniformly for tgx<¢7 and 
En2(y)—> E,(y) uniformly for y 20, as m tends to infinity. Thus, 


lim sup | ~ |E,(x)— Eys(x)|dx S$ 2[1 — Ey(t)+ E(7)] 
ma? © 0 
for o<t<7r<o. Let to and r—> o& to obtain 


i) ” Uh ia)~Rudabie-0 mow, 


By (28), (29) and (31), 
|A-A,,||-+0 a8 m-> oo, 
\|S—35,,|| +o as m—> o. 
Thus, 5,,(7)—> S(7) uniformly for oS r< & as m—> oo. 
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Since the functions £, and E,,,, are known explicitly, (29) furnishes an estimate 
for the deviation of 3,, from 3. However, it is not satisfactory to have to 
evaluate the integrals for each m. ‘These integrals can be estimated by means 
of the standard error formula for the Gauss quadrature formula (cf. Kopal 
[(3), p. 372]). Unfortunately, this procedure yields estimates which tend to 
infinity as m-» «©. Another approach, which involves the Weierstrass approxi- 
mation theorem and utilizes the fact that the Gauss formula of order 2m inte- 
grates correctly polynomials of orders up to 4m — 1, leads to a sequence of estimates 
which does converge to zero as m-»o. However, this sequence of estimates 
is of doubtful value because the convergence to zero is very much slower than 
experience would lead one to predict. Efforts to obtain practical numerical 
estimates are continuing. 

Exchange the pairs (A, 3) and (A,,, 5,,) in the analysis leading to (26) to 
obtain the ‘‘ relative ’’ estimates, 


- Wy 
\|5 B,alI < 1—™, IA All 13,1], (34) 
0 


\|5 : 3,1] = —o [Sul | |E,(x) i Eni(x)|dx. (35) 
I—- Wy 0 
Whenever ||3,,|| is known or available in good approximation, (34) and 
(35) may be used. Since there is reason to think that the first inequality in 
(24), which was used in the derivation of (27), is a poor estimate of ||5||, it is 
likely that the ‘‘ relative ’’ estimates for ||3—5,,|| derived above are better than 
the ‘‘ absolute ’’ estimates, (27) and (29). 
5. Convergence of &,, to & for the conservative case.—In this section assume 
that w,=1. 
The average intensity in the classical transfer problem is a constant multiple 
of ++4q(r), where q— Aq= 4s, 
q=4 ¥ AWE, (36) 
n= 
and o<gqg<1. In ‘Theorem XII of (7), Hopf expressed the source function for 
the problem of diffuse reflection in terms of g. With appropriate notational 
changes, Hopf’s formula (182’) becomes. 


8(r) = P| me + "ent a(tat | 
. E enti M7) _ =| ¢ meeg( tat |. (37) 


Ho Ho Jo 
Let q,, denote the Chandrasekhar approximation to q or order m, which is 
given by formula (46) on page 75 of (1). According to results reported in (5), 
Im ~ And s bE ns» 
din _ 4 > A,,"E,,3) (38) 


nO 
and o<qg,,<1. ‘The same reasoning by which Hopf derived (37) leads to the 
analogous equation, 


8,,(r) =3F [ 4 [- e ‘dt 


) 


E a e-tee 4, Sah?) 2, [ean (tt |. 
Ho Po Jo 
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It was proved by this writer in (6) that ||q—4q,,||->0 as m—> oo. It therefore 
follows from (37) and (39) by routine arguments that 
||3 —G,,||--0 as m-> o. (40) 
This procedure also yields an estimate of ||3—9,,|| in terms of ||g—g,,||. Since 
no estimates of ||q—gq,,|] have been obtained, there is as yet no ‘‘absolute’”’ 
estimate of || —9,,,|| in the conservative case. 
A consequence of (40) is that 
Bin ( a) >3( 2) as m-> OO, (41) 
where 3() and S,,(«) are the (positive) limits of 3(r) and 9,,(7) as r-> o. 


m 


6. Final remarks._We have proved that ||3—5,,||->0 as m-> oo whenever 
O<m,S1. By (2) and (6), 3—53,,=a,(J —/J,,). Therefore, 


\|7 —J,, || +o as m-> ow. (42) 


m 


The expressions in (10) and (13) for / and /,, in terms of 5 and 9,, can be used 
to prove that, as m-> oo. 


I,,(7, #) > I(r, ») uniformly for oS r< oo and —1SpS1. (43) 


Since the proof of (43) is formally identical with that of Theorem 2 of (5), it is 
omitted. The estimates derived for ||3—5,,|| imply corresponding estimates 
for the deviations of J,, and J, from J and J. 

This completes the proof of the theorem stated in Section I. 
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